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FOREWORD 

An  inhouse  study  was  initiated  in  March  1973  to  investigate  static 
and  dynamic  tests  and  associated  data  reduction  and  analysis  computer 
programs  that  could  potentially  assist  the  Reference  Systems  Software 
Group  in  accomplishing  its  mission:  development  and  evaluation  of 
aided-inertial  navigation  software.  Study  emphasis  has  been  in  the 
areas  of  error  modeling  to  support  the  development  of  the  navigation 
Kalman  Filter  and  the  definition  of  tests  and  development  of  data 
reduction  and  analysis  computer  programs  that  have  the  capability  of 
separating  hardware  and  real-time  software  errors  and  isolating  these 
errors  within  their  respective  classes. 

The  I MU  calibration  technique  reported  herein  is  the  first  in  a 
series  of  six  tests  defined  in  the  study.  Development  of  the  calibration 
methodology,  computer  simulation,  and  a data  reduction  and  analysis 
computer  program  was  accomplished  in  the  June  to  December  1973  time  period. 
Development  of  the  real-time  test  software,  integration  of  the  test 
equipment,  and  experimental  verification  of  the  calibration  technique 
were  accomplished  during  the  September  to  December  1974  time  period. 

Several  refinements  to  the  test  software,  equipment,  and  data  analysis 
computer  program  were  made  in  1975. 

The  author  acknowledges  the  assistance  of  Mr,  D.  Kaiser  who  was  the 
Task  Engineer  responsible  for  the  integration  of  the  laboratory  test 
equipment  and  development  of  the  real-time  test  software.  Capt  P.  Howe 
and  Lt  H.  Vassar  also  contributed  to  the  development  of  the  real-time 
test  software.  Mr.  R.  Clark  refined  and  assembled  the  equipment 
associated  with  the  I MU  gimbal  angle  measurements.  TSgt  E.  Edgecoinb 
conducted  many  of  the  calibration  tests. 
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SECTION  1 
INTRODUCTION 


1.1  BACKGROUND  AND  NEED 

As  part  of  the  Air  Force  Avionics  Laboratory's  Reconnaissance  and 
Weapon  Delivery  Division,  the  Reference  Systems  Branch  mission  is 
concerned  with  conducting  exploratory  and  advanced  development  of 
navigational  reference  subsystems,  analyzing  new  concepts  for  reference 
subsystem  data  processing  and  software  development,  and  developing  new 
techniques  for  subsystem  integration.  The  Branch's  Reference  Systems 
Software  and  Evaluation  Groups  support  its  mission  by  conducting  in-house 
analysis  and  evaluation  of  navigation  reference  subsystems,  with  emphasis 
in  the  areas  of  inertial  navigation  software  and  subsystem  integration 
via  Kalman  Filtering.  The  Evaluation  Group  maintains  a laboratory  and 
instrumented  Greyhound  bus  for  both  static  and  dynamic  hardware/scf tware 
evaluation. 

In  accomplishing  this  mission,  the  Software  and  Evaluation  Groups 
have  and  continue  to  encounter  numerous  mundane  problems  associated 
with  the  design,  development,  test,  and  evaluation  of  aided-inertial 
navigation  software.  Some  of  the  more  notable  problems  are:  How  is  an 
inertial  measurement  unit  ( I MU)  error  model,  suitable  for  Kalman  Filter 
des.gn  and  comprehensive  subsystem  convariance  error  analysis,  developed 
and/or  verified?  Considering  that  inertial  navigators  are  among  the 
most  precise  of  all  el ectro-mecnanical  devices,  how  is  it  readily 
determined  that  the  IMU  is  operating  within  its  specification  prior  to 
conducting  software  evaluation  tests?  At  the  beginning  of  a reference 
subsystem  test,  what  is  the  IMU  platform  alignment  accuracy  considering 
the  wind  buffet  and  personnel  loading  perturbating  motion  effects  on  the 
greyhound  bus?  How  are  the  potential  sources  of  aided-inertial  navigation 
subsystem  errors:  hardware,  alignment  and  software  separated  at  the 
conclusion  of  an  evaluation  test?  How  can  the  various  classes  of 
software  errors,  coding,  equations,  inexact  algorithms,  etc.,  be 
distinguished  so  that  software  development  projects  can  progress  in 
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an  orderly  fashion?  What  is  the  best  approach  for  establishing  a position/ 
velocity/attitude  reference  for  the  aided- inertial  navigation  subsystem 
under  evaluation:  by  an  on-line  real-time  or  by  an  off-line  after- 
the-fact  data  reduction/analysis  computer  program  approach  using 
recorded  data? 

1.2  OBJECTIVES  OF  THE  STUDY 

From  the  preceding  discussion,  it  is  clear  that  aided- inertial 
navigation  software  projects  can  require  extensive  analyses,  test 
procedures,  arid  supporting  data  reduction  and  analyses  computer  programs 
if  comprehensive  development  and  evaluation  efforts  are  desired. 

Recognizing  this  need,  the  Reference  Systems  Software  Group  undertook 
the  development  of  a series  of  tests  and  supporting  data  reduction/ 
analyses  computer  programs  to  facilitate  its  software  development  and 
evaluation  projects. 

Table  1-1  indicates  the  types  of  tests  that  are  currently  planned 
and  are  under  development.  All  tests  have  several  major  points  in 
common.  When  a test  is  conducted,  test  data  is  stored  on  a magnetic 
tape  for  subsequent  processing  by  a large  off-line  scientific  computer. 

The  concept  of  off-line  non-real  time  data  reduction/analyses  computer 
programs  introduces  a time  delay  in  obtaining  test  results  but  affords 
computer  programs  that  are  more  comprehensive  and  easier  to  develop  and 
change  than  those  resultant  from  on-line  real-time  approaches. 

The  performance  analysis,  error  diagnostic  and  evaluation  concept 
for  each  data  reduction/analysis  program  is  a large  dimensional 
Kalman  Filter  operating  in  parallel  with  deterministic  logic  and 
navigation  software  programs.  This  is  a state-of-the-art  approach 
that  ensures  the  optimum  use  of  the  test  data,  minimum  test  time,  and 
utmost  error  estimation  accuracy.  The  extent  to  wnich  Kalman  Filters 
can  be  used  for  error  source  detection  and  isolation  is  a major  finding 
to  be  determined  by  the  study. 
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TABLE  1-1 

AIDED-INERTIAL  NAVIGATION  SOFTWARE  DEVELOPMENT  AND 
EVALUATION  TESTS 

Kalman  Filter 

lest  States  Test  Location 


(1) 

IM’J  Calibration 

68 

SEL 

(2) 

MEL  Suspension 

35 

MEL 

(3) 

Scorsby  Table  Platform  Alignment 

39 

SEL 

(4) 

Tilt  Table  Static  Navigation 

37+ 

SEL 

(5) 

MEL  Perturbed  Base  Alignment 

39+ 

MFL 

(b) 

MEL  Dynamic  Navigation 

43+ 

MEL 

SEL  - Software  Evaluation  Laboratory 

MEL  - Mobile  Evaluation  Laboratory 
+ - Final  Filter  State  Size  Will  Be  Somewhat  Larger 

The  major  objectives  of  each  of  the  tests  are  as  follows.  The  IMU 
calibration  test  addresses  the  problems  associated  with;  IMU  error  model 
development  and  validation,  calibration  of  random  bias  type  IMU 
hardware  errors,  estimation  of  time-varving  error  sources,  and  IMU 
hardware  status  assessment  prior  to  software  developmental  testing. 

The  MEL  Suspension  Test  utilizes  three  orthogonal  strapdown 
accel erometers  mounted  in  the  greyhound  bus  (MEL)  to  obtain  data  relevant 
to  the  perturbed  base  platform  alignment  problem:  MEL  suspension 
physical  characterist ics  including  mass,  spring  and  shock  absorber  data, 
IMU  relative  mounting  data,  manloading-engine-wind  perturbating  motion 
magnitudes,  and  statistical  modeling  of  these  disturbances. 

The  Scorsby  Rafe  Table  Platform  Alignment  Test  is  designed  to 
evaluate  Kalman  filter  alignment  software  in  an  environment  where  the 
frequency  and  magnitude  of  base  motion  perturbations  are  accurately 
contr  ol  led. 
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The  Tilt  Table  Static  Navigation  Test  addresses  initial  Kalman 
Fi 1 ter/inertial  navigation  software  and  IMU/computer  hardware  integration 
problems  in  a highly  controlled  environment.  Emphasis  is  placed  on  the 
ability  of  the  data  reduction/analysi s computer  program  to  detect  and 
isolate  the  potential  errors.  Special  tests  such  as  determining 
navigation  performance  as  a function  of  gyro  heading  sensitivity  can  be 
conducted. 

Ttie  MEL  Perturbed  Base  Alignment  Test  assesses  the  ability  of  the 
alignment  software  to  both  align  the  platform  axes  and  estimate  gyro 
drifts  for  multi-position  platform  alignment  schemes  in  the  presence  of 
real  world  MEL  base  perturbations  and  inertial  instrument  temperature 
transients. 

The  MEL  Dynamic  Navigation  Data  Reduction/Analysis  Program  is 
largely  a composite  of  the  previous  programs  and  utilizes  much  of  the 
data  generated  by  them.  The  overall  objective  is  to  assess  aided-inertial 
navigation  subsystem  performance  and  separate  the  associated  errors  into 
their  respective  classes:  hardware,  alignment  and  software.  Hardware 
and  alignment  errors  are  separated  by  the  inherent  operation  of  the 
Kalman  Filter.  Software  errors  are  further  separated  into:  coding, 
equation,  inexact  algorithm,  inertial  navigation  and  Kalman  Filtering 
by  "ideal  inertial  navigation  equations"  operating  with  logic  in 
conjunction  with  the  large  dimensional  Kalman  Filter  of  the  data 
reducticn/analysi s program.  1 

Each  data  reduction/analysis  computer  program  has  been  developed 
and  exercised  against  a computer  simulation  of  the  appropriate  test 
subsystem  in  its  expected  environment.  Simulation  results  have  been 
encouraging  especially  in  the  area  of  error  detection  and  isolation. 

The  JMU  calibration  data  reduction/analysi s computer  program  has 
additionally  been  exercised  against  real  world  I MU  test  data  as  collected 
in  the  Evaluation  Group's  Software  Evaluation  Laboratory  (SEL).  The 
experimental  results  obtained  clearly  indicate  the  value  of  the  Kalman 
Filter  as  an  effective  tool  for  detecting,  isolating,  and  accurately 
estimating  the  error  sources  in  complex  navigation  subsystems. 
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1.3  REPORT  OVERVIEW 

This  report  is  concerned  with  the  design,  development,  analysis, 
and  laboratory  test  results  of  the  I MU  calibration  technique.  The 
remaining  five  data  reduction/analysis  computer  programs  in  the  aided- 
inertial  navigation  software  development  and  evaluation  test  series  are 
to  be  reported  after  they  are  exercised  against  read-world  data. 

In  Section  2,  the  ground  rules  and  versatility  and  accuracy  goals 
for  the  I MU  calibration  technique  are  discussed.  The  IMU  calibration 
methodology  formulated  in  response  to  the  restrictions  and  goals  is 
presented. 

In  Section  3,  an  error  model  for  the  IMU  is  given  and  the  error 
equations  governing  its  operation  during  the  calibration  test  are  derived. 

Section  4 contains  a discussion  of  the  conceptual  approach  used  for 
separation  of  the  IMU  error  sources. 

Section  5 contains  a description  of  the  Kalman  Filter  mechanized  for 
estimation  and  calibration  of  the  IMU  error  sources. 

Section  6 presents  and  discusses  theoretical  results  obtained  from  a 
CDC-6600  computer  simulation  of  the  IMU,  Kalman  Filter,  and  calibration 
test. 

Section  7 discusses  the  laboratory  hardware,  computer  software, 
and  the  test  procedure  used  in  testing  the  IMU. 

Section  8 contains  the  "lessons  learned"  from  the  initial  IMU 
calibration  tests  and  the  subsequent  IMU  error  model  refinements. 

Section  9 presents  a pictorial  tabulation  of  the  results  obtained 
from  the  IMU  calibration  tests.  IMU  error  source  calibration  accuracy, 
as  obtained  from  the  Kalman  Filter's  covariance  matrix,  is  given  in 
addition  to  the  recovered  error  source  magnitudes. 
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Section  10  points  out  several  applications  for  which 
I MU  calibration  technique  appears  to  be  well  suited. 

Appendix  A contains  the  CDC-6600  computer  FORTRAN  IV 
to  assist  those  attempting  to  develop  a Kalman  Filter  and 
algorithms  for  a similar  application. 


the  developed 

program  listing 
associated 
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SECTION  II 

I MU  CALIBRATION  METHODOLOGY 


2.1  GROUND  RULES 

The  I MU  used  by  the  Reference  Systems  Software  and  Evaluation 
Groups  in  conducting  aided-inertial  navigation  software  development  and 
evaluation  tests  is  a Si nger-Kearfott  KT-73  four-gimbal , all-attitude 
unit  which  contains  gyroscope  and  accelerometer  inertial  instruments  and 
associated  electronics.  This  IMU  is  also  a part  of  the  A-7D  aircraft 
avionics  suite  and  is  calibrated  at  the  operational  bases  by  the  Singer- 
Kearfott  ASM-375  testset.  The  IMU  is  maintained  by  the  Aerospace 
Guidance  and  Metrology  Center  (AGMC),  Heath  AFB,  Newark,  Ohio.  AGMC 
employs  numerous  digital  IMU  (DIMLl)  testsets  for  calibrating  the  units. 

Initially,  serious  consideration  was  given  to  periodically 
calibrating  the  IMU  by  one  of  the  above  testsets.  This  approach  was 
subsequently  discarded  for  several  reasons.  Funding  was  not  available 
nor  was  it  considered  cost-effective,  procurement  costs  of  the  ASM-375 
testset  were  roughly  one  million  dollars  each,  to  buy  and  maintain  a 
testset  for  the  periodic  calibration  of  one  IMU.  Although  acceptable, 
turnaround  times  for  having  the  IMU  calibrated  at  Heath  AFB  were 
undesirable.  Finally,  the  sixteen  parameters  calibrated  by  the  testsets 
were  considered  insufficient  for  comprehensive  IMU  error  modelling 
studies . 

The  ground  rules  established  for  developing  the  calibration  technique 
were  necessitated  by  the  manpower  skills  and  equipment  available  in  the 
Evaluation  Group  and  its  Software  Evaluation  Laboratory.  The  available 
skills  and  equipment  precluded  disassembling  the  IMU  for  special  testing. 
"No  disassembly"  implies  that  only  easily  accessible  system-level  IMU 
quantized  velocity  pulses  and  gimbal  synchro  angle  output  measurements 
could  be  used.  No  specialized  AGE  was  to  be  required  other  than  that 
already  available  in  the  laboratory:  digital  computers,  tilt  and 
scorsby  tables,  a Si ngar-Kearfott  ASN/90(V)  iner+ial  navigation  system 
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test  equipment  rack,  etc.  A final  ground  rule  was  t.nat  the  solution  to 
the  calibration  problem  be  practically  i nip  1 etrien table  on  a general 
purpose  digital  computer. 

2.2  PtRFORMANCE  AND  VERSATILITY  GCALS 

A major  limitation  of  contemporary  I MU  calibration  concepts  is  their 
limited  diagnostic  capabilities.  Present  test  concepts  assume  a priori 
the  major  sources  of  I MU  errors,  usually  ten  to  twenty,  with  the  test 
being  designed  to  sequentially  determine  the  magnitudes  of  the  assumed 
error  sources.  A deficiency  of  this  approach  is  that  the  test  results 
cannot  provide  any  indication  as  to  whether  or  not  the  test  article 
conforms  to  the  assumed  error  model.  It  is  entirely  possible  for  an  IMIJ 
to  contain  significant  error  sources  other  than  those  assumed  and  for 
them  to  go  undetected  at  the  conclusion  of  the  calibration  test.  Adequate 
Kalman  Filter  design  and  comprehensive  navigation  subsystem  covariance 
error  analysis  are  dependent  upon  having  an  accurate  IMU  error  model. 

One  goal  of  the  study  is  therefore  to  devise  an  IMU  calibration  concept 
that  affords  a capability  for  IMU  error  model  development  and/or 
validation.  Simply  stated,  the  calibration  test  results  should  provide 
some  indication  of  the  degree  that  the  test  article  conforms  to  the 
assumed  error  model . 

It  is  desirable  that  the  calibration  concept  be  versatile  in  regards 
to  both  IMU  mechanizations  and  data  processing  approaches.  The  present 
study  is  largely  concerned  with  the  gimballed  KT-73  IMU  although  future 
Reference  Systems  Branch  projeccs  may  involve  a strapdown  IMU.  A 
versatile  methodology  will  permit  adaption  to  the  strapdown  mechanization 
with  a minimum  of  modifications  to  the  test  procedures  and  data  reduction/ 
analysis  computer  programs. 

Off-line  after-the-fact  and  on-line  real-time  data  processing 
approaches  are  both  desirable.  An  off-line  approach  to  IMU  calibration 
greatly  facilitates  the  development  of  the  Data  Reduction/Analysis 
Computer  Program  since  it  can  be  written  in  a higher  order  language  and  is 
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not  subject  to  the  numerous  real-time  software  development  problems. 
Off-line  data  reduction/analysis  programs  hosted  on  large  scientific 
computers  can  also  be  more  sophisticated  than  on-line  versions  due  to 
the  large  computer  resources  available.  On-line  real-time  approaches 
offer  the  advantages  of  real-time  I MU  closed  loop  control  and  provide 
test  results  immediately  upon  conclusion  of  the  test. 

I MU  calibration  performance  goals  are  to  provide  a comprehensive 
calibration  with  optimum  accuracy  and  to  minimize  I MU  test  time.  Since 
static  IMU  hardware  error  sources  can  be  compensated  for  in  the  operational 
flight  computer  program,  it  is  desirable  to  estimate  each  significant  IMU 
error  source  to  an  accuracy  such  that  the  navigational  accuracy  acnievable 
with  the  IMU  would  be  limited  only  by  uncompensatable  random  disturbances. 
Total  test  time  of  a few  hours  or  less  is  desirable  to  allow  aided-inertial 
navigation  software  tests  to  be  conducted  on  the  same  days  as  the 
cal ibration  tests . 

2.3  FORMULATED  CALIBRATION  METHODOLOGY 

The  performance  and  versatility  goals  dictate  that  a generalized  IMU 
calibration  methodology  be  formulated  which  is  based  only  on  the 
essential  functions  an  IMU  performs.  These  functions  are  to  permit 
the  establishment  and  maintenance  of  a reference  platform  coordinate 
frame  and  measurement  of  the  instantaneous  specific  force  vector.  Noting 
that  calibration  is  an  error  estimation  process,  the  IMU  error  analysis 
methodology  formulated  can  be  understood  with  the  following  rationale. 

If  the  orientation  of  the  platform  of  an  error-free  IMU  were  initially 
precisely  known  with  respect  to  the  physical  earth  rate  and  gravity 
vectors,  then  both  the  orientation  of  the  platform  with  respect  to  these 
vectors  and  the  platform  coordi uatized  specific  force  vector  could  be 
deterministically  calculated  as  a function  of  time.  Any  differences 
between  the  ideal  time-dependent  platform  orientation  and  accelerometer 
outputs  and  those  of  a real  IMU  must  be  caused  by  either  IMU  error 
sources  or  an  initial  platform  misalignment.  This  fact  is  true 
regardless  of  initial  platform  orientation  and  can  be  employed  for  the 


9 


AFAL-TR-77-75 


derivation  of  system-level  I MU  error  equations  relating  measurable 
IMU  performance  parameters  to  1 MU  and  platform  misalignment  angle 
error  quantities. 

Application  of  the  generalized  metnodology  to  a fixed-base  gitnballed 
IMU  can  be  illustrated  by  the  IMU  line  schematic  of  Figure  2-1.  The 
fundamental  function  of  the  gyros,  gimbals,  and  servo  feedback  control 
system  is  to  maintain  a reference  platform  coordinate  frame  relative 
to  an  inertial  frame.  If  the  gyro  triad  of  an  ideal  IMU  were  commanded 
to  precess  in  equality  with  the  earth's  angular  rotation  vector,  then 
the  angular  velocity  of  the  platform  with  respect  to  an  earth  fixed 
reference  frame  would  be  identically  zero.  This  ideal  IMU  performance 
would  result  in  a fixed  angular  orientation  of  the  gimbals  relative  to 
each  other  and  constant  system-level  synchro  gimbal  angle  and 
accelerometer  measurement  outputs.  The  exact  angular  orientation  of  the 
platform  frame  with  respect  to  the  earth  rate  and  gravity  vectors  must 
be  known  for  a priori  computation  of  the  required  gyro  precession  rate 
torquing  magnitudes.  In  the  real  world,  inexact  initial  platform 
attitude  information  and  IMU  error  sources  preclude  the  above  ideal 
IMU  performance  and  synchro  gimbal  angle  and  accelerometer  measurements 
exhibit  a time  dependency. 
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SECTION  III 

ERROR  EQUATION  DEVELOPMENT 


3.1  INTRODUCTION 

The  essential  task  involved  ,n  developing  the  error  equations  is 
the  derivation  of  analytical  expressions  relating  system-level  I MU 
gimUal  synchro  angle  and  velocity  measurements  to  IMU  and  platform-to- 
navigation  frame  misalignment  angle  error  quantities.  These  derivations 
are  relatively  straightforward  and  follow  directly  from  the  physical 
vector  equations  expressing  the  physical  rela tionshi ps  implicitly  implied 
by  the  generalized  methodology. 

A prerequisite  for  these  derivations  is  an  IMU  model  explicitly 
containing  each  IMU  error  source's  functional  dependency.  As  stated  in 
Section  2.1,  a Si nger-Kearfott  KT-73  IMU  was  chosen  for  the  calibration 
methodology  application  largely  because  it  is  used  in  the  MEL  as  part  of 
the  Reference  Navigation  System  and  requires  periodic  status  validation 
and  calibration.  Reference  1 contains  an  error  model  for  the  IMU  which 
reportedly  contains  all  pertinent,  error  sources. 

To  achieve  the  stated  performance  goals,  the  "finalized  IMU  error 
model"  must  contain  all  IMU  error  sources  that  are  significant  relative 
to  the  gyro  and  accelerometer  inertial  instrument  random  disturbances. 

Trie  phrase  "finalized  IMU  error  model"  is  used  to  emphasize  that  the  error 
model  given  in  Reference  1 need  not  be  strictly  accurate  in  the  sense 
that  it  contains  all  error  sources  that  are  significant  relative  to  the 
inertial  instrument  random  disturbances . The  error  model  of  Reference  1 
needs  only  to  represent  a good  starting  point  since  a major  objective 
of  this  study  is  to  develop  a conceptual  approach  to  IMU  error  model 
development  and  validation. 

The  conceptual  approach  to  IMU  error  model  development  and  validation 
that  will  be  pursued  in  the  study  contains  two  basic  premises:  (1)  Esti- 
mation of  the  IMU  error  sources  in  a simultaneous  or  parallel  fashion 
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rather  than  the  sequential  or  serial  manner  of  contemporary  techniques 
will  permit  di screpancies  between  the  assumed  and  actual  error  model  to  be 
detected.  And,  (2)  by  adding  error  model  residuals  or  additional  error 
sources  possessing  very  broad  functional  dependencies  to  the  assumed 
error  model,  it  should  be  possible  to  determine  specific  functional 
dependencies  through  correlation  analysis.  The  method  of  implementing 
this  conceptual  approach  will  be  described  in  the  sequel. 


3.2  GIMBAL  SYNCHRO  ANGLE  MEASUREMENT  EQUATIONS 

The  gimbal  synchro  angle  measurement  equations  express  the  time 
rates  of  change  of  the  IMU  gimbals  relative  to  each  other.  The  key 
point  that  must  be  recognized  in  order  to  derive  these  expressions  is 
that  the  gimbals  rotate  in  such  a manner  as  to  isolate  the  IMU's 
platform  angular  velocity  from  that  of  its  case.  Before  these  relation 
ships  can  be  determined  it  is  however  necessary  to  define:  various 
reference  frames,  the  relationship  between  their  angular  velocities, 
gyro  torquing  rates,  and  a gyro  drift  model. 


3.2.1  Reference  Frames  and  Their  Angular  Velocities 

The  platform  frame  is  defined  as  the  innermost  IMU  gimbal  upon  which 
the  gyro  and  accelerometer  inertial  instruments  are  mounted.  The  inertial 
frame  is  considered  as  a frame  that  is  nonrotating  relative  to  the  stars 
and  has  its  origin  at  the  earth's  center.  The  angular  velocity  vector  of 
the  platform  frame  relative  to  the  inertial  frame  coordinatized  in 
platform  frame  coordinates  is  given  as: 


P q 

c«  v. 


where 

p . 

Cg  is  a small-angle  nonorthogonal  error  transformation  matrix  that 
accounts  for  the  physical  misalignments  of  the  gyro  input  axes 
from  the  platform  reference  frame  axes.  The  elements  of  this 
matrix  are  shown  in  Figure  3-1. 
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w|  is  the  angular  velocity  vector  of  the  gyro  frame  relative  to  the 
inertial  frame,  coord inati zed  in  gyro  frame  coordinates. 


Wj  is  composed  of  two  basic  vectors.  The  first  represents  the 
angular  velocity  at  which  the  gyros  are  commanded  to  precess  relative 
to  inertial  space  WT  The  second  represents  gyro  frame  angular  velocity 
relative  to  inertial  space  due  to  all  gyro  imperfections  and  is  commonly 
referred  to  as  gyro  drift  Wj 


+ 


To  maintain  the  gyro  frame  stabilized  at  approximately  a constant 
attitude  with  respect  to  the  navigation  frame,  it  is  necessary  to  torque 
the  gyros  corresponding  to  the  best  a priori  information  of  the  coordini- 
zation  of  the  earth  rate  vector  in  the  gyro  frame. 


,P  PPG 

'PG  kS 


RG  — N 
C W 
CASE  I,N 


where 

[Kp]  is  a 3x3  diagonal  matrix  of  the  loaded  gyro  scale  factors.  The 
true  I MU  gyro  scale  factors  [Kpy]  differ  from  the  scale 
factors  loaded  into  the  digital  navigation  computer  [Kp]  by 
an  error  designated  as  [AKp], 

PRF  is  a three-state  vector  of  the  gyro  torquing  pulse  repetition 
frequencies. 

wj  N is  the  angular  velocity  vector  of  the  navigation  frame  relative 
to  inertial  space  coordinated  in  the  navigation  frame. 
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'l  is  the  magnitude  of  the  earth  rate  vector. 

L is  the  latitude. 

C(-)  and  $(•)  represent  the  trigonometric  cosine  and  sine  functions. 
RG 

Ccase  is  the  IMU  case  to  roll  (outer)  gimbal  frame  coordinate 
transformation  matrix. 
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Ipg  is  the  roll  gimbal  to  pitch  (middle)  gimbal  frame  coordinate 
transformation  matrix. 
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is  the  pitch  gimbal  to  azimuth  (inner)  gimbal  or  platform 
frame  coordinate  transformation  matrix. 


PG 


eft) 

-»(*.} 


sft) 

eft) 


L o 


o 


G 

0 

3 


and  6z  are  the  nominal  gimbal  synchro  angle  readings  which 
are  obtained  by  subsequent  rotations  about  the  roll,  pitch, 
and  azimuth  gimbal  axes. 
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The  K.T-73  I MU  is  a four-gimba'l  unit  with  an  outer  to  inner  gimbal 
progression  of  outer  roll,  pitch,  inner  roll,  and  azimuth.  The  inner 
roll  gimbal  is  riot  instrumented  with  a synchro  device  and  is  slaved  to 
coincide  with  the  outer  roll  gimbal  by  a reluctance  E-bridge.  Tne 
inner  roll  gimbal  rotates  out  of  coincidence  with  the  outer  roll  gimbal 
to  prevent  gimbal-lock  when  the  pitch  gimbal  exceeds  t70  degrees  from 
its  neutral  position.  The  inner  roll  gimbal  is  not  considered  in  the 
above  gimbal  frame  coordinate  transformations  because  its  angular 
orientation  cannot  be  measured  at  system-level.  The  only  penalty 
associated  with  this  approach  is  that  the  I MU  cannoc  be  tested  with  pitch 
gimbal  more  than  ±70  degrees  from  its  neutral  position. 

3.2.2  Gyro  Drift  Model 

The  gyro  drift  vector  is  an  error  quantity  composed  of  all  gyro 
error  sources  that  are  significant  relative  to  and  including  the  random 
disturbances  or  gyro  noise. 

Ko  - *A  + "b  + * *0  + *T  + *RL.  + *RS 

W.  is  the  platform-to-IMU  case  attitude  dependent  gyro  bias  drift 

vector.  This  erro"  source  is  not  listed  in  Reference  1 by  the 
IMU  manufacturer,  but  has,  however,  been  added,  as  discussed 
in  Section  3.1,  to  account  for  the  possibility  of  a significant 
discrepancy  between  the  Reference  1 and  the  true  gyro  error 
model.  Note  that  the  error  source  has  no  specific  functional 
dependency  other  than  attitude  relative  to  the  IMU  case. 

Wp  is  the  constant  gyro  bias  drift  vector, 
is  the  gyro  nsass  unbalance  drift  vector. 
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Sx  v 
0 - M< 


MSx’  MSy’  an<^  MSz  are  t^ie  9-yr0  sp^n  axes  mass  u 1 ance  coefficients. 

M,  ,MT  , and  Mr  are  the  gyro  input  axes  mass  unbalance  coefficients. 
Gyro  input,  output  and  spin  axes,  and  mass  unbalance 
orientations  used  for  deriving  the  mass  unbalance  matrix 
shown  above  are  given  in  Figure  3-2.  The  vector  equation 
equating  the  cross  products  of  the  mass  unbalance  position 
and  specific  force  vectors  to  the  gyro  angular  drift  and 
angular  momentum  vectors  establishes  the  sign  conve.tions. 

f?  is  the  specific  force  coordinatized  in  the  gyro  frame. 


..  r8  rP  -PG-RG  -.CASE  jN 
t,  ’ °PGCRG  °CASE  N 


CASE 

is  an  orthogonal  small  angle  navigation  to  IMU  case  frame 
coordinate  transformation  matrix. 


„CASE 


I T7  -ft 

'i  'y 

1 -’i,  1 \ 

v,  -v,  i 


is  the  specific  force  coordinated  in  the  navigation  frame. 


rN 

*5=0 
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g is  the  magnitude  of  the  earth's  gravity  vector. 

W0  is  the  gyro  anisoelastic  drift  vector. 

o,  • 'l,  • 

9 

Oy  ' fS*  * 

0,  * fS*  * 

KqX>  and  KQZ  are  the  anisoelasticity  coefficients  for  the 
x,  y,  and  z gyros. 

is  the  gyro  thermal  transient  drift  vector. 

-('-To) 

wr  * wr0  • 

WT  is  the  initial  gyro  warm-up  drift  vector. 

'o 

Tq  is  the  time  from  power  turn-on  to  the  start  of  fine  gyro- 
compass alignment. 

t is  the  time  from  power  turn-on. 

x t is  the  gyro  warm-up  drift  time  constant. 

WR^  is  the  long  correlation  time  exponentially  correlated  gyro 
random  drift  vector. 

WR~  is  the  short  correlation  time  exponentially  correlated  gyro 
random  drift  vector. 

- (0  rj  _ 

"r"1 t,  + °g  V-^-N(t) 
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i „ and  'i  are  the  long  and  short  gyro  random  noise  correlation 
g!L  gs 

times . 

a . and  a are  the  lonq  and  short  correlation  time  random  drift 
g i gs 

standard  deviations. 

N(t)  is  unity  variance  white  noise. 

3.2.3  Gimbal  Angular  Velocity  Isolation 

The  foregoing  equations  for  the  I MU  error  sources,  IMU  case-to- 
navigation  frame  misalignment  angles  and  gyro  frame  commanded  angular 
velocity  permit  calculation  of  the  angular  velocity  of  the  platform  frame 
relative  to  an  inertial  frame.  The  IMU  gimbals  and  closed-loop  servo- 
control  mechanism  function  so  that  this  angular  velocity  is  isolated 
from  that  of  the  IMU  case.  Assuming  that  the  IMU  is  rigidly  mounted  to 
a static  test  fixture,  IMU  case  angular  velocity  is  simply  that  of  the 
rotating  earth.  To  achieve  IMU  pi atform-to-case  angular  velocity 
isolation,  the  gimbals  must  rotate  relative  to  one  another  in  a manner 
dictated  by  the  following  equation. 


_p  _p  p r^PG 

WI,P  5 WPG,P  + CPg[^RG,PG 


sK"  . J*G  fCASE 
fC  ASF.,  RG  cCASE  N 


\*)] 


where 

- RG 

'vCASr  RG  1:5  an9u^ar  velocity  vector  of  the  roll  gimbal  frame 
relative  to  the  IMU  case  coordi natized  in  the  roll 
gimbal  frame. 
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RG.PG 


is  the  angular  velocity  vector  of  the  pitch  girnbal 
relative  to  the  roll  gimbal  frame  coordi nati zed 
pitch  gimbal  frame. 


frame 
in  the 


PG 

WRG,PG 


is  the  angular  velocity  vector  of  the  platform  frame  relative 
to  the  pitch  gimbal  frame  coordinated  in  the  platform  frame. 


and  o are  the  time  rates  of  change  of  the  roll,  pitch, 
and  azimuth  gimbal  angles. 


RG  PG  P 

CCASE’  irKl  cpg  are  ttle  three-roll,  pitch,  and  azimuth  gimbal 

frame  coordinate  transformation  matrices  as  defined 
previously  except  that  the  trigonometric  function  arguments 
are  now  the  G , and  o^z  gimbal  angles  rather  than 

the  nominal  gimbal  synchro  angle  readings. 


The  roll,  pitch,  and  azimuth  synchros  provide  a measurement  of  the 
correspond < ng  gimbal  angles  relative  to  each  other  or  the  1MU  case.  The 
absolute  accuracy  of  the  KT-73  IMU  synchro  devices  is  characterised  by  a 
6-arc-minute  standard  deviation.  The  synchros  are,  however,  capable  of 
5-arc-second  precision  for  sufficiently  small  deviations  about  any 
nominal  gimbal  angle.  To  fully  utilize  the  synchro  angle  measurement 
information,  each  gimbal  angle  is  subdivided  into  three  components. 

0g  (t)  = 8 + &6  + d 0{t) 
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where 

6g(t)  is  the  total  time-varying  gimbal  angle. 

8 represents  the  nominal  synchro  gimbal  angle  used  in  orientating 
the  platform  relative  to  the  I MU  case.  0 is  considered 
fixed  after  initial  platform  orientation. 

A8  represents  a six-arc-minute  bias  angular  error  inherent 
between  0 and  the  actual  gimbal  orientation,  due  to 
synchro  measurement  accuracy  limitations. 

cl 0 ( t)  represents  the  differential  time-varying  gimbal  angle 

component  measurable  to  5-arc-second  accuracy  by  using  the 
inherent  precision  of  the  synchros. 

3.2.4  Linearized  Error  Equations 

The  foregoing  system  of  equations  given  in  Sections  3.2.1  through 
3.2.3  can  be  manipulated  to  yield  a set  of  three  simultaneous  nonlinear 
stochastic  differential  equations  which  relate  I MU  error  sources  and 
platform-to-navigation  frame  misalignment  angles  to  system-level  gimbal 
synchro  angle  measurable  quantities.  The  nonlinear  nature  of  the 
equations  preclude  a rigorous- solution  and  for  this  reason  a first-order 
linearization  of  the  equations  is  assumed  to  be  an  accurate  approximation. 
This  is,  of  course,  a major  assumption  that  will  be  justified  in  the 
sequel.  The  manipulation  and  linearization  of  the  system  of  equations 
by  the  use  of  partial  derivatives  is  a straightforward  process  but  is 
involved  and  tedious,  and  for  these  reasons  is  not  illustrated  here. 

The  set  of  three  simplified  simultaneous  stochastic  differential 
equations  that  result  from  this  process  is  given  below  for  the 
differential  roll,  pitch,  and  azimuth  gimbil  angles. 
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d*(9x  = C{dz)  [ wAx  + wBx+  wTx+wRLx  + wRS  J /C(0y) 
-S^z)[wAy  + WBy  + WTy  + WRLy+WRSy]/C(0y) 
+ C(0?)[-gxMSjt  + gzMjJ/C{6y) 

+ S(02)[gyMSy-gzM,y]  /C(0y> 

^ ^y) QyQzKQx/C lOy) ~S(Qz) gxgz  K Qy / C ( 0y ) 

+ C(0z)PRFxAKFx/C(0y)-S(0z)PRFyAKFy/C(0y) 

+ C(0z)[-WCy£xz+WCz£xy]  /C(6?y) 

+ S(0z)[-wCx  /3yz  + Wqz  /3yxj  / aBy) 
-X2T(0y)S(0x)[£;{L)(qx  + A0x  + d0x)-C{Ui!z] 

+ ^[T(0y)C(0x)C(L)  + S(L)]gy 
+ {r(0y)[c(x)WCx-S(0z)WCy]  /C(0y) 

+ &C(0x)S{L)/CZ(0y)}  [A0y  + d0y] 

“ [s(0z)  WCx  + C(0Z)  WCy] [ A0z  + d0z]  /C{Gy) 

d8y  =S(0z)[wAx  + wBx+wTx+wRU+wRSx] 

c<0z)[wAy+wBy  + wTx+wRLx+wRSx] 

+ S(0z)[-gxMSx  + gzM,x] 

•+C((9z)[-gyMSy+gzM|y] 
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+ S(0z)gygzKOx  + C(Sz)gxgzKOy 
+ S(0z)PRFxAKFx  + C(#z)PRFyAKFy 
+ S(0z)[-WCy/3xz+WCz^xy] 

+ C <02>[  wCx/3xz-wCz^yx] 

-X2  C(0X)  [s(L){77x  + A 8%  + d Qx ) - C ( L)  77  2 J 
-Cl  S(#k)C(L)77x 

+[c  ( ez)  wCx  - S(0Z!  WCy  ] [ A0Z  + d£z] 

d0z  = -T(^y)C(0z)[wAx  + WBx+WTx  + WRLx+WRSx] 

+ T(dy)S(0z)[wAy  + WBy  + WTy  + WRLy+WRSy] 

+ [wAz  + WBz+WTz+WRLz  FWRSzj  fgzMSz-gyM|z 
+ T(0y)C(0z)^gxMgx  -gz  M|  x j 
+ T(0y!S(02)[-gyMSy+gzM|y] 

- T(0y)C(0z)gygzKOx  + T(0y)S(0z)gxgz  K0y 
+ 3x9y»<Or-T^y)C(0i)  PRFxAKPx 
+ T(0y)S(i9z)PRFy  AKFy  >PRFzAKFz 
+ T(0y)C(5I)[wCy,8xz-WCz/3xy] 
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+ T(£y)S(2z)[wCx/3y,-wCz/3yx] 

" ^Cx/^zy  + ^Cy 

+ Q.S{9x)[s{L)(t)x+£9x  + d ^x) -- C ( L) 77 z ] /C(£y) 

-flC(0x)C(L)77y/C(0y) 

-{[c,(02)WCx-S(0z)WCy]  /C2{0y) 

-&T(0y)C(0x>S(L)/C(0y)}[A0y  + d£y] 

+ T(0yi[s(0z)WCx  + C(02)WCy][A02  + d02] 

where  S(  ),C(  ) ond  T(  ) represent  the  trigonometric 
sine, cosine , and  tangent  functions. 

gx  = g[c{0„)  S(0y)C(0z)-S(0x)S(0z)] 

gy  :-~g[c<0x)S(0y)S(0z)  + S(0x)C(0z)] 

gz  = -gC(0x)C(0y) 

wCx  = &{c(L)C(0z,C(0y) 

-S(L)[c(  0Z)  S(0y'  C(0X)~  S(0Z)  S (0X)  j 

WCy  = >Q,{-C(US(6,I)C(0y) 

+ S(  L) [s( 0Z)  S(0y)C(0x)  + C (02)  S(  0X)  j } 

WCz  = &jc(L)S(0y)  + S(L)C(0x)C(0y)] 
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Sign  conventions  for  the  error  equation  development  assume  that  the 
navigation  frame  x,  y,  and  z axes  are  originated  north,  west,  and  up, 
respectively.  The  I MU  ro1!,  pitch,  and  azimuth  gimbal  neutral  positions 
are  coincident  with  the  q y,  and  z navigation  frame  axes,  respectively. 

All  gimbal  rotations  have  positive  directions  as  specified  by  the  right- 
hand  rule. 

3.3  VELOCITY  MEASUREMENT  EQUATIONS 

3.3.1  The  Specific  Force  Equation 

Derivation  of  the  equations  relating  system- level  IMU  velocity 
measurements  to  IMU  error  sources  and  piatforn-to-navigation  frame 
misalignment  angles  follows  directly  from  the  specific  force  equation 
coordinatized  to  the  accelerometer  frame.  Since  the  T M1 " s case  will  be 
rigidly  fixed  to  the  earth  during  calibration  testing.,  the  navigation 
frame  is  considered  fixed  to  the  earth  and  the  appropriate  specific  force 
equation  is  given  as  follows. 

f*A  = ®A'  (PN  RNA  ) " (2*IN  X Pn  RNa) 

where 

is  the  position  vector  from  the  navigation  frame  origin  to  the 
accelerometer  test  mass. 

P^l  is  the  derivative  with  respect  to  time  tafen  relativ°  to  the 
navigation  frame. 

The  last  two  terms  on  the  right-hand  side  of  this  equation  represent 
tangential,  centripetal,  and  coriolis  acceleration  specific  force 
components  produced  by  a finite  accelerometer  test  mass  displacement 
from  the  origin  of  the  accelerometer  frame  and  both  an  angular  acceleration 
,rd  velocity  of  this  frame  relative  to  the  earth  fixed  navigation  frame. 
These  components  can,  however,  be  shown  to  be  insignificant  relative  to 
the  accelerometer  inertial  instrument  random  disturbances  ror  the  KT-73  *MU. 
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Under  these  conditions,  the  accelerometers  provide  a measure  of  the  earth's 
gravity  vector  coordina ti zed  in  the  accel erometer  frame. 


C 


p rPGrRG  -CASE  -N 
PtTRG^CASE  N 9 


where 

C’p  is  the  platform-to-a'..  ?1  erometer  frame  nori-orthogonal  transformation 
matrix  to  account  for  physical  misalignments  of  the  acceler- 
ometer sensitive  axes  from  the  platform  axes. 


4 


1 

y -r 

KZ  Xy 

-/ 

! y 

'y* 

y# 

r 

-y  f 

zx 

its  of 

this  matrix 

given  in  Figure  3-1  for  the  gyre  to-platform  mi  sal ignmerit  matrix.  The 
first  subscript  refers  to  the  misaligned  axis  ana  the  second  subscript 
refers  to  the  axis  about  which  a small  angular  rotation  has  occurred. 


3.3.2  Accelerometer  error  Model 

In  addition  to  ideally  providing  a measu"e  of  the  specific  force 
vector,  the  acce (erometer  outputs  are  corrupted  by  various  hardware 
induced  errors.  The  composite  accelerometer  output  signals  are  net, 
however,  measurable  at  systetc- level  for  the  KT-73  I MU  In  the  hardware 

mechanization  of  this  IMU,  the  accelerometer  output  signals  are 
integrated  in  capacitive  reset  integrators  which  ir.de;.  velocity  registers 
w'th  delta  velocity  pulses.  The  accumulated  velocity  pulses  are 
subsequently  multiplied  by  airborne  computer  loaded  accelerometer  scale 
factors  to  obtain  total  velocity  in  the  desire.)  units. 


The  digitization  of  the  accelerometer  output  signal  into  a train  of 
pulses  whose  rate  represents  the  magnitude  of  the  accelei afion  int  duces 
a velocity  quantization  level  which  represents  the  change  in  velocity  per 
pulse.  Since  analog  velocity  data  is  not  continuously  available  at  any 
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instant,  of  time,  the  quantization  process  introduces  an  error.  If 
the  quantization  effects  are  treated  as  errors  in  the  system- level  I MU 
velocity  measurement: , then  the  accelerometer  output  signals  can  be 
useo  to  define  an  artificial  time  derivative  velocity  vector  modeled 
as  fo flows. 


V 


+ Aa  + ar  + 


A * 
TJ 


[K^]  is  a 3x3  diagonal  matrix  of  the  airborne  computer  loaded 

acce lerometer  scale  factors.  The  true  accelerometer  scale 
factors  [Kpj]  differ  from  the  loaded  ones  by  the  sum  of 


[AKa]  is  considered  as  the 


two  errors,  [AK^’l  anc;  [AK^a ] . 
nominal  accelerometer  scale  factor  error  as  specified  in 
Reference  1.  [AK^]  is  an  additional  error  source  showing 
a platform  attitude  dependency  which  has  been  added  to  the 
Reference  1 error  model  to  account  for  potential  discrepancies 
between  the  assumed  and  actual  accelerometer  error  models. 
Since  [AK.^]  can  acquire  new  valuta  each  time  the  platform 
changes  attitude  with  respect  to  the  I MU  case,  uthcr 
assumptions  regarding  the  origin  of  this  error  model 
residual  would  be  equally  acceptable. 

is  the  acceler ometer  bias  error  vector. 

is  the  exponential ly  correlated  accelerometer  random  noise 
vector. 


- A (t) 

adO)  = ^ 


N(») 


is  the  random  noise  correlation  time, 
is  the  random  noise  standard  deviation. 


N(t)  is  unity  variance  white  noise. 
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A 


T 


i 

o 

t 

tAT 


is  the  accelerometer  thermal  warm-up  bias  vector. 

(» • To) 


is  the  initial  accel erometer  warm-up  bias  vector. 

is  the  time  from  power  turn-on  to  the  start  of  fine  gyro- 
compass alignment. 

is  the  time  from  power  turn-on. 

is  the  gyro  warm-up  drift  time  constant. 


A special  hardware  feature  of  tne  KT-73  I MU  is  the  availability  of 
toth  high  and  low  gain  accel erometer  scale  factors.  The  nigh  gain 
accelerometer  scale  factor  of  100,000  pulses-per-second  per  "g"  (pps/g) 
is  used  for  fine  platform  alignment  tc  redu  e velocity  quantization  errors 
and  thereby  shorten  the  necessary  alignment  imes.  The  low  gain  scale 
factor  of  1,000  pps/g  is  used  for  navigation  since  high  accelerations 
would  saturate  tha  high  gain  capacitive  reset  integrator's  maximum 
OL'cput  of  10,000  pfs.  High  gain  scale  factor  modes  are  mechanized  only 
as  required  on  the  two  level-axis  accelerometers  and  separate  accelerometer 
bias  errors  are  associated  with  the  nigh  and  low  gain  scale  factor  modes. 
This  TMU  mechanization  results  in  four,  two  bias  and  two  scale  factor, 
additional  accelerometer  errors  tc  be  calibrated. 


3.3.3  Linearized  t>ror  Equations 


In  contrast  to  the  gimbal  angle  time  rates  of  change,  the 
accelerometer  outputs  dc  not  hove  an  error-free  valuQ  cf  zero.  A zero 
mean  accelerometer  error  vector  equation  can,  however,  be  defined  as 
follows. 


1“  - P PG  RG  _N 

Ay  - V - C C C Q 

V PG  PG  CASE  * 


All  terms  arc  as  previously  defined  with  the  gimbul 
transformations  being  computed  as  functions  of  the  initia 


synchro  angle  settings 


frame  coordinate 
1 6y,  0 , and 
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The  system  of  equations  given  in  3.3.1  through  3.3.3  can  be  solved 
to  yield  a set  of  three  simultaneous  nonlinear  stochastic  differential 
equations  which  relate  the  I MU  error  sources  and  platform- co-navi gation 
frame  misalignment,  angles  to  system-level  I MU  velocity  measurements. 
These  equations  are  manipulated  and  linearized,  in  a manner  similar  tc 
that  for  the  gimbal  synchro  measurement  equations,  such  that,  a solution 
can  be  obtained.  The  resultant  error  equations  for  the  x,  y and  i 
platform  axes  velocity  measurements  are  given  below. 

Avx  = - g[cv£z)S(6y)S(0x)  + C{£x)l  77r 

+ g[c(02)c(0y)j 

- g[s{62)C(0x)  + C(02)S(9/’S{&X)][ A0x  + d£x] 

+ g [c<  <9Z  J c(0y ) C(  + d£y] 

- g[s(02)s(£y)c(£x)  +c {Bz)  s(£x)][a0z  + n 9Z\ 

+ g[c(6':)S{£y5c(£x)  - 3(£z)s»£x)][akAx  + £kaa<] 

- g[s«9z)s(£y/C(£x)  t c(0z)s(0x)]  y%t 
~ g j^ct£y)c(£x) j 

+ ABs  + ARx  + ATk 
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AVy  = g[3(^)S(0y)S(ex)  - C(02)C{Ox)]t;x 

- g[s{^)C(a/)]  ^ 

+ g[si02)s{0y)s(0x)-c(02)G(6x)][A0x  + d#,] 

- g[s{02)C{0y)c(0x)][A0x  + d£x] 

+ e[s(02)S(0x)-C(02)S(0y)C(0x)][A02  + d02] 

- g[c(02)S(0J  + S(a2)S(ay)C{0x)][AKAy  + AkaaJ 
+ g [cva2)  ${8,)ci6n)  - s($2)  s(ax)]  yn 
-g[c(ay)c(ax)]yyx 

+ ABy  + ARy  + AIy 

Av2  = g [c(0y)  S(0X)]  77x 
+ g S(ay)^y 

+ g[c(0y)s(0x)][A0x  + d0x] 

+ g[s(ay)c(ax)][Aay  + d0y] 

-g[c<ay)c(ax)][AKAl  + AkAAz] 

+ g[c(02)S{0y)C(0x)  -S(02)  S(0„)]y 

- g [s(02)  S[9y)  C(0X)  + C(0t)S(ax)]yM 
+ ABi  + ARz  + ATz 


A. 
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The  functional  dependencies  for  the  level-axes  high  and  low  gain 
accelerometer  scale  factors  and  correspondi ng  biases  are,  of  course,  the 
same.  In  the  above  equations,  only  one  error  source  is  explicitly  shown 
but  is  considered  to  represent  either  the  high  or  low  gain  error, 
depending  upon  which  mode  is  being  used  at  a particular  instant  during 
the  test. 

Note  that  while  the  velocity  error  equations  appear  somewhat  simpler 
than  those  for  the  corresponding  gimbal  angles  in  the  above  formulation, 
the  velocity  error  equations  are  functions  of  the  differential  gimbal 
angles  and,  therefore,  functions  of  all  the  I MU  and  platform  alignment 
error  sources. 
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SECTION  IV 

SEPARATION  OF  IMU  ERRORS 

4.1  CONCEPTUAL  APPROACH  TO  ERROR  SEPARATION 

The  preceding  derivations  have  resulted  in  a set  of  six  general 
simultaneous  equations  relating  the  six-system-level  IMU  measurable 
quantities  to  roughly  sixty  error  sources.  Although  a maximum  of  six 
measurements  and,  hence,  equations  are  available,  IMU  calibration  tests 
could  be  conducted  where  only  some  of  the  measurements  were  utilized. 

If,  for  example,  the  gimbal  synchro  angle  measurements  were  not  made,  then 
only  the  three- vel oci ty-measurement  equations  could  be  employed  to  solve 
for  the  error  sources.  Since  a set  of  six  or  fewer  specific  simultaneous 
stochastic  differential  equations  cannot  be  solved  for  sixty  unknowns, 
a means  functionally  separating  the  error  quantities  had  to  be  formulated. 
In  view  of  the  stochastic  nature  of  the  equations  and  IMU  measurement 
noia,  stochastic  estimation  of  the  error  quantities  is  the  best  solution 
possible. 

The  conceptual  approach  to  error  source  separation  is  to  generate 
multiple  sets  of  specific  simultaneous  differential  equations  from  the 
basic  set  of  six  general  ones.  If  not  all  the  available  IMU  measurements 
are  to  be  used,  then  the  basic  set  of  general  equations  would  be 
correspondingly  reduced.  By  generating  enough  multiple  sets  of  specific 
equations  such  that  the  total  number  of  equations  is  equal  to  the 
number  of  unknowns,  accurate  stochastic  estimation  of  IMU  error  sources 
should  be  possible. 

The  generation  of  the  multiple  sets  of  specific  equations  is 
accomplished  by  angularly  orientating  the  platform  at  various  attitudes 
relative  to  the  earth  rate  and  gravity  vectors.  Note  that  a given 
platform  attitude  relative  to  these  vectors  is  obtained  by  a specific 
set  of  gimbal  angle  rotations.  The  specific  gimbal  angles  used,  when 
inserted  into  the  set  of  general  equations,  results  in  a set  of  specific 
equations.  This  approach  follows  directly  from  the  observation  that  the 
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measurable  I MU  parameters  are  functions  of  both  the  I MU  error  sources 
and  the  coord i ni zation  of  the  above  vectors  in  the  platform  and  inertial 
instrument  reference  frame. 

In  developing  this  approach,  it  is  convenient  to  group  the  error 
quantities  into  three  groups: 

(1)  The  first  group  consists  of  the  I MU  error  sources  which  are 
modeled  as  random  biases  and  instrument  thermal  transients  which  are 
commonly  denoted  as  calibration  parameters. 

(2)  The  second  group  consists  of  those  angular  misalignment  errors 
associated  with  platform-to-navigation  frame  alignment. 

(3)  The  third  group  consists  of  the  stochastic  processes  modeling 
the  inertial  instrument  random  disturbances. 

4.2  I MU  CALIBRATION  PARAMETERS 


A necessary,  but  not  necessarily  sufficient,  condition  for  separation 
of  the  IMU  error  sources  is  that  the  multiple  sets  of  system- level  IMU 
measurement  equations  must  be  solvable  for  the  IMU  calibration  parameters 
after  the  equations  are  simplified  to  contain  only  the  calibration  error 
sources.  This  simplification  is  accomplished  by  omitting  the  error 
sources  associated  with  the  above  second  and  third  groups  from  the 
general  error  equations.  Performing  this  simplification  reduces  the 
set  of  six  simultaneous  stochastic  differential  equations  to  a system 
of  six  simultaneous  ordinary  differential  equations. 

This  system  of  ordinary  differential  equations  can  be  indefinitely 
expanded  by  positioning  the  platform  at  new  attitudes  as  pointed  out 
Section  4.1.  A deteririi 'istic  solution  to  the  multiple  sets  of  the 
ordinary  differential  equations  can  be  attempted  by  application  of 
Laplace  transforms  and  using  the  initial  and  final  IMU  measurements  for 
each  attitude. 
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To  confirm  the  conceptual  approach  to  error  source  separation  and 
to  determine  the  most,  relevant  IMU  platform  test  attitudes  for  conducting 
the  calibration  tests,  a solution  to  the  simplified  system  of  equations 
was  attempted  as  described  above.  To  facilitate  the  analysis,  several 
additional  error  sources  ■■""e  omitted.  The  gyro  and  accelerometer  thermal 
transient  warm-up  errors  are  deleted  since  they  decay  rapidly  with  time. 
Only  low  gain  accel erometer  scale  factor  and  bias  errors  were  considered 
since  the  errors  associated  with  the  high  gain  mode  could  be  separated 
in  a similar  model.  By  necessity,  the  gyro  and  accelerometer  error 
model  residuals,  which  had  oeen  added  to  allow  for  the  possibility  that 
the  actual  IMU  error  model  is  significantly  different  from  that  assumed, 
were  deleted  since  they  assume  a new  value  for  each  new  platform  attitude 
used  and  would  thus  prevent  solving  the  system  of  equations. 

Table  4-1  illustrates  how  the  IMU  calibration  error  sources  become 
separable  as  a function  of  the  cumulative  number  of  platform  test 
attitudes  employed.  Both  sets  of  velocity  and  aimbal  synchro  angle 
measurements  were  assumed  to  be  available.  The  x's  indicate  the  earliest 
point  at  which  the  generated  multiple  sets  of  equations  could  be  deter- 
ministically solved  for  a particular  error  source.  All  of  the  33  IMU 
calibration  errors  could  be  determined  after  7 platform  attitudes  times 
6 specific  equations  per  attitude  or  42  equations  were  generated.  Since 
only  33  equations  should  ideally  be  required,  not  all  new  equations 
generated  are  independent  of  previous  ones.  The  specific  attitudes 
shewn  were  determined  by  a combination  judgment  and  trial -and-error 
process. 

4.2.1  Girnbal  Angle  Measurement  Requirements 

The  KT-73  IMU  is  mechanized  with  synchro  control  transmitters  and 
transformers  having  an  accuracy  of  6 arc  minutes  standard  deviation. 
Contemporary  synchro-to-digital  converters  have  a 12  to  13  bit  conversion 
capability  which  results  in  a precision  on  the  order  of  1 to  3 arc 
minutes.  Since  the  precision  of  the  synchro-to-digital  converter  is 
greater  than  the  accuracy  of  the  synchro,  the  equipment  combination  is 
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generally  suitable  for  most  inertial  navigation  applications.  To 
accurately  measure  platform  drift  rates,  however,  a synchro  angle 
measurement  precision  on  the  order  of  the  precision  of  the  synchro  device 
itself  of  5 arc  seconds  is  desirable.  Since  a special  measurement 
technique  would  nave  to  be  developed  to  achieve  this  precision,  it  would 
be  desirable  to  calibrate  the  I MU  by  use  of  only  the  velocity 
measurements . 

As  discussed  in  Section  3.3.3,  the  velocity  measurement  equations 
contain  the  differential  gimbal  angle  errors  and  thus  all  the  IMU 
calibration  error  parameters.  To  address  the  question  of  the  feasibility 
of  calibrating  the  IM'J  by  use  of  only  the  velocity  measurements , the 
Laplace  transform  analysis  was  repeated  using  only  these  measurements. 
Results  of  this  analysis  showed  that  all  calibration  errors  could  be 
determined  except  the  gyro  spin  axes  mass  unbalance  errors. 

The  Laplace  analysis  results,  however,  raise  a question  as  to  the 
reason  for  the  inseparability  of  the  spin  axes  mass  unbalance  errors.  Is 
the  inseparability  due  to  failure  to  select  the  required  combination  of 
platform  test  attitudes  or  is  there  a fundamental  limitation  inherent  in 
the  conceptual  approach?  To  answer  this  question,  the  basic  form  of  the 
combined  velocity  and  gimbal  synchro  angle  measurement  equations  was 
examined  for  linear  independence.  To  facilitate  the  analyses,  all  error 
sources  but  the  spin  axes  mass  unbalance  terms  were  omitted.  One 
representation  of  the  resultant  set  of  three  simultaneous  equations  is 
as  follows: 


AV«  =’2AC(8r)c(«,)K,  + »„] 

Avy  = g2ec(8,)c(9>)[„s<+  „SJ 

Ms,] 
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where 

A =[s(®i)s(Sy)c(fi«)  +c(e,)s(e*)] 

B = [c(fi,)s(9,)c(9,)-s(9J)s(8,)] 

These  equations  relate  the  rate  of  change  of  the  specific  force 
vector,  as  coordinati zed  in  the  accelerometer  frame,  to  the  spin  axes 
mass  unbalance  error  coefficients.  The  equations  a>-e  not  independent, 
as  can  be  shown  by  combining  the  first  two  to  yield  an  equation  directly 
proportional  to  the  third. 

bAVx  - A A Vy  = 2g2  AB  [ms„  - Mgx  J 

Several  conclusions  can  be  reached  as  a result  of  this  analysis. 
Without  a gimbal  synchro  angle  measurement  or  its  equivalent,  separation 
of  the  gyro  spin  axes  mass  unbalance  errors  is  not  possible.  Only 
combinations  of  these  errors 

(MS,  + M5i)  • (MS,  + MSl)  °"d  (MS,  - MS«) 

can  be  determined.  A single  piece  cf  relevant  information  is  required 
in  addition  to  the  velocity  measurements  to  separate  the  spin  axes  mass 
unbalance  errors . This  information  can,  however,  consist  of  one  gimbal 
synchro  angle  measurement  for  a particular  test  attitude  rather  than 
require  all  three  synchro  gimbal  angle  measurements  for  all  test  attitudes. 

4.3  PLATFORM  ALIGNMENT  ERRORS 

The  analyses  of  Section  4 . 2 showed  that  all  IMU  calibration  parameters 
could  be  determined  if  the  second  and  third  groups  of  errors  defined  in 
Section  4.1  were  omitted  from  the  measurement  equations.  The  first  order 
effect  of  introducing  the  platform-to-navigatiori  frame  alignment  errors 
of  the  second  group  into  the  solution  for  the  separation  of  the  IMU 
calibation  oarameters  of  the  first  group  can  be  understood  with  reference 
tc  figure  4-1.  For  a small  platform- to-navigation  frame  misalignment,  the 
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effect  is  precisely  the  classical  problem  of  separating  west  gyro  drift 
from  the  west  platform  axis  commanded  angular  velocity  error  due  to  the 
azimuth  misalignment  angle  i p. 

To  null  the  angular  velocity  of  an  ideal  I MU  platform  with  resoect 
to  the  navigation  frame,  the  gyros  must  be  commanded  to  process  at 
exactly  the  magnitude  of  the  earth  rate  vector  coordinatized  in  the  gyro 
frame.  Since  ip  is  a zero  mean  random  variable  whose  magnitude  is  net 
known  a priori,  the  west  platform  axis  would  not  be  commanded  precess. 

For  this  condition,  the  angular  velocity  of  die  west  platform  axis  with 
respect  to  the  navigation  axis  is  approximately  given  by: 

"n.p,  a *(•-)  + 

Given  the  standard  deviation  of  the  azimuth  misalignment  angle  cm  and  a 
perfect  measurement  of  the  standard  deviation  ow  of  the  W.  ^ 

estimation  error  can  be  determined  a„: 

® c ( L ) 

This  inability  to  separate  the  west  platform  commanded  angular 
velocity  error  from  the  composite  west  gyro  drift  does  not  permit 
separation  of  all  the  IMU  calibration  errors  as  predicted  in  Table  4-1. 
The  error  separation  problem  can,  however,  t?  alleviated  by  realizing 
that  the  north  and  up  platform  axes  do  not  encounter  this  first  order 
limitation  due  to  a different  commanded  angular  ve'ocity  attitude 
dependency.  The  solution  ;s  tc  add  additional  platform  test  attitudes 
until  all  IMU  calibration  parameters  become  separable.  It  was  determined 
by  the  Laplace  analysis  Lnat  by  adding  three  now  test  attitudes,  similar 
to  the  first  three  in  Table  4-1  except  that  the  x and  y platform  axes 
were  interchanged,  the  calibration  parameters  became  separable. 

4.4  INERTIAL  INSTRUMENT  RANDOM  DISTURBANCES 

incorporating  the  inertial  instrument  random  disturbances  of  +he 
third  group  of  IMU  error  sources  into  the  system  of  measurement  equations 
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i har.qas  the  problem  description  from  one  of  ordinary  differential  equations 
to  stochastic  differential  equations.  The  nature  of  the  stochastic 
equations  and  the  presence  of  I MU  measurement  noise  precludes  a deterministic 
solution  to  the  separation  of  I MU  errors  and  dictates  the  use  of  some 
form  of  stochastic  estimation  if  optimal  results  are  to  be  obtained 

Separation  of  some  of  the  IMU  calibration  parameters  from  the  inertia! 
instrument  random  disturbances  will  have  to  be  made  on  the  basis  of  time 
rather  than  platform  attitude  since  it  is  the  only  difference  in 
functional  dependency  exhibited. 
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SECTION  V 


KALMAN  FILTER 


5.1  STOCHASTIC  ESTIMATION 

Selection  of  the  technique  for  estimating  the  I MU  errors  wa r,  Made  in 
consideration  of  several  major  criteria. 

(1)  The  technique  should  be  optimum  in  the  mini mum  variance  sense. 

(2)  The  technique  should  provide  a statistical  Measure  of  the 
accuracy  to  which  the  I MU  hardware  errors  can  be  call  ora  ted. 

(3)  The  technique  should  provide  a convenient  means  cf  validating 
or  justifying  the  first  order  measurement  equation  1 inearizatie ns  and 
other  simplifying  a.-. sumptions. 

(4)  The-  technique  should  provide  a comprehensive  solution  by 
evaluating  all  error  sources  in  u parallel  or  simultaneous  mariner  to 
enhance  the  possibility  of  defecting  unmodeled  I MU  hardware  errors . 

(5)  The  technique  should  be  capable  of  being  readily  modified  to 
include  additional  error  sources  or  facilitate  conoucting  correlation 
studies. 

Although  many  optimal  stochastic  estimation  techniques  are  available, 
the  Kalman  Filter  algorithm  with  its  associated  covariance  matrix  was 
considered  to  be  the  most  applicable  approach.  Conversely,  the  formulated 
problem  sat.sfier.  the  linearity,  a priori  system  dynamics,  and 
Gauss-Markov  stochastic  process  modeling  restrictions  under  which  the 
Kalman  Filter  algorithm  was  derived. 

5.2  KALMAN  FILTER  IMPLEMENTATION 

5.2.1  Filter  Equations 

The  error  equation  development  results  of  Section  3 are  ce  <s  is  tent 
with  state  space  analysis  notation  which  is  a convenient  means  ot 
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mathema c i <~a  1 ly  representing  a physical  system  in  the  time  domain.  The 
description  of  a continuous  linear  system  by  state  variables  is  done  by 
means  of  n first-order  simultaneous  differential  equations.  By  using 
vectors  and  matrices,  the  I MU  calibration  problem  can  be  defined  oy  the 
given  and  derived  error  equations  formulated  in  a vector  differential 
creation. 

i {♦)  = K(t)  x (t>  + G { t ) w {)} 

where 

x(t)  is  the  system  n dimension  state  vector 
F(t)  is  the  (nxn)  dirner.sion  system  dynamics  matrix 
G^t)  is  the  (nxp)  noise  distribution  matrix 
w(t)  is  the  p dimension  disturbance  vector 

During  the  IMIJ  calibration  tests,  a set  of  discrete  measurements 
will  be  made  at  time  t . The  corresponding  measurement  process  can  be 
described  as  a vector  algebraic  equation;  thus, 

S(M  s H(,n)  * (tn)  + v(»n) 

where 

z(t  ) is  the  m dimension  measurement  vector 

H(t  ' is  the  (mxn)  measurement  distribution  matrix 

v(tn)  is  the  m dimension  measurement  noise  vector 

The  mean  and  covariance  of  the  assumed  white  noise  processes  w(t) 

and  v(t  ) are  given  as: 
ri  3 

E [*(')]  =° 

E [«(t)  wT(t  + r)]  • 0{l)a{r) 
e [*(»«)]  *° 
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where 

3(t)  and  a(t  ~t  ) are  dirac  delta  functions 
Q(t)  is  che  (pxp)  system  noise  matrix 
R(t  ) is  the  (mxm)  measurement  noise  matrix 

The  Kalman  Filter,  References  2 and  3,  provides  an  optimal  (minimum 
variance)  estimate  of  x(t),  which  is  denoted  as  x(t).  The  optimal 
estimate  of  the  state  vector  x(c)  and  its  (nxn)  covariance  matrix  P(t) 
are  propagated  netween  measurements  by  the  equations. 

$(•>  - F(t)S(t) 

p(0  * f'(t)  p(i)  +■  + «<»)  o(i)gt(i) 

With  the  initial  conditions  at  time  zero  t.0, 

* M = E[*(,o)] 

P(t0)  - e{[*(t0)-*(t0)][*(t0)  - x(t0)]  } 

At.  measurement,  times  t , the  Kalman  Gain  K(t  ) is  determined  and 

n n 

the  state  vector  optimal  estimate  and  its  covariance  matrix  updated 
according  to  the  expressions, 

*(»n)  - p/V)HT(»n)[H{»n)'r  (tn')HT(^n)  + R(»n)] 

*(C)  = *(»n"}+»<(tn)[t(fn)-H(tn)^(tn")] 

P(tn+)=[t-K(tn)H(tn)]p{tn“)[t  -K(tn)H(t„)]  +K(t,,)R(tr,)KT(tn) 

5.2.2  Filter  Vectors  and  Matrices 

Formulation  of  the  error  equations,  given  or  derived  in  Section  3, 
in  state  variable  vector  and  matrix  notation  is  a straightforward 
procedure.  Initially,  both  the  gyro  and  accelerometer  short  correlation 
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time  random  disturbances  were  included  in  the  Kalman  Filter  state  vector. 
Subsequent  computer  simulation  studies  showed,  however,  that  the  filter 
could  neither  reduce  their  uncertainty  nor  estimate  their  magnitudes 
acxurately  due  to  their  relatively  short  correlation  times  as  compared 
to  the  measurement  time  intervals.  These  errors  were,  therefore,  deleted 
from  the  state  vector  and  incorporated  as  disturbances  on  the  differential 
gimbal  angle  and  velocity  states.  The  mechanized  Kalman  Filter  68  state 
vector  is  shown  in  Table  5-1.  The  location  of  individual  states  within 
the  state  vector  has  no  particular  significance  and  is  partly  a result  of 
the  iterative  development  process.  The  error  sources  associated  with 
filter  states  4 through  6 and  68  were  not  included  in  Section  3 because 
they  were  discovered  during  actual  laboratory  testing  of  t(ie  KT-73  IMU. 
These  errors  will  be  discussed  further  in  Section  8. 

The  great  majority  of  the  nonzero  68  x 68  system  matrix  F elements 
are  located  in  the  F matrix  rows  correspond!' ng  to  the  system-level 
differential  gimbal  angle  do  and  pseudo  velocity  AV  measurement  states. 

Due  to  the  rather  large  dimensionality,  this  matrix  is  not  readily 
pictorially  illustrated.  Specific  F matrix  elements  are,  however, 
readily  identified  from  the  linearized  differential  gimbal  angle  and 
velocity  error  equations  tabulated  in  Sections  3.2.4  and  3.3.3. 

With  the  exception  of  the  three  states'  modeling  IMU  case-to- 
navigation  frame  misalignment  angles,  all  state  vector  elements  modeled 
as  random  constants  are  considered  as  IMU  calibration  parameters,  and 
do  not  contribute  any  nonzero  F matrix  elements.  The  F matrix  contains 
the  negative  reciprocals  of  the  correlation  times  and  time  constants 
as  diagonal  elements  in  the  rows  corresponding  to  the  states,  modeled 
as  exponentially  correlated  random  noise  and  random  constants  with  an 
exponential  decay.  The  F maxtrix  elements  are  independent  of  time  but 
some  do  change  in  magnitude  each  time  the  IMU  platform  is  re-oriented  due 
to  the  nominal  synchro  angle  dependency. 
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TABLE  5-1 

IMU  CALIBRATION  KALMAN  FILTER  STATE  VARIABLES 


MODEL/ 

STATE  INITIAL  la  VALUE 


G-INSENSITIVE 

GYRO  DRIFT 

RANDOM 

CONSTANTS 

1 . 

WBx 

X 

Gyro  Drift  Rate 

0.1o/hr 

2. 

Wn 

By 

Y 

Gyro  Drift  Rate 

O.r/hr 

3. 

Lj 

dz 

Z 

Gyro  Drift  Rate 

0. 5°/hr 

G-SENS ITI VE  GYRO  DRIFT  COEFFICIENTS 

RANDOM 

CONSTANTS 

4. 

Wn 

Qx 

X 

Gyro  Output  Axis  Quadrature  G-Sensi ti vi ty 

0.4°/hr/g 

5. 

wQy 

Y 

Gyro  Output  Axis  Quadrature  G-Sensi ti vi ty 

0.4°/hr/g 

6. 

V 

Z 

Gyro  Output  Axis  Quadrature  G-Sensi ti vi ty 

0.05°/hr/g 

GYRO 

LONG 

CORRELATION  TIME  RANDOM  DRIFT 

EXPONENTIALLY 

CORRELATED 

7. 

WRLx 

X 

Gyro  Random  Drift  Rate 

0.005°/hr 

8. 

WRLy 

Y 

Gyro  Random  Drift  Rate 

0.005°/hr 

9. 

WRLz 

Z 

Gyro  Random  Drift  Rate 

0.005°/hr 

G-SENSITIVE  GYRO  DRIFT  COEFFICIENTS 

RANDOM 

CONSTANTS 

10. 

MSx 

X 

Gyro  Spin  Axis  G-Sensi ti vi ty 

1.0°/hr/g 

11  . 

MIx 

X 

Gyro  Input  Axis  G-Sensi ti vi ty 

0.01°/hr/g 

12. 

MSy 

Y 

Gyro  Spin  Axis  G-Sensi ti vi ty 

1 .0°/hr/g 

13. 

Miy 

Y 

Gyro  Input  Axis  G-Sensi ti vi ty 

0.01°/hr/g 

14. 

MSz 

Z 

Gyro  Spin  Axis  G-Sensi ti vi ty 

0.5°/hr/g 

15. 

Mlz 

Z 

Gyro  Input  Axis  G-Sensitivi ty 

0.01o/hr/g 
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TABLE  5-1  (Contd) 

g2-sensitive 

GYRO  DRIFT  COEFFICIENTS 

RANDOM 

CONSTANTS 

16. 

K0x 

X 

2 

Gyro  Output-Spin  Axes  G Sensitivity 

0.0r/hr/q2 

17. 

K0y 

Y 

2 

Gyro  Output-Spin  Axes  G Sensitivity 

0.01o/hr/g2 

18. 

K0z 

z 

2 

Gyro  Output-Spin  Axes  G Sensitivity 

0.01o/hr/g2 

GYRO 

SCALE 

FACTOR  ERRORS 

RANDOM 

CONSTANTS 

19. 

AKFx 

X 

_3 

Gyro  Scale  Factor  Error  2.8x10 

arc  sec/pulse 

20. 

AKFy 

Y 

_3 

Gyro  Scale  Factor  Error  2.8x10 

arc  sec/pulse 

21 . 

akFz 

Z 

-3 

Gyro  Scale  Factor  Error  2.8x10 

arc  sec/pulse 

GYRO 

INPUT 

AXES  MISALIGNMENTS 

RANDOM 

CONSTANTS 

22. 

Bxz 

X 

Gyro  Input  Axis  Misalignment  About  z Axis 

1 .0  arc  min 

23. 

Bxy 

X 

Gyro  Input  Axis  Misalignment  About  y Axis 

1.0  arc  min 

24. 

Byz 

Y 

Gyro  Input  Axis  Misalignment  About  z Axis 

1.0  arc  min 

25. 

Byz 

Y 

Gyro  Input  Axis  Misalignment  About  x Axis 

I . 0 arc  min 

26. 

Bzy 

Z 

Gyro  Input  Axis  Misalignment  About  y Axis 

1.0  arc  min 

27. 

Bzx 

z 

Gyro  Input  Axis  Misalignment  About  x Axis 

1 .0  arch  min 

IMUJ 

CASE  AXES 

MISALIGNMENTS 

28. 

nx 

Case  Misalignment  About  x Navigation  Axis 

4.0  arc  min 

29. 

ny 

Case  Misalignment  About  y Navigation  Axis 

4.0  arc  min 

30. 

nz 

Case  Misalignment  About  z Navigation  Axis 

4.0  arc  min 
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TABLE  5-1  (Contd) 


DIFFERENTIAL 

GIMBAL  ANGLE 

ERRORS 

DYNAMIC 

31 . 

d0 

X 

Differential 

Roll  Gimbal  Angle  Error 

0.0  arc 

mi  n 

32. 

de 

V 

Di f ferential 

Pitch  Gimbal  Angle  Error 

0.0  arc 

miri 

33. 

de 

z 

Differential 

Azimuth  Gimbal  Angle  Error 

0.0  arc 

mi  n 

ACCELEROMETER  LOW  GAIN  MODE  BIAS 


RANDOM 

CONSTANTS 


34. 

An 

Bx 

X 

Accel erometer 

Low  Gain  Bias 

2,500  ug 

35. 

V 

Y 

Accelerometer 

Low  Gain  Bias 

2,500  pg 

36. 

ABz 

Z 

Accelerometer 

Low  Gain  Bias 

2,500  ug 

ACCELEROMETER 

THERMAL  WARM-UP  TRANSIENTS 

RANDOM  CONSTANT 

F XPONENi  IAL  DECAY 

37. 

At 

Tx 

X 

Accelerometer 

Thermal  Warm-Up  Error 

5.0  pg 

38. 

S 

Y 

Acce1 erometer 

Thermal  Warm-Up  Error 

5.0  pg 

39. 

ATz 

Z 

Accel erometer 

Thermal  Warm-Up  Error 

200  pg 

ACCE! 

.EROMETER 

LOW  GAIN  MODE 

SCALE  FACTOR  ERRORS 

RANDOM 

CONSTANTS 

40. 

AK. 

Ax 

X 

Accelerometer 

Low  Gain  Scale  Factor  Error 

0.03  numeric 

41 . 

aKAy 

Y 

Accelerometer 

low  Gain  Scale  Factor  t\ 

a or 

0.03  numeric 

42. 

AKAz 

Z 

Accel erometer 

Low  Gain  Scale  Factor  Error 

0.03  numeric 

ACCELEROMETER 

INPUT  AXES  MISALIGNMENTS 

RANDOM 

CONSTANTS 

43. 

^xz 

X 

Accelerometer 

Input  .-.x'ls  Misalignment 

About 

z 1 .0  arc  min 

44. 

yxy 

X 

Accel erometer 

Input  Axis  Misalignment  About 

y 1 .0  arc  min 

45. 

yyz 

Y 

Accelerometer 

Input  Axis  Misalignment 

About 

z 1 .0  arc  min 

46. 

Yyx 

Y 

Accelerometer 

Input  Axis  Misalignment  About 

x 1.0  arc  min 

47. 

rzy 

Z 

Accelerometer 

Input  Axis  Misalignment  About 

y 1.0  arc  min 

43, 

yzx 

7 

Accel erometer 

Input  Axis  Misalignment, 

About 

x 1 . 0 a re  mi  n 
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TABLE  5-1  (Conte!) 

PSEUDO  VELOCITY  ERRORS 


49 

• AVW 

X 

Platform  Axis 

Velocity  Error 

50 

. AV 

y 

y 

Platform  Axis 

Velocity  Error 

51 

. AV 

2 

Platform  Axis 

Velocity  Error 

G- 

iNSENST 

TIVE 

GYRO  DRIFT 

52 

■ WAx 

X 

Gyro  Attitude 

Dependent.  Gyro 

Drift 

Rate 

53 

’ WAy 

Y 

Gyro  Attitude 

Dependent  Gyro 

Drift 

Rate 

54 

■ WAz 

Z 

Gyro  Attitude 

Dependent  Gyro 

Drift 

Rate 

GIMBAL  SYNCHRO  ANGLE  MEASUREMENT  ERROR 


55. 

, A0x 

Roll  Synchro  Measurement  Lias 

56 . 

. aey 

Pitch  Synchro  Measurement  Bias 

57. 

. A6Z 

Azimuth  Synchro  Measurement  Bias 

ACCELEROMETER 

HIGH  GAIN  MODE  BIAS 

58. 

ABHx 

X 

Accelerometer  High  Gain  Bias 

59. 

ABHy 

Y 

Accelerometer  High  Gain  Bias 

ACCELEROMETER 

HIGH  GAIN  MODE  SCALE  FACTOR  ERRORS 

60. 

AKAHx 

X 

Accelerometer  High  Gain  Scale  Factor 

Error 

61 . 

AKAHg 

Y 

Accelerometer  High  Gain  Scale  Factor 

Error 

DYNAMIC 

0.0  ft/' sec 

0.0  ft/sec 

0.0  ft/ sec 

RANDOM 
CONSTAN i S 

0.005c/hr 

O.OOSVhr 

O.OOE°/hr 

RANDOM 

C0NSTAM1 

6.0  Ci  c min 

6.0  a*"C  min 

6.0  arc  min 

RANDOM 

CONSTANTS 

2,500  yg 

2,500  yg 

0.03  numeric 
0.03  numeric 
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TABLE  5-1  (Contd) 

GYRO 

THERMAL. 

WARM-UP  DRIFT  TRANSIENTS 

RANDOM  CONSTANT 
EXPONENTIAL  DECAY 

62. 

WT  X 

Tx 

Gyro  Thermal  Warm-Up  Drift  Rate 

0,014°/hr 

63. 

WTy  Y 

Gyro  Thermal  Warm-Up  Drift  Rate 

0.014°/hr 

64. 

WTz  Z 

Gy-'o  Thermal  Warm-Up  Drift  Rate 

0.024°/hr 

ACCELEROMETER  ATTITUDE  DEPENDENT  SCALE  FACTOR 

ERRORS 

RANDOM 

CONSTANTS 

65. 

AKAAx  X 

Accelerometer  Attitude  Dependent 
Factor  Error 

Scale 

10  ppm 

66. 

AKAAy  Y 

Accelerometer  Attitude  Dependent  Scale 
Factor  Error 

10  ppm 

67. 

AKAAz  1 

Accelerometer  Attitude  Dependent 
Factor  Error 

Scale 

30  ppm 

MEASUREMENT  TIMING  ERROR 

EXPONENTIALLY 

CORRELATED 

68. 

4> 

Measurement  Timing  Error 

1 ppm 

F matrix  elements  corresponding  to  the  gyro  quadrature  G drift 
coefficients  and  the  measurement  timing  error  discovered  during  the 
laboratory  tests  are  giver,  next  for  completeness. 
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F(3I,4)  = 
F(32,4)  = 

F(33,4)  = 
F{3I,5)  = 
F(32,5)  = 
F(33,5)  = 
F(33,6)  = 

F (49,68) 
F (50,68) 
F (51,68) 
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The  68x10  noise  distribution  matrix  G is  illustrated  in  Figure  5-1. 
This  matrix  contains  ten  nonzero  elements,  four  of  which  are  located  in 
the  rows  correspond]' no  to  the  long  correlation  time  exponentially 
correlated  random  noise  gyro  drift  rate  and  measurement  error  stochastic 
processes.  These  nonzero  elements  are  readily  identified  from  the 
stochastic  differential  equations,  modeling  these  noise  processes,  given 
in  Section  3.2.?.  The  remaining  six  nonzero  elements  are  located  in  the 
rows  corresponding  to  the  dynamic  gimbal  angle,  31  through  33,  and 
velocity,  49  through  51,  measurement  states.  This  modeling  accounts 
fur  the  uncertainty  of  the  measurement  states  produced  by  the  short 
correlation  time  gyro  and  accelerometer  random  noise  errors.  These  errors 
were  omitted  *-'rom  the  Kalman  Filter  state  vector  due  to  the  filter's 
inability  to  estimate  them.  The  magnitude  of  these  G matrix  elements 
was  determined  s:  ch  that  the  variance  of  the  dynamic  measurements  states 
is  the  same  as  when  the  short  correlation  time  raniom  noise  errors  were 
included  as  stages  in  the  filter. 

The  diste. nance  vector  w(t)  is  a 10x1  input  vector  whose  elements 
are  zet o m. an,  unity  variance  white  noise  processes.  These  white  noise 
processes  are  assumed  to  be  uncorrelated  with  each  othfr  resulting  in 
the  10x10  system  noise  matrix  Q being  an  identity  matrx. 

The  6x68  measurement  matrix  contains  up  to  a maximum  of  six  nonzero 
elements  of  unity  magnitude  in  the  columns  correspond! ng  to  the  rows  of 
the  measurement  states,  if  all  three  velocity  and  gimbal  angle  measure- 
ments are  utilized,  the  measurement  matrix  requires  six  rows.  The 
measurements  eventually  decided  upon  consist  of  the  three-velocity  but 
only  the  azimuth  gimbal  synchro  angle  measurement,  requiring  only  four 
rows  for  the  measurement  matrix. 

The  measurement  z and  measurement,  noise  v vectors  require  from  four 
to  six  states,  depending  on  whether  the  roll  and  pitch  gimbal  synchro 
angle  measurements  are  used.  The  first  one,  or  three,  elements  of  the 
v vector  are  zero  mean  Gaussian  white  random  variables,  with  5-arc-second 
standard  deviations,  and  model  the  differential  gimbal  synchro  angle 
measurement  noise.  The  last  three  elements  have  similar  properties  and  a 
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standard  deviation  representative  of  the  accel eromete**  quantization 
process  (quantization  level/3.46)  and  model  the  artificial  velocity 
measurement  noise.  The  assumption  that  elements  of  the  measurement  noise 
vector  are  statistically  independent  of  each  other  results  in  the  6x6 
measurement  noise  matrix  R being  a diagonal  matrix.  Diagonal  elements  of 
the  R matrix  are  variances  of  the  corresponding  measurement  noise  vector 
elements. 

r^e  66;'.63  covariance  matrix  P is  initially  a diagonal  matrix 
composed  largely  of  the  IMU  error  source  variances.  The  initial  uncertainty 
of  all  the  68  Kalman  Filter  states,  in  terms  of  their  standard  deviations 
as  determined  by  tailing  the  square  root  of  the  corresponding  initial 
covar’ance  matrix  diagonal  elements,  is  given  in  Table  5-1,  The  initial 
uncertainty  of  many  of  tnesc  error  states  is  given  in  Reference  1.  Since 
Reference  1 values  correspond  to  inertial  instrument  biases  and  scale 
factors  which  have  bean  trimmed  for  navigation,  while  the 'Kalman  Filter 
mechanization  does  not,  these  instrument  errors  were  increased  to  their 
untrimmed  values,  in  those  instances,  gyro  spin  axic  mass  unbalance 
coefficients,  where  IMU  laboratory  calibration  test  results  showed  that 
the  actual  errors  were  significantly  larger  than  those  given  in 
Reference  1 the  initial  uncertainties  were  increased  to  reflect  the 
test  results.  T.ht  initial  variances  of  the  velocity  and  differential 
gimbal  angle  measurement  states  are  identically  zero  due  to  tne  manner 
in  which  these  states  were  defined  The  initial  variances  of  the  IMU 
case-to-navigation  frame  misalignment  angles  depend  on  the  laboratory 
test  table  or  mounting  fixture  used  and  the  manner  in  which  the  IMU 
case  is  aligned  to  the  navigation  frames.  This  topic  is  further 
discussed  in  Section  7. 

A partial  covariance  matrix  *'e-initiaIization  is  required  each 
time  the  platform  frame  is  ref  *nted,  for  several  reasons.  The 
variances  of  the  velocity  and  deferential  gimbal  angle  measurement 
states,  having  increased  during  testing  at  a given  p’atform  attitude, 
must  be  reset  to  zero.  The  gimba1  synchro  angle  bias  variances,  having 
decreased  over  the  previous  test  attitude,  must  be  increased  to  the 
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nominal  6-arc -minute  uncertainty  to  reflect  that  new  biases  are  being 
encountered.  Covariance  matrix  elements  associated  with  the  gyro  bias 
and  accel erometer  scale  factor  platform  attitude  dependent  error  model 
residual  states  are  reset  to  reflect  the  attitude  dependency.  In 
addition  to  resetting  the  variances  of  these  states,  any  covariances 
that  have  built  up  between  other  error  states  must  be  re- i ni tial i zed  to 
zero. 


In  addition  to  those  error  sources  included  in  the  Kalman  Filter 
state  vector,  it  is  desirable  to  computer  optimal  estimates  and  estimation 
errors  for  several  composite  error  quantities.  Inertial  navigation 
system  parameters,  such  as  total  gyro  and  platform  drift  and  platform 
and  accelerometer-to-navigation  frame  misalignment  angles,  provide 
additional  insight  to  th~  Kalman  Filter's  capability  to  separate  the 
IMIJ  calibration  parameters  and  estimate  them  accurately.  Optimal 
estimates  for  these  composite  parameters  are  computed  as  functions  of  the 
sum  of  appropriate  state  vector  element's  optimal  estimates.  Optimal 
estimation  errors  for  these  composite  error  quantities  are  computed  by 
use  of  the  expectation  operator  and  corresponding  elements  of  the 
covariance  matrix. 

5.2.3  Integration  Algori thms 

The  majority  of  the  storage  and  computational  time  computer  resources 
required  to  implement  and  exercise  the  Kalman  Filter  algorithm  on  a 
digital  computer  are  associated  with  propagating  the  Kalman  Filter  state 
vector  and  covariance  matrix.  Computational  time  requirements  can  easily 
become  burdensome,  when  we  consider  that  the  state  vector  and  covariance 
matrix  are  propagated  for  time  periods  of  one  to  two  hours.  The  necessary 
computer  resources  can,  however,  be  kept  within  reasonable  limits  by 
carefully  selecting  the  integration  algorithms  required  to  propagate 
these  quantities. 

The  final  version  of  the  .Kb  calibration  Kalman  Filter  was  hosted 
on  a CDC -6600  digital  compute"  bav’ng  a "word  length  of  60  bits.  The  state 
vector  was  propagated  with  laboratory  test  data  stored  on  a magnetic. 
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tape  correspond? ruj  to  off-line  non-real-time  operation.  Total  !MU  test 
time  was  69  minutes  with  measurement:,  being  nominally  taken  every 
"'C  seconds.  CDC-6600  computer  resources  required  to  execute  the  job 
were  900  seconds  for  compilation  and  execution  time  and  221,000  octal 
words  for  memory. 

In  developing  the  CDC-6600  Kalman  Filter  computer  program,  it  was 
found  that  the  most  efficient  operation  of  acceptable  accuracy  was 
accomplished  by  using  relatively  large  integration  step  sizes  accompanied 
by  sophisticated  integration  algorithms.  The  covariance  matrix  is 
propagated  in  one  step  between  the  measurements,  integration  step  size 
of  30  seconds,  according  to  the  fifth  order  Kutta-Mersion  integration 
algori thm. 

1 = P(t) 

o 

».  3 ». + *(».) 

>*=’<,-  -g~  Ow + pw] 

»3  3 »o  + -81-  [P(*c)  + 3P(»*)] 

\ 3 >»  + -f1-  PW ' PW  * 2AT  p(s) 

p(,  + At)  = , + -£-[*(,,)  + P(5)]  * f(»s) 

where  AT  is  the  integration  step  size 

T y 

P{ »)  = t)  P<f)  * [r(t)P(t)]  + 6(t)Q(t)6(t) 

The  covariance  matrix  derivative  P(t)  is  determined  using  the  above 
expression.  Due  to  tie  sparse  nature  of  the  system  dynamics  matrix  F(t), 
the  P(t)  algorithm  was  coded  in  a manner  that  eliminated  the  bulk  of 
the  zero  multiples. 
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Propagation  of  the  Kalman  Filter  state  vector  x(t)  requires 
substantially  less  computer  resources  than  those  necessary  for  propagation 
of  the  covariance  matrix.  Modified  Eule"'  integration  is  used  to 
propagate  the  state  vector  with  an  integration  step  size  of  one  second. 

* (t  + At)  = x (t)  + (0  + *(»  + At)] 

where 

»(»)  - F(')t(0 
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SECTION  VI 

THEORETICAL  COMPUTER  SIMULATION  RESULTS 

5.1  SIMULATION  OBJECTIVES 

A CPC -6600  computer  simulation  of  the  KT-73  IMU,  Kalman  Filter  and 
laboratory  calibration  test  was  developed  in  consideration  of  several 
major  objectives 

(1)  To  validate  the  concept  of  tie  precise  postulated  for  separation 
o:‘  the  IMU  calibration  parameters. 

/ 

(2)  To  determine  the  accuracy  to  whicn  the  Kalman  Filter  can 
estimate  the  IMU  calibration  parameters,  considering  the  effects  of 
inertial  instrument  random  disturbances,  measurement  noise,  and  limited 
test  time. 

(3)  To  perform  limited  analysis  in  regard  to  the  power  of  the  iMlJ 
calibration  methodolc jv  and  tne  usefulness  of  the  Kalman  Filter  as  a tool 
for  conducting  IMU  error  model  development  and  validation  studies. 

(4)  To  justify  the  first-oi der  equation  linearizations  made  for 
facilitating  a solution  to  the  calibration  equations. 

(5)  To  assist  in  the  preparation  of  a computer  software  specification 
and  test  procedure  for  the  ensuing  IMU  calibration  laboratory  development 
program. 

5 ?.  VALIDATION  OF  CALIBRATION  METHODOLOGY  AND  ASSUMPTIONS 

To  obtain  the  highest  level  of  confidence  possible  by  the  use  of 
computer  simulations,  the  concept  validation  philosophy  of  demonstrating 
strict  convergence  of  the  Kalman  Filter  optimal  estimates  to  the  error’ 
sources  of  a simulated  KT-/3  1 HU  was  adopted.  Strict  sense  convergence 
imp'ie.  that  I he  ensemble  average  of  the  simulat'd  estimation  errors 
conform  ar  the  zero  mean  gaussian  processes  having  time-dependent 
variances  given  by  the  filter's  covariance  matrix. 
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In  keeping  with  the  concept  validation  philosophy,  the  developed 
KT-73  I MU  simulation  contains  all  the  error  sources  given  in  Reference  1. 
The  I MU  simulation  is  based  on  the  whole  number,  nonlinear  equations 
develooed  above,  and  utilizes  random  number  generators  to  simulate  the 
random  bias  and  stochastic  process  error  sources.  Digital  gyro  torquing 
ar.d  accelerometer  quantization  effects  are  included.  The  adequacy  to 
which  the  Reference  1 error  sources  model  the  real  KT-73  I MU  will  be 
addressed  in  the  laboratory  development  program  and  also  in  conducting 
the  aided-inertial  navigation  software  development  and  evaluation  tests 
in  1 abl e 1 -1 . 

Figure  6-1  illustrates  the  Kalman  Filter's  ability  to  accurately 
estimate  I MU  platform  drift  rate  and  pseudo  velocity  errors  in  terms  of  the 
standard  deviation  of  the  estimation  error.  Platform  drift  rate  is  a 
composite  error  source  containing  up  to  eighteen  error  sources  per  axis 
for  the  platform  and  navigation  frames  approximately  aligned  situation 
shown  in  the  figure.  Similar  remarks  also  apply  to  the  pseudo  velocity 
errors.  Estimation  error  standard  deviations  for  the  psuedo  velocity 
errors  are  obtained  by  taking  the  square  root  of  the  appropriate  diagonal 
element  of  the  filter's  covariance  matrix.  Since  total  platform  drift 
rate  error  is  not  a state  variable  in  the  Kalman  Filter,  corresponding 
estimation  accuracies  had  to  be  obtained  by  utilizing  the  expectation 
operator  in  conjunction  with  on  and  off  diagonal  elements  of  the 
covariance  matrix. 

Each  graph  in  Figure  6- i contains  two  sets  of  data,  correspondi ng  to 
whether  or  not  the  system-level  IMU  velocity  and  giinbal  angle  measurements 
are  taken.  Velocity  measurement  errors  corresponding  to  the  low-gain 
accelerometer  mode  quantization  levels  were  assumed.  The  x,  y,  and  z 
platform  axes  symbols  indicate  the  time  instants  the  six  simultaneous 
IMU  measurements  were  taken  and  are  uniformly  spaced  at  thirty  second 
intervals.  The  line  segments  between  measurement  points  indicate  how 
the  estimation  uncertainty  propagates  in  time.  The  primary  information 
conveyed  by  these  graphs  is  that  the  Kalman  Filter  is  capable  of 
estimating  composite  system- level  IMU  errors  to  an  accuracy  essentially 
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limited  by  the  inertial  instrument  random  disturbances  in  a five  to  ten 
minute  I MU  test.  The  length  of  time  required  to  achieve  this  level  of 
accuracy  is  heavily  influenced  by  the  magnitude  of  the  measurement  errors. 

Although  the  Kalman  Filter  can  accurately  estimate  the  composite 
error  quantities  during  a single  platform  orientation  test,  it  cannot 
separate  or  accurately  observe  the  individual  I MU  calibration  parameters. 
The  filter  can,  however,  during  the  single  platform  orientation  test, 
heavily  correlate  certain  groups  of  error  sources  from  an  initially 
uncorrelated  condition.  To  calibrate  the  1MU,  the  error  sources  must  be 
decorrelated  and  accurately  observed.  The  means  by  which  this  is 
accomplished  is  by  testing  the  I MU  with  its  platform  at  numerous 
orientations  with  respect  to  the  earth  rate  and  gravity  vectors  as 
discussed  in  Section  IV. 

Figure  6-2  is  a graphic  illustration  of  the  error  separation  process. 
This  figure  shows  an  increasing  Kalman  Filter  capability  to  estimate 
the  y-gyro  bias  drift  rate  error  component  as  additional  platform  test 
attitudes  are  introduced.  Platform  attitudes  No.  1 through  5,  10,  and  11 
correspond  to  those  of  Table  4-1.  Platform  attitudes-No.  7 through  9 
are  the  additional  attitudes  introduced  to  permit  error  source  separation 
in  the  presence  of  the  west-gyro-commanded  angular  velocity  error  due  to 
the  azimuth  platform-to-navigation  frame  misalignment.  Attitudes  No.  6 
and  12  are  identical  to  No.  1 and  illustrate  that  repeated  use  of  the 
same  test  attitude  does  not  significantly  improve  error  source  observation 
accuracy.  These  attitudes  are  synonomous  with  the  navigation  frame  axes 
of  north,  west,  and  up. 

The  IMU  is  tested  for  five  minutes  per  attitude  as  suggested  by  the 
estimation  accuracy  versus  measurement  time  characteristic  of  Figure  6-1. 
Total  IMU  test  time  would,  however,  be  somewhat  longer  than  the  one  hour 
measurement  time  shown  due  to  the  several  minutes'  platform  slew  time 
required  to  periodically  re-orient  the  platform.  The  sequence  of  platform 
test  attitudes  shown  was  selected  such  that  the  most  significant  IMU  error 
sources  are  calibrated  first  and  allow  the  option  of  omitting  latter  test 
attitudes  in  favor  of  reducing  test  time. 
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The  upper  and  lower  bounds  superimposed  about  the  simulated  y-gyro 
bias  error  represent  three  standard  deviations  of  the  optimal  estimation 
error  stochastic  process.  The  standard  deviations  used  in  determining 
these  bounds  were  obtained  by  taking  the  square  root  of  covariance  matrix 
diagonal  element  corresponding  to  the  y-gyro  bias  error  as  a function  of 
time.  Optimal  estimation  theory  predicts  that  the  Kalman  Filter's 
optimal  estimates  should  converge  to  the  simulated  y-gyro  bias  error  to 
within  these  bounds  a certain  percentage  of  the  time  on  an  ensemble 
average  basis.  For  Gaussian  processes  and  three  standard  deviation 
bounds,  this  percentage  is  99.5%  of  the  time.  Inspection  of  Figure  6-2 
shows  that  all  120  or  100%  of  the  filter's  optimal  estimates  fall  within 
the  three  standard  deviation  bounds,  thus  satisfying  the  criterion. 
Examination  of  the  computer  simulation  results  for  convergence  of  the 
filter's  optimal  estimates  to  the  other  simulated  I MU  error  sources 
showed  that  the  criterion  was  met  also. 

Several  conclusions  can  be  reached  as  a result  of  the  filter's 
optimal  estimates  converging  to  the  simulated  IMU  error  sources  to 
within  the  accuracy  predicted  by  optimal  estimation  theory.  The  major 
conclusion  is  that  the  gimballed  IMU  calibration  problem,  as  formulated, 
can  be  adequately  represented  b>  a system  of  first-order  linear 
equations.  A second  conclusion  is  that  treatment  of  velocity  quantization 
errors  as  measurement  noise,  omission  of  digital  gyro  pulse  torquing 
effects,  and  other  simplifying  assumptions  are  justifiable.  A final 
conclusion  is  that  the  Kalman  Filter  is  a suitable  analytical  technique 
for  calibrating  IMUs. 

6.3  IMU  CALIBRATION  ACCURACY  LIMITATIONS 

Examination  of  the  estimation  accuracy  bounds  depicted  in  Figure  6-2 
shows  that  the  estimation  accuracy  increases  rapidly  for  the  first  four 
or  five  platform  test  attitudes.  A slight  increase  in  estimation  accuracy 
occurs  during  the  latter  platform  test  attitudes  due  to  the  further 
separation  of  the  y-gyro  bias  error  from  the  other  IMU  calibration 
parameters  modeled  as  random  biases.  The  y-gyro  bias  estimation  error 
does  not,  however,  approach  zero  at  the  conclusion  of  the  platform  test 
attitude  sequence. 
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To  determine  the  reasons  for  the  limitation  in  I MU  calibration 
accuracy,  the  analytical  results  shown  in  Figure  6-3  were  computed  by 
use  of  the  Kalman  Filter  simulations.  The  upper  graph  shows  the  optimal 
estimation  errors  for  all  three  gyro  bias  drift  rate  components.  This 
representation  contains  the  same  estimation  accuracy  information  as 
Figure  6-2  but  is  plotted  as  one  standard  deviation  on  an  absolute  scale. 
This  data  is  obtained  by  taking  the  square  root  of  the  appropriate 
diagonal  element  of  the  updated  covariance  matrix. 

The  lower  graph  in  Figure  6-3  shows  the  estimation  accuracy  for  the 
random  variables  defined  as  the  sum  of  the  gyro  bias  and  long  correlation 
time  random  gyro  drift  rate  error  components.  Comparison  of  the  upper 
and  lower  graphs  illustrate  the  gyro  bias  calibration  accuracy 
limitations  imposed  by  the  gyro  random  noise  components.  Each  gyro  is 
modeled  with  two  exponentially  correlated  disturbances  having  1 and 
3600  second  correlation  times  and  0.003  and  0.005  degree  per  hour  drift 
rates,  respectively.  Reference  to  the  upper  graph  shows  that  the  long 
correlation  time  disturbance  imposes  itself  as  an  approximate  lowe1"  bound 
on  gyro  bias  calibration  accuracy.  The  reason  that  the  long  correlation 
time  random  gyro  disturbance  acts  as  lower  bound  for  the  gyro  bias  error 
estimation  accuracy  is  that  time  is  the  only  difference  in  the  functional 
dependency  of  these  two  error  sources.  Accurate  measurements  taken 
frequently  over  an  approximate  correlation  time  period  are  required  to 
separate  them.  This  condition  is  not  satisfied  during  the  calibration 
test  due  to  those  platform  test  attitudes  where  platform  drift  can  only 
be  inaccurately  inferred  from  the  differential  gimbal  angle  measurements. 
This  limitation  is  illustrated  in  the  lower  graph  in  Figure  6-3  by  the 
filter's  ability  to  estimate  the  sum  of  these  two  error  sources  to  an 
accuracy  now  limited  by  the  short  correlation  time  disturbance. 

Estimation  of  the  I Mil  calibration  parameters  to  an  accuracy  limited 
by  the  inertial  instrument  random  disturbances  is  the  desired  result 
from  the  navigational  accuracy  standpoint  Conceptually,  this  result 
implies  that  the  I MU  calibration  error  sources  could  be  software 
compensated  for  with  the  pure  inertial  navigation  accuracy  achievable 
oith  the  I MU  being  essentially  limited  by  the  inertial  instrument  random 
disturbances . 
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CYRO  BIAS  CALIBRATION 
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figure  6-3.  Gyro  Bias  and  Gyro  Bias  Plus  Random  Noise 
Calibration  Accuracy 
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6.4  DYNAMIC  ERROR  ESTIMATION 

The  Kalman  Filter's  ability  to  estimate  a dynamic  composite  1 MU  error 
as  opposed  to  the  static  I MU  calibration  parameters  is  illustrated  by 
the  total  gyro  drift  rate  characteristic  shown  in  Figure  6-4.  Total 
gyro  drift  is  defined  as  a random  variable  composed  of  the  sum  of  the 
long  and  short  correlation  time  disturbances,  thermal  warm-up,  constant 
and  attitude-dependent  biases,  input  :.nJ  spin  axis  mass  unbalance,  and 
anisoelastic  drift  rate  error  components.  Total  platform  drift  would 
additionally  contain  the  gyro  scale  factor,  gyro-to-platform  misalignment 
angle,  and  commanded  angular  velocity  error  sources. 

The  simulated  x-gyro  drift  rate  is  plotted  with  straight  line  segments 
connecting  the  30-second  measurement  intervals  and  does  not  show  the 
fine  grain  structure  due  to  the  random  disturbances.  Convergence  of  ttie 
Kalman  Filter's  optimal  estimates  to  the  simulated  gyro  drift  rate  error 
is  illustrated  along  with  the  one  standard  deviation  estimation  error 
upper  and  lower  Dounds.  Assuming  gaussian  statistics,  68‘T  of  the  filter's 
optimal  estimates  should  fall  within  these  boundaries.  Optimal  estimation 
errors  appear  consistent  with  this  criterion. 

6.5  IMU  ERROR  MODEL  REFINEMENTS 

A major  objective  of  the  IMU  calibration  study  is  to  develop  a 
methodology  and  companion  analytical  tool  that  affords  some  capability 
for  conducting  IMU  error  model  development  and  validation  efforts.  An 
essential  feature  on  any  such  approach  is  the  ability  to  start  with  a 
functional ly  correct  IMU  error  model,  but  with  incorrect  statistical 
information  regarding  the  magnitudes  of  the  individual  error  sources, 
and  through  iterative  testing  determine  the  correct  error  source 
magnitudes.  The  problem  of  refining  a functionally  incorrect  IMU  error 
model  is  addressed  in  Section  8. 

To  address  the  problem  of  refining  t^e  a priori  IMU-error-mode I 
statistical  magnitude  data,  the  analysis  depicted  in  figure  6-5  was 
conducted.  The  upper  graph  shows  that  the  nominal  computer  loaded 
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Figure  6-6.  Gyro  Scale  Factor  Frtunation 
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gyro  scale  factor  errors  do  not  contribute  significant  platform  drift 
rate  components.  Since  these  error  sources  do  not  contribute  significant 
platform  drift  rate  components  relative  to  the  gyro  disturbances,  the 
filter  is  precluded  from  reducing  the  uncertainty  associated  with  them. 

Gyro  scale  factor  errors  are  shown  in  terms  of  the  platform  drift  rate 
error  components  they  produce  to  emphasize  the  calibration  accuracy 
limitations  imposed  by  the  gyro  random  noise  errors.  The  lower  graph 
shows  the  filter's  performance  against  a simulated  IMU  containing  a 
x-gyro  scale  factor  error  an  order-of-magni tude  larger  than  assumed 
in  the  initial  covariance  matrix. 

The  optimal  estimates  for  two  simulation  runs  are  plotted  in  the  lower 
graph.  The  plus  symbols  correspond  to  the  first  run  and  an  initial 
covariance  matrix  based  on  the  nominal  a priori  error  model  data.  The 
nominal  x-gyro  scale  factor  error  assumed  for  the  initial  covariance 
matrix  is  a factor  of  ten  less  than  that  simulated  in  the  IMU.  The  filter 
begins  to  estimate  an  abnormally  large  error,  optimal  estimates  compared 
with  the  estimation  error  standard  deviation,  during  the  fourth  platform 
attitude  test.  It  is  interesting  to  note  that  the  Kalman  Filter  begins 
to  estimate  the  abnormally  large  scale  factor  error  during  the  same 
platform  test  attitude  that  the  e^ror  source  separation  data  of  Table  4-1 
predicts.  The  filter's  optimal  estimates  do  not,  however,  strictly 
converge  to  the  simulated  error  in  accordance  with  the  three  sigma 
bounds.  This  failure  of  ihe  filter's  optimal  estimates  to  converge  to 
the  simulated  error  is  due  to  the  non-conformance  of  the  initial  covariance 
matrix  with  the  simulated  error. 

A major  consideration  at  the  end  of  the  first  simulation  run  is  the 
behavior  of  the  Kalman  Filter's  estimates  for  all  the  other  IMU  calibration 
parameters.  If  the  filter  estimates  abnormally  large  errors  for  some  of 
the  other  calibration  parameters,  this  might  preclude  isolation  of  ^he 
problem  to  its  true  source.  Examination  of  filter's  estimates  for  all 
other  error  sources  indicates,  however,  that  these  estimates  converge 
to  the  simulated  errors  to  within  the  accuracy  predicted  by  the  filter's 
ronvariance  matrix. 
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The  second  simulation  run  was  made  to  verify  that  filter  convergence 
rould  be  obtained  by  increasing  th°  initial  covariance  matrix  in 
accordance  with  the  x-oyro  scale  factor  optimal  estimates  obtained  from 
the  first  run.  Specifically!  the  initial  x-gyro  scale  factor  variance 
was  increased  ^rom  that  corresponding  to  a platform  drift  rate  component 
of  0.0035  to  0.024  degrees  per  hour  as  indicated  by  the  filter's  final 
first  run  estimate.  The  small  circle  symbols  represent  the  second  run 
optimal  estimates  and  show  a convergence  which  is  consistent  with  the 
estimation  uncertainty  obtained  from  the  first  run.  This  result 
indicates  that  calibration  accuracy  is  essentially  determined  by  the 
inertial  instrument  random  disturbances  and  is  relatively  insensitive 
to  the  initial  covariance  matrix. 
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SECTION  VI) 

IMU  LABORATORY  TEST  REQUIREMENTS 
7.1  SOFTWARE  ALTERNATIVES 

Digital  computer  computational,  storage,  and  associated  programming 
requirements  for  implementing  the  IMU  calibration  tests  and  estimating 
the  error  sources  can  vary  greatly  in  accordance  with  the  specific 
application  and  software  approach.  All  software  approaches  must,  however, 
satisfy  three  functional  requirements. 

(1)  The  platform  must  be  commanded  tc  precess  at  earth  rate  during 
each  IMU  test  sequence. 

(2)  System  level  IMU  psaudo  velocity  and  gimbal  angle  measurements 
must  be  taken  and  recorded. 

• (3)  IMU  error  source  estimates  must  oe  computed. 

The  first  two  requirements  involve  online  in  real-time  s iftware 
and  can  be  met  in  a similar  manner  for  various  approaches  with  software 
complexity  beiny  a function  of  the  desired  automation.  The  manner  of 
accomplishing  the  third  function  can  be  either  on-line  in  real-time  or 
off-line  af ter-the-fac*  and  greatly  distinguishes  the  software  approaches. 

7.1.1  Kalman  Filter  Software  Approach 

Two  approaches  to  the  Kalman  Filter  error  source  estimation  software 
were  considered  for  implementation  of  the  laboratory  IMU  tests.  The 
first  approach  consisted  of  ari  off-line  parameter  estimation  scheme  in 
which  the  error  source  estimates  were  computed  after  the  IMU  test  was 
completed.  The  Kalman  Filter  mechanization  and  IMU  test  procedure 
formulated  purposely  omitted  rea'i-time  closed-loop  platfmn  control  in 
favor  of  the  off-line  software  flexibility.  In  this  approach,  the  IMU 
test  wou!d  be  conducted  with  use  of  a small  general  purpose  digital 
computer/,  $)  with  the  IMU  measurements  being  recorded  on  a magnetic  lope. 


7? 


AFAL-TR- 77-75 


\ 

The  recorded  I MU  measurements  would  serve  as  input  data  to  the  Kalman 
Filter  software  program  hosted  on  a large  scientific  digital  computer. 
This  approach  permits  use  of  the  Kalman  Filter  software  program  developed 
for  the  CDC-6600  computer  simulation  studies. 

The  second  Kalman  Filter  software  approach  considered  utilized  a 
small  general  purpose  digital  computer  and  stored  Kalman  gain  matrices 
to  estimace  the  error  sources  on-line  in  real-time.  This  technique  is 
possible  because  trie  filter's  covariance  matrix  can  be  propagated  and 
updated  ^independently  of  the  1M‘J  measurements , This  allows  prior 
computation  of  the  covariance  and  related  Kalman  gain  matrices  on  a 
large  digital  computer.  The  precomputed  Kalman  gain  matrices,  corre- 
sponding to  specified  I MU  measurement  times,  would  be  stored  on  a 
magnetic  tape  and  read  by  the  small  tost  computer  as  required  during 
the  test.  Estimation  of  the  error  sources  is  well  within  a small 
digital  computer's  capabilities  and  is  accomplished  by  multiplying  the 
measurement  vector  by  the  Kalman  gain  matrix  after  each  IMU  measurement 
and  propagating  the  optimal  estimates  between  measurements.  The 
software  program  could  be  written  in  assembly  or  a higher  order  language 
depending  on  the  support  software  available  for  the  small  digital 
computer. 

The  first  approach  of  off-line  error  sources  estimation  was  selected 
for  the  initial  IMU  calibration  tests  for  several  reasons.  Real-time 
software  development  was  required  only  for  recording  data  and  not 
additionally  for  error  source  estimation  as  in  the  second  approach. 

The  ease  of  modifying  the  CDC-6600  Fortran  !V  Kalman  Filter  computer 
program  and  the  ability  to  exercise  the  modified  versions  against  the 
mag  tape  stored  IMU  measurement  data  without  continuously  repeating  the 
iMU  tests  was  expected  to  greatly  facilitate  the  development  effort. 

This  consideration  was  especially  important  i.i  view  of  the  IMU  error 
model  development  anJ  validation  objective  of  the  study. 

The  major  disadvantage  of  the  selected  approach  is  the  turn-around 
time  associated  with  the  CDC-6660  computer  system.  It  was  thought  that 
if  aided- inertial  navigation  software  development  efforts  demonstrated 
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a need  for  having  I MU  calibration  test  results  prior  to  conducting  the 
software  development  tests,  the  real-time  error  estimation  scheme  could  be 
implemented  for  the  production  calibration  tests.  This  need  has  not, 
however,  materialized. 

7.1.2  Real  Time  Test  Software  Approach 

To  keep  the  real-time  software  development  effort  to  a minimum, 
maximum  utilization  was  made  of  the  Si nger-Kearfott  AN/ASN-90  Inertial 
Measurement  Set  and  associated  equipment,  including  the  SKC-2000  digital 
computer  and  the  operational  computer  program  (OCR)  documented  in 
Reference  4.  A functional  block  diagram  of  the  hardware  used  in  conducting 
the  IMU  calibration  tests  is  given  in  Figure  7-1.  The  Singer  equipment 
constitutes  a self-contained  inertial  navigation  capability  including 
IMU  ground  alignment  and  a two-position  north  and  west  gyro  auto- 
cal ibration. 

The  philosophy  for  obtaining  the  necessary  real-time  test  software  was 
to  utilize  applicable  portions  of  the  Singer  OCP  with  minor  software 
modifications  being  made  as  required.  References  to  Singer's  OCP 
development  specification  indicated  that  this  approach  was  feasible 
and  would  entail  a modification  of  the  ground  align  and  gyro  auto- 
calibration software  routines. 

Two  versions  of  the  OCP  real-time  test  software  were  developed. 

The  first  version  retained  the  15  sequence  ground  align  and  gyro  auto-cal 
software  routines  so  that:  (1)  The  IMU  case  could  be  properly  aligned 
to  the  navigation  frame,  and  (2)  The  gyro  turn-on  biases  recovered  by 
the  OCP  could  be  compared  to  those  estimated  by  the  Kalman  Filter. 
Compatibility  between  the  OCP  and  Kalman  Filter  gyre  turn-on  bias 
estimates  was  ensured  by  conducting  the  12  IMU  calibration  platform 
attitude  test  sequences  immediately  after  the  15  OCP  ground  align  and 
gyro  auto-cal  test  sequences.  Comparison  of  the  OCP-recovered  gyro 
biases,  as  well  as  the  high  and  low  gain  accelerometer  bias  and  scale 
factor  and  gyro  scale  factor  ASM-375-test-set  calibrated  errors,  with 
those  recovered  by  the  Kalman  Filter,  was  considered  essential  in  the 
development  of  the  calibration  technique. 
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After  the  development  of  the  calibration  technique  was  well  underway, 
the  second  version  of  the  real-time  test  software  was  developed  to 
eliminate  a scheduling  problem  inherent  in  the  above  versions.  The 
15-sequence  ground  align  and  gyro  auto-cal  routines  required  24  minutes 
to  execute  and  thus  prohibited  obtaining  the  I MU  measurements  necessary 
to  estimate  the  inertial  instrument  warm-up  transients  corresponding  to 
the  A-7D  aircraft  fine  ground  alignm-;nt  time  interval  of  7 to  10  minutes 
after  power  turn-on.  The  second  version  was  similar  . j the  first,  except 
that  the  15-sequence  ground  align  and  gyro  auto-cal  routines  were 
omitted.  Proper  IMU  case  alignment  to  the  navigation  frame  was  obtained 
by  use  of  the  original  OCP  if  realignment  was  made  necessary  due  to 
dismounting  and  remounting  the  IMU  on  the  tilt  table,  as  discussed  in 
Section  7.2.1. 

Initiating  arid  terminating  the  different  IMU  platform  test  sequence^ 
and  recording  the  IMU  measurements  on  a magnetic  tape  could  not  readily 
be  accomplished  with  the  Singer  AN/ASN-90  and  associated  equipment. 
Referring  to  Figure  7-1,  the  teletyoe  and  PDP-11/45  computers  wore 
added  for  these  purposes.  The  PDP-1 1/45  computer  and  its  peripheral 
equipment,  is  used  largely  for  bookkeeping  purposes.  Nominal  earth  rate 
torquing  values  are  stored  and  sent  to  the  SKC-2000  computer  at  the  start 
of  each  new  platform  test  sequence.  IMU  measurements  are  transferred 
from  the  SKC-2000  computer  and  recorded  on  magnetic  tape.  Although  the 
initiation  of  new  platform  test  sequences  is  controlled  from  the  teletype, 
the  exact  instants  of  time  the  IMU  measurements  are  taken  is  determined 
by  the  SKC-2000  computer's  internal  clock  and  associated  real-time 
software. 

Input/ output  data  precision  and  IMU  measurement  timing  accuracy 
requirements  were  determined  such  that  errors  in  these  quantities  do 
not  contribute  any  significant  IMU  calibration  errors.  Thi  criterion 
implies  that  the  errors  are  on  the  same  order  as  the  accelerometer  and 
gyro  random  disturbances.  For  gyro  torquing,  the  net  gyro  torquing  PRF 
required  to  command  the  platform  at  earth  rate  must  be  entered  into  the 
SKC-2000  computer  with  an  accuracy  of  50  ppm.  The  magnitude  of  the  time 
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interval  between  I MU  pseudo  velocity  measurements  is  not  critical; 
nominally,  both  10  and  30  seconds  are  used.  This  time  interval  must 
be  constant  to  simplify  execution  of  the  Kalman  Filter  algorithm  and 
precisely  known  to  0.1  msec.  The  total  number  of  pseudo  velocity  pulses 
generated  each  measurement  time  interval  are  recorded  with  the  least 
significant  bit  representing  0.0322  and  0.000322  ft/sec  for  the  low  and 
high  gain  accelerometer  scale  factor  modes,  respectively.  Gimbal  angle 
measurements  are  discussed  in  Section  7.2.2,  but  are  not  automatically 
recorded  via  the  real-time  software. 

7.1.3  Recorded  Pseudo  Velocity  Pulses 

The  finalized  1 MU  platform  test  attitudes  used  in  conjunction  with 
the  second  version  of  the  0CP  real-time  software  are  shown  in  Figure  7-2. 

A sample  of  the  total  I MU  pseudo  velocity  pulses  recorded  by  the  real-time 
software  during  a calibration  test  is  illustrated  in  Figure  7-3.  in  this 
figure,  the  left  hand  column  represents  time  in  seconds  and  from  left  to 
right  the  three  columns  per  test  attitude  represent  the  total  number 
of  accelerometer  pulses  generated  per  measurement  time  interval. 

Data  recording  for  IMU  platform  test  attitude  No.  1 begins  seven 
minutes  after  power  turn-on  and  continues  for  ten  minutes,  with  velocity 
pulses  being  recorded  every  10  seconds.  The.  x and  y level  axis 
accelerometers  are  employed  in  the  high  gain  mode  (the  z accelerometer 
does  not  have  a high  gain  mode),  in  an  attempt  to  accurately  estimate 
the  gyro  and  accelerometer  warm-up  transients. 

Data  recording  periods  are  nominally  five  minutes  long  for  platform 
test  attitudes  No.  2 through  12,  with  measurements  being  recorded  every 
30  seconds.  The  low  gain  navigation  accelerometer  scale  factor  mode  is 
used  for  the  first  five  minutes  of  testing  Tor  each  of  these  test 
sequences.  Accel eromecer  sea1?  factors  are  switched  from  the  low  to  the 
high  gain  mode  at  the  end  of  the  five  minutes  for  test  sequences  No.  3 and 
10,  with  testing  being  continued  for  an  additional  one  and  a half  minutes. 
This  approach  permits  the  Kalman  Filter  to  accurately  observe  the  high 
gain  < celerometer  scale  factor  and  bias  errors  without  adding  two 
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Figure  7-2.  I MU  Platform  Test  Attitudes  at  ■!  C or r^spondi ng 
Gimbal  Angle  Rotation;. 
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additional  test  sequences  specifically  for  this  purpose.  Platform  tilts 
are  limited  to  approximately  six  degrees  to  prevent  saturation  of  the 
capacitive  reset  integrators  (10,000  pulses  per  second  maximum)  when  in 
the  high  gain  mode.  As  can  be  seen  from  Figure  7-3,  I MU  calibration  does 
not  require  an  abnormally  large  quantity  of  measurements  to  be  taken. 

7.1.4  Data  Recording  Problems 

Two  data  recording  problems  were  encountered  during  the  IMU 
calibration  testing.  Eventually  one  problem  was  attributed  to  the  IMU 
hardware  and  the  other  to  the  OCP  r 1-time  test  software,  both  of  these 
problems  caused  improper  operation  of  the  Kalman  Filter  and  were  only 
correctly  isolated  as  data  recording  problems  after  scanning  the  raw 
data  lists  as  depicted  by  Figure  7-3.  The  recorded  velocity  pulse  raw 
data  computer  printouts  proved  to  be  extremely  valuable  during  the 
development  cycle,  particularly  when  the  data  was  massaged  to  show  the 
number  of  accelerometer  pulses  generated  in  each  measurement  interval, 
as  in  Figure  7-3,  rather  than  the  cumulative  time  total  which  constitutes 
the  Kalman  Til  to*-  measurements  . 

The  data  recording  anomaly  attributed  to  the  IMU  hardware  was  termed 
an  accelerometer  dropout  problem.  Anomalous  data  was  found  to  exist  on 
2 of  the  total  34  IMU  calibration  tests  conducted  and  took  trie  form  of  a 
zero  accelerometer  pulse  count  for  one  to  two  minutes  of  the  total 
68  minutes  IMU  measurement  time.  Only  one  dropout  time  interval  was 
found  to  occur  for  each  of  the  two  ca  ibration  tests  that  the  problem 
manifested  itself.  Although  this  aromalous  behavior  was  never  verif-ied 
as  being  caused  by  the  IMU  hardware,  discussions  with  personnel  at  the 
Aerospace  Guidan"1  and  Metrology  Center  (AGMC)  indicated  they  had  observed 
similar  problems  with  the  KT-73  IMUs.  AGMC  personnel  thought  the  problem 
was  due  to  burnt  or  corroded  slip  rings.  Valid  Kalman  Filtei  estimates 
for  the  IMU  calibration  parameter  could  still  be  obtained  by  smoothing 
over  the  zero  pulse  count  anomalous  recordings  with  interpolated  data. 
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The  second  delta  recording  problem  wa  , traced  rc  .•  t r--;n.:  i n-,  tne ; ' : * 
in  the  real-time  OCR  software.  Instead  of  t i ! ,n.:  • in-  s <v  ■ r • no 

pseudo  velocity  pulses  to  the  requ : red  accuracy  of  -1.1  -.sec , errors  as 
large  as  i5  msec,  were  found  to  exist.  The  effect  of  these  timing  errors 
on  I MU  calibration  accuracy  is  discussed  in  Section  8. 

Th°  timing  instability  problem  was  also  isolated  by  scanning  the 
raw  velocity  pulse  data  as  depicted  in  Figure  7-3.  For  the  accelet  omefo> 
random  noise  magnita  •-»  given  by  the  i Mu  error  model  of  Reference  1 and  a 
timing  accuracy  of  0.1  msec,  the  difference  in  the  number  of  velocity 
pulses  generated  over  a measurement  interval  should  never  exceed  one 
when  an  a :< . elereneter  is  oriented  parallel  to  the  gravity  vector  The 
difference  in  velocity  pulse  counts  was  found  to  be  as  high  as  lb  pul_.es 
for  adjacent  measurement,  intervals.  For  a low  gain  accelerometer-  scale 
factor  of  1,000  pulses  per  second  pe1  g,  15  pulses  implies  a 15  msec 
iiieasurc,v''rt  timing  error.  Longer  term  timing  variations  of  greater 
than  ib  ulus  were  also  found  to  exist  over  che  entire  5-minute  test 
i ntc i 1 ' a I 

; he  cause  of  the  timing  errors  was  tiaced  to  tin-  "interrupt.' 
prior  itie.  of,  and  tire  fa 'lure  to  remove  all  .mins  es'.urv  Software 
■on'  ties  'com,  trie  Oi.P  roui-time  software  prog-am.  T 'it-  timing 
ins  tuM  r’.y  problem  was  never  completely  solved  but  the  ms  tub'  i i ties 
.-re  » edo  ed  to  roughly  one  msec  in  magn 1 fade . s i ru.  e a timing  error 
of  thi,  magnitude  is  not  negligible  relative  to  the  inertial  instrument 
random  di  si urbaru.es , mud i I i ations  si  the  Kalman  fr'fer  had  to  be  made 
*o  .Huiu'i  tor  its  presence.  The  short  term  uncon  elated  ti"iing  in- 
stability ,1a  a • i-Oi.'ii tea  for  by  adding  a disturbance  to  the  affected  or 
pseudo  ve . oi.  i f.y  *rror  states.  The  longer  term  .„un  luted  , irr.ing  tnstabi  : i 
was  aeeouiited  for  in  a first  order  manner  by  adding  an  additional  clock 
(useudo  accelerometer  scale  ‘‘actor)  ern  r to  t ne  Kalman  filter  state 
vector,  trror  modeling  of  the  clock  error  state  is  characterized  by  a 
standard  deviation  of  one  ppm  and  a correlation  time  of  ?50  seconds.  The 
net  result  of  these  timing  errors  s a very  slight  decrease  i n the  a ura, 
to  which  accelerometer-  errors  can  be  calibrated.  Gyro  calibration 
accuracy  is  virtually  unaffected. 
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7.2  TEST  HARDWARE  CONSIDERATIONS 

In  addition  to  the  general  hardware  capabilities  provided  by  the 
Singer  AN/ASN-90  inertial  measurement  set  and  associated  system  test 
rack  equipment,  two  additional  hardware  capabilities  are  required  to 
conduct  the  I MU  calibration  tests.  These  hardware  cons iderations  involve 
accurately  aligning  the  IMU  case  to  the  navigation  frame  and  precisely 
measuring  gimbal  angle  movements. 

7.2.1  IMU  Case  Alignment 

The  methodology  developed  for  IMU  calibration  assumes  that  the 
attitude  of  the  platform  relative  to  the  navigation  frame  is  accurately 
known.  The  IMU  case  frame  was  defined  synonymous  with  the  gimbal  neutral 
frame  defined  by  coincidence  of  the  neutral  IMU  roll,  pitch,  and  azimuth 
gimbal  axes  with  the  x (north),  y (west),  and  z (up)  navigation  frame 
axes.  Rotating  the  navigation  frame  through  the  [n]  coordinate  trans- 
formation produces  the  IMU  case,  and  gimbal  neutral  frames  and  rotating 
the  gimbal  neutral  frame  through  the  gimbal  anyle  [y]  coordinates 
transformation  produces  the  platform  frame. 

If  given  that  the  IMU  case  is  aligned  to  the  navigation  frame,  the 
twelve  relative  platform-to-navigation  frame  angular  orientations  ; 

required  for  separation  of  the  IMU  error  sources  shown  in  Figure  7-2  can 
be  obtained  by  gimbal  angle  rotations,  thus  obviating  any  need  to  rotate 
the  IMU  case.  There  are,  however,  several  advantages  to  mounting  the 
IMU  on  a rotating  tilting  precision  table.  (1)  It  assists  in  aligning 
the  case  to  the  navigation  frame.  (2)  It  allows  for  an  independent 
calibration  of  the  gimbal  synchro  angle  measuring  devices  foi  comparison 
with  the  Kalman  Filter  estimates.  And  (3)  additional  special  testing 
including  heading  sensitivity  can  readily  be  conducted.  For  these 
reasons,  an  IMTCO  notary  tilting  table  was  used  throughout  all  IMU 
calibration  testing. 
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Aligning  the  I MU  case  to  the  navigation  frame  is  a straightforward 
procedure.  The  I MU  case  is  mounted  on  the  tilt  table  with  use  of  the 
Singer  mounting  bracket.  The  necessary  three  angular  degrees  of 
freedom  a^e  provided  by  this  arrangement,  about  the  gimbal  neutral  IMU 
pitch  and  azimuth  axes  by  the  tilt  table  and  about  the  gimbal  neutral 
roll  axis  by  the  mounting  bracket.  The  IMU  is  brought  up  and  Singer's 
OCP  software  routines  are  executed  through  ground  alignment.  This  aligns 
the  x,  y,  and  z platform  frame  with  the  north,  west,  and  up  navigation 
frame.  The  gimbal  neutral  and  hence  IMU  case  frame  is  then  aligned  to 
the  navigation  frame  by  nulling  the  synchro  angle  readouts  using  tilt 
table  and  mounting  bracket  angular  rotations.  The  IMU  case  alignment 
accuracy  is  essentially  the  accuracy  of  the  synchro  measuring  devices 
or  6 arc  minutes.  Assuming  that  tilt  table  and  mounting  brackets  are 
unperturbed  or  can  be  accurately  returned  to  their  original  position, 
the  IMU  can  be  installed  for  additional  testing  without  repeating  the 
alignment  process. 

7.2.2  Gimbal  Angle  Measurements 

As  discussed  in  Section  4.2.1,  a gimbal  angle  measurement  or  its 
equivalent  is  required  for  separation  of  the  gyro  spin  axes  mass 
unbalance  error  sources.  Although  the  gimbal  synchro  angle  measuring 
devices  are  capable  of  measuring  small  differential  gimbal  angular 
positions  to  5 arc  seconds,  no  hardware  capability  existed  in  the 
Singer  AN/ASN-90  IMS  or  system  test  rack  equipment  to  obtain  measurements 
to  this  precision.  Initially,  low  precision  manual  gimbal  angle 
measurements  were  taken  with  the  visual  readout  indicators  in  the 
Singer  test  rack  equipment  to  support  the  IMU  calibration  development. 

Two  high  precision  synchro  measurement  technique  developments  were 
undertaken  that  would  hove  afforded  the  capability  of  automatic  data 
recording.  These  developments  were  never  completed.  A high  precision, 
manual  synchro  measurement  technique  was  eventually  perfected  for  use  in 
the  final  calibration  tests.  These  topics  are  discussed  in  the  jtquel . 
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The  Singer  system  test  rack  equipment  contains  roll,  pitch,  and 
azimuth  gimbai  synchro  angle  readouts  having  1 arc  minute  graduations. 
These  meters  can  be  read  to  roughly  0.2  arc  minute  precision.  Repeated 
I MU  testing  showed  that  the  azimuth  gimbai  angle  drifted  about  0.3  arc 
minutes  in  ten  minutes  when  torqued  at  earth  rate  in  test  attitude  No.  1. 
This  data  implies  a vertical  or  z-axis  platform  drift  rate  angular 
velocity  of  0.03  degrees  per  hour. 

Knowledge  of  the  actual  z-axis  platform  drift  rate  error  is 
sufficient  additional  information  to  enable  the  Kalman  Filter  to  separate 
the  gyro  mass  unbalance  error  sources.  To  utilize  this  data  in  the 
Kalman  Filter,  it  must  be  recognized  that  z-gyro  bias  and  soin  axis  mass 
unbalance  error  sources  completely  dominate  z-axis  platform  drift  rate 
error  when  the  axis  is  oriented  vertically.  This  fact  implies  that  these 
two  error  sources  are  statistically  correlated  such  that  the  sum  of  the 
platform  drift  rates  they  produce  has  an  initial  uncertainty  equal  to 
that  determined  by  the  above  measurements.  This  information  is  utilized 
by  the  Kalman  Filter  in  the  form  of  an  initial  covariance  between  the 
error  sources  having  a magnitude  given  by: 

P ( *Bz  • MSz  ) = [P ("Bz^Bz  ) + «* ' p(mSi . MSl ) - (O.Oif ]/2, 
p(  »Sz-*Bz)'  p(wBz.MSx) 

Entering  these  off-diagonal  covariances  into  the  Kalman  Filter's 
initial  covariance  matrix  permits  the  filter  to  estimate  each  of  the 
three  gyro  spin  axes  mass  unbalance  errors  to  an  accuracy  determined  by 
the  z-axis  platform  drift  rate  magnitude  or  roughly  0.03  degrees  per 
hour  per  g.  This  technique  of  separating  the  gyro  spin  axes  mass 
unbalance  errors  worked  quite  well  resulting  in  a reduction  of  30  of 
the  initial  uncertainty  of  1.0  degrees  per  hour  per  g.  Since  this 
technique  does,  however,  have  the  disadvantages  of  the  estimation 
accuracy  being  limited  by  z-axis  platform  drift  rate  and  the  measurements 
not  being  automatically  recorded  via  the  real-time  software,  additional 
gimbai  angle  measurement  techniques  were  investigated. 
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The  first,  attempt,  at  achieving  high  precision  automatic  gimbal 
angle  measurements  involved  a modification  of  the  Si nqer-Kearfott 
synchro-to-di gi tal  converter  obtained  from  the  TRW  INI  program.  This 
converter  contained  three  synchro  channels  with  a 12  bit  plus  sign  con- 


version capability  resulting  in  a measurement  precision  of  50  arc 
seconds.  The  philosophy  for  upgrading  this  converter  to  a measurement 
precision  of  approximately  one  arc  second  was  to  increase  the  precision 
of  the  device  at.  the  expense  of  dynamic  range.  This  was  achieved  by 
inserting  a variable  voltage  source  in  series  with  the  synchro  demodulator's 
output  to  define  a nominal  operating  point  of  zero  volts.  Amplification 
of  this  voltage  by  a factor  of  50  immediately  prior  to  anal og-to-di gi ta 1 
conversion  increased  the  precision  to  the  desired  one  arc  second.  The 
dynamic  range  of  the  gimbal  angle  measurements  about  the  nominal 
operating  point  was,  however,  reduced  from  180  degrees  to  69  arc  minutes. 
This  dynamic  range  is  about  an  order  of  magnitude  larger  than  required 
for  the  calibration  tests. 

The  Singer  synchro-to-digi tal  converter  was  modified  to  include  the 
required  electronic  circuitry.  Laboratory  tests  were  conducted  with  a 
synchro  standard  which  verified  the  one  arc  second  measurement  precision. 
Integration  testing  with  the  IMU  was  being  conducted  to  determine  the 
magnitude  of  potential  noise  problems  when  the  converter  failed.  The 
Singer  converter  was  a specialty  item  employing  non-standard  parts  which 
required  redesigning  a failed  printed  circuit  board  and  fabricating  a 
new  one  to  keep  the  device  in  service.  After  a second  failure,  a 
management  decision  was  made  that  the  converter  could  not  be  practically 
maintained  and  the  measurement  scheme  abandoned. 

The  second  attempt  at  achieving  high  precision  automatic  gimbal  angle 
measurements  involved  the  use  of  a production  Burr-Brown  anal og-to- 
frequency  electronic  converter.  The  operation  of  this  circuitry  is 
analogous  to  those  used  to  quantize  analog  accei erorneter  outputs  into 
velocity  pulses.  A second  synchro  was  used  in  conjunction  with  the  IMU's 
synchro  to  complete  the  control  transmitter  and  transformer  pair. 
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Unfortunately,  this  measurement  scheme  was  also  abandoned  when  it  became 
obvious  that  the  I/O  interface  electronics  and  real-time  data  collection 
software  could  not  be  completed  in  accordance  with  the  management  schedule 
for  completing  the  IMU  calibration  development  effort. 

The  final  attempt  at  obtaining  high  precision  gimbal  angle 
measurements  was  successful  and  is  depicted  in  Figure  7-1.  The 
measurements  were  manually  recorded  with  synchronization  between  the 
measurements  and  the  velocity  pulses  being  accomplished  with  a buzzer 
tied  into  the  teletype.  The  measurement  scheme  consisted  of  a serial 
arrangement  of  the  IMU  synchro  control  transmitter,  scott  tee  transformer, 
electronic  demodulator,  and  a Leeds  and  Northrop  bridge, 

A parametric  covariance  analysis  study  conducted  with  use  of  the  IMU 
calibration  Kalman  Filter  showed  that  only  the  azimuth  gimbal  angle 
measurement  was  required  during  the  first  platform  test  attitude  for 
accurate  estimation  of  the  gyro  spin  axes  mass  unbalance  error  coefficients. 
Employment  of  the  roll  and  pitch  and  even  the  azimuth  gimbal  angle 
measurements  during  platform  test  attitudes  2 through  12  did  not  signifi- 
cantly increase  the  filter's  estimation  accuracy.  Actual  gimbal  angle 
measurements  taken  during  the  IMU  calibration  tests  were  consistent  with 
these  findings. 

Calibration  of  the  azimuth  gimbal  angle  measurement  equipment's 
scale  factor  and  verification  of  its  high  resolution  capability  was 
accomplished  with  the  IMU  calibration  Kalman  Filter.  Referring  to 
Figure  7-2,  the  platform  orientation  in  test  attitudes  No.  5 and  No.  12 
is  such  that  the  Kalman  Filter  is  capable  of  accurately  estimating  the 
differential  azimuth  gimbal  angle  about  the  same  operating  point  as  test 
attitude  No.  1 (zero  dearees)  by  virture  of  the  mcasu;  ements  obtained 
from  the  level  accelerometers.  By  normalizing  the  manually  recorded 
Leeds  and  Northrup  data  to  the  corresponding  Kalman  Filter  estimates 
for  the  differential  azimuth  gimbal  angles,  a scale  factor  calibration 
is  obtained  from  the  same  data  used  for  calibrating  the  IMU. 
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Calibration  data  for  the  azimuth  gimbal  angle  measurement  technique's 
scale  factor  covering  the  last  eight  I MU  calibration  tests  is  shown  in 
Figure  7-4.  The  tests  were  conducted  by  positioning  the  azimuth  gimbal 
synchro  readout  indicator  located  in  the  system  test  rack  slightly  below 
the  nominal  setting  used  in  test  attitude  No.  1 and  then  allowing  the 
platform  to  drift  through  this  nominal  setting.  This  technique  ensures 
that  the  scale  factor  derived  from  the  No.  5 or  No.  12  platform  test 
attitude  corresponds  to  the  nominal  operating  point  used  in  test 
attitude  No.  1.  Scale  factor  determination  accuracy  is  a function  of 
both  the  Kalman  Filter  differential  azimuth  gimbal  angle  estimation 
accuracy  and  the  angular  excursion  of  the  azimuth  gimbal  during  the  test 
sequence.  Nominal  estimation  accuracies  of  1.4  arc  seconds'  one  sigma 
and  azimuth  gimbal  angle  excursions  of  3 arc  minutes  results  in  a scale 
factor  calibration  error  of  roughly  one  percent. 

The  calibration  data  of  Figure  7-4  shows  that  the  scale  factor 
differs  somewhat  from  day-to-day  tests  and  also  slightly  during  the 
nominal  one  hour  time  period  between  the  No.  5 and  Mo.  12  platform  test 
attitude  calibrations.  The  causes  for  these  scale  factor  variations  are 
not  completely  understood.  Replacement  of  400  Hz  demodulator  line  voltage 
reference  by  a precision  Hewlett  Packard  oscillator,  component  electro- 
static shielding  and  circuit  filtering  were  found  to  be  helpful  in 
reducing  the  scale  factor  variations  from  those  initially  encountered. 
Scale  factor  variations  for  the  No.  5 and  No.  12  test  attitudes  are  seen 
to  vary  only  a few  percent  for  the  same  I MU  calibration  test  over  the 
last  six  tests  conducted.  The  average  of  the  No.  5 and  No.  12  test 
sequence  scale  factors  is  used  in  conjunction  with  the  No.  ,1  test 
sequence  azimuth  gimbal  angle  measurement  for  the  formal  Kalman  Filter 
measurement.  Assuming  that  this  average  scale  factor  differs  from  the 
true  one  for  test  attitude  No.  1 by  this  same  small  percentage,  the 
scale-factor-induced  measurement  error  is  on  the  order  of  two  arc  seconds 
and  well  within  the  5 arc  second  noise  assumed  for  this  measurement. 

The  precision  of  the  azimuth  gimbal  angle  measurement  scheme  around 
the  zero  degree  operating  point  is  shown  in  Figure  7-5.  This  data  was 
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Figure  7-4.  Azimuth  Synchro  Scale  Factor  Calibration  Results 
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plotted  by  normalizing  the  Leeds  and  Northrup  bridge  data  to  the  first  and 
last  Kalman  Filter  differential  azimuth  gimbal  angle  estimates.  The 
Kalman  Filter  estimates  at  the  interior  data  points  serve  as  a reference 
with  an  accuracy  of  1.4  arc  seconds  one  sigma.  This  data  shows  that  the 
precision  of  measurement  technique  is  consistent  with  the  assumed  5-arc- 
secorid-measurement  noise. 
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SECTION  VIII 

I MU  ERROR  MODEL  REFINEMENTS 

As  stated  in  Section  2.2,  a major  goal  of  this  study  is  to  devise 
and  demonstrate  an  I MU  calibration  methodology  that  affords  a significant 
capability  fo"  conducting  IMU  error  model  development  and  validation 
efforts.  The  approach  taken  to  achieve  this  goal  has  been  to  formulate 
a methodology  employing  an  optimal  estimation  technique  that  estimates 
all  IMU  error  sources  in  a simultaneous  rather  than  the  sequential  manner 
used  in  contemporary  IMU  calibration  methods.  It  is  anticipated  that 
simultaneous  estimation  of  all  IMU  errors  sources  will  permit  IMU  error 
model  inadequacies  to  be  readily  detected.  Correlation  analyses  conducted 
in  conjunction  with  IMU  error  model  residual  Kalman  Filter  states,  added 
to  the  assumed  IMU  error  model  to  account  for  potential  error  model 
omissions,  should  permit  determining  the  functional  dependency  of 
detected  unmodeled  IMU  errors.  By  adding  the  detected  and  modeled  error 
source  to  the  initial  IMU  error  model  and  iterating  the  process,  develop- 
ment of  a complete  and  functionally  correct  error  model  appears  feasible. 

Kalman  Filter  processing  of  the  system  level  IMU  measurement,  data 
generated  during  the  laboratory  calibration  tests  resulted  in  the  detection 
of  both  an  unmodeled  gyro  quadrature  gee  am.,  accelerometer  attitude 
dependent  scale  factor  error.  Neither  of  these  error  sources  had  been 
included  in  the  Reference  1 error  model  developed  by  the  IMU  manufacturer . 
The  detection  and  modeling  of  these  error  sources  was  accomplished  as 
indicated  above  and  is  the  subject  of  this  section. 

8.1  INTI AL  KALMAN  FILTER  ESTIMATION  RESULTS 

Examination  of  the  initial  Kalman  Filter  IMU  error  estimation  results 
reguVes  a judgment  in  regard  to  the  adequacy  of  the  IMU  error  model 
incorporated  by  the  Kalman  Filter.  Two  approaches  were  found  useful  in 
arriving  at  this  judgment.  The  first  approach  involved  processing  the 
IMU  measurements  with  the  Kalman  Filter,  embodying  only  those  IMU  error 
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sources  thought  to  exist  The  degree  of  inconsistency  between  the 
Kalman  Filter  estimates  and  known  hardware  manufacturing  tolerances  for 
the  1 MU  error  sources  comprised  the  criterion  for  the  judgment.  The 
second  approach  involves  iterative  use  of  the  assumed  IMU  error  model 
residual  states  to  determine  the  magnitude  of  the  discrepancies.  This 
approach  is  well  suited  to  both  detecting  and  isolating  unmodeled  errors 
and  is  the  topic  of  Section  8.2. 

Two  examples  of  the  inconsistencies  obtained  between  the  Kalman 
Filter  estimates  and  the  known  manufacturing  tolerances  for  the  KT-73  IMU 
error  sources  are  shown  in  Figure  8-1.  The  left  and  right-hand  portions 
of  this  figure  contrast  the  Kalman  Filter  estimates  against  the  fi  gyro 
frame  misalignment  angle  and  gyro  long  time  constant  correlated  random 
drift  rate  manufacturing  tolerances.  Manufacturing  tolerances  are  indi- 
cated by  the  dashed  horizontal  lines.  In  addition  to  those  shown,  many 
other  similar  inconsistencies  were  observed. 

KT-73  IMh  manufacturing  methods  are  such  that  the  alignment  errors 
for  the  six  gyro-to-platform  misalignment  angles  are  nominally  held  to 
one  arc  minute.  Figure  8-1  illustrates,  however,  that  as  soon  as  the  IMU 
has  been  tested  with  its  platform  oriented  at  enough  different  attitudes 
to  permit  the  Kalman  Filter  to  separate  the  misalignment  angles  from  the 
other  error  sources,  recovered  estimates  for  these  angular  errors  are  a 
factor  of  twenty  times  the  manufacturing  tolerances.  Similar  results 
were  obtained  for  tne  correlated  gyro  random  drift  rate  error  with  recovered 
values  exceeding  production  tolerances  by  an  oruer  of  magnitude. 

The  conclusion  to  be  drawn  as  a result  of  the  Kalman  Filter 
estimates  for  many  of  the  IMU  errors  exceeding  their  manufacturing 
tolerances  by  an  order  of  magnitude  or  more  is  that  the  IMU  contains 
significant  errors  other  than  those  modeled  in  the  filter.  This  judgment 
was  particularly  easy  to  make  since  many  different  IMU  errors  exhibiting 
different  functional  dependencies  were  estimated  as  being  abnormally 
large.  As  discussed  in  Section  6.5,  had  only  one  error  source  been 
recovered  as  abnormally  large,  there  would  be  a strong  likelihood  that  a 
hardware  anomaly  existed. 
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Having  "lade  f he  judgment  that  a s i gn  i t icant.  IMII  error  model 
dim  repair  ■/•;•.)  existed,  it  be;:oi'°>  ne<  ev.sarv  to  determine  both  the  magni- 
tude ot  the  onntteo  error-  source  and  its  functional  dependency.  Suc- 
cessfully accomplishing  this  tasl  will,  o'  course,  confirm  the  initial 
I M| I error  model  deficiency  iudgienf.  The  manner  of  accomplishing  this 
task  is  by  the  employment  of  IMII  error  model  residual  Kalman  filter  states 
in  con  innct ion  with  iterative  processing  of  the  same  recorded  I Ml  i measure- 
ment data. 


The  initial  obiecfive  ot  iterative  Kalman  filter  processing  of  the 
1 Ml i measurement  data  is  to  obtain  consistency  between  the  initial  uncer- 
tainties assumed  for  the  filter's  initial  covariance  matrix  and  the 
filter's  estimates  of  the  corresponding  error  source.  Consistency  implies 
that  the  Kalman  filter  estimates  for  the  random  bias  error  states  are 
within  three  standard  deviations  of  the  uncertainties  assumed  for  the 
initial  covariance  matrix.  Three  standard  deviations  is  a rule-of-thumb 
used  in  conjunction  with  the  Gaussian  distribution  and  includes  roughly 
99.5-  of  all  experiment  trials.  Consistency  is  sought  for  filter  estimates 
corresponding  to  each  platform  test  attitude  as  well  as  final  filter 
estimates  corresponding  to  the  final  test  sequence. 

Obtaining  consistency  between  the  Kalman  filter's  estimates  and  its 
initial  covariance  matrix  is  required  for  proper  or  accurate  estimation 
by  the  Kalman  filter  algorithm.  Without  consistency,  the  filter's  esti- 
mation accuracy  is  not  giver,  by  its  covariance  matrix  and  estimation  errors 
are  unknown.  Ry  adding  the  IMU  error  model  residual  states  to  the  Kalman 
filter  state  vector  and  determining,  by  iteration,  their  uncertainties 
for  insertion  in  the  initial  covariance  matrix,  consistency  can  be  obtained. 
Having  obtained  consistency,  Kalman  Filter  error  source  estimation  errors 
are  now  accurately  kr.cwn  and  given  by  the  filter’s  covariance  matrix. 
Consistency  permits  the  Kalman  filter  estimation  results  to  be  used  with  a 
high  degree  of  confidence  in  conjunction  with  correlation  studies  to  deter- 
mine the  magnitude  and  functional  dependency  of  unmooeled  IMU  errors. 
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An  illustration  of  obtaining  Kalman  filter  optimal  estimate  consis- 
tency through  the  use  of  error  model  residual  states  is  shown  in  Figure 

The  left  and  right-hand  portions  of  this  figure  contain  similar 
data  for  the  x and  y-gyros,  respectively.  Three  gyro  error  model  residual 
states  were  included  in  the  Kalman  Filter  state  vector,  a drift  rate  error 
for  each  gyro.  These  drift  rate  errors  were  modeled  as  random  biases  and 
were  reset  to  their  assumed  initial  uncertainty  at  the  start  of  each  plat- 
form test  sequence.  In  effect,  this  resulted  in  a 3 errors  per  test 
sequence  times  12  test  sequences  or  36  additional  error  sources  for  the 
Kalman  Filter  to  estimate.  These  additional  36  error  sources  do  not,  of 
course,  permit  the  Kalman  Filter  to  separate  and  accurately  estimate  the 
other  1 MIJ  error  sources  when  the  assumed  initial  uncertainties  of  the 
residual  states  are  large  compared  with  the  other  errors. 

The  error  model  residual  states  depicted  in  Figure  3-2  could  also 
be  termed  as  gyro-attitude-dependent  bias  errors,  since  they  are  modeled 
as  being  constant  for  a given  platform  attitude  relative  to  the  I MU  case 
but  assume  a new  value  anytime  the  platform's  orientation  is  changed. 

The  Kalman  Filter's  estimates  for  these  attitude-dependent  biases  is 
sh"wn  for  each  of  the  twelve  platform  test  orientations  The  upper  and 
lower  brackets  about  each  estimate  represent  the  standard  deviation  of 
the  filter's  estimation  error  as  obtained  from  the  covariance  matrix  at 
the  end  of  each  platform  test  sequence.  Since  filter  consistency  or 
proper  convergence  has  been  obtained,  these  estimation  errors  are  accurate. 

Filter  consistency  was  obtained  by  iteratively  increasing  the 
initial  uncertainty  of  the  attitude-dependent  bias  error  states  in 
conjunction  with  reprocessing  of  the  same  recorded  I MU  measurement  data 
until  the  filter's  estimates  for  these  errors  were  within  three  standard 
deviations  of  their  initial  uncertainty.  Initial  trial  values  were 
chosen  to  be  the  same  as  the  correlated  gyro  random  noise  of  0.005  degree 
per  hour  with  about  two  iterations  being  required  to  reach  the  consistency 
value  of  0.2  degree  per  hour  shown  in  Figure  8-2.  New  trial  values  were 
taken  to  be  the  maximum  of  the  Kalman  Filter's  estimates  for  the  attitude- 
dependent  biases  on  the  previous  trial.  Figure  8-2  shows  that  maxmium 
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Kalrnan  Filter  estimates  for  the  x and  y-gyro  attitude-dependent  biases 
are  roughly  0.4  deqree  per  hour,  which  is  within  two  standard  deviations 
of  the  initial  uncertainty  of  0.2  degree  per  hour.  Two  standard  deviations 
include  about  95%  of  the  trial  outcomes  for  Gaussian  statistics. 

Examination  of  the  Kalman  Filter's  estimation  errors  for  the  other 
I MU  error  sources  at  the  end  of  the  twelfth  platform  test  sequence  indi- 
cates that  the  filter  is  unable  to  accurately  calibrate  all  of  them  due 
to  the  presence  of  the  atti tude-deoendent  biases  with  large  initial 
uncertainties.  In  effect,  a solution  is  precluded  because,  although  each 
platform  test  sequence  restilts  in  three  add i ti onal  equations  for  seDaration 
of  the  I MU  calibration  parameters,  three  additional  unknowns  are  added 
also.  The  filter's  estimates  for  all  I MU  error  sources  are,  however, 
within  three  standard  deviations  of  their  initial  uncertainties.  This 
consistency  is  obtainable  only  through  the  use  of  error  model  residual 
states . 

In  addition  to  a desire  to  upgrade  the  KT-73  IMU  error  model,  the 
exact  functional  dependency  of  the  errors  producing  the  0.4  deqree  per 
hour  gyro  drift  rate  for  some  platform  attitudes  must  be  determined  so 
that  it  can  be  explicitly  included  in  the  Kalman  Filter  system  dynamics 
matrix  to  enable  the  filter  to  accurately  calibrate  the  other  IMU  errors. 
Determining  the  exact  functional  dependency  of  error  model  residual  states 
amounts  to  performing  a trial  and  error  analysis  to  determine  the  para- 
meter, if  any,  which  varies  in  magnitude  or  vector  direction  with  changirg 
platform  attitude  with  which  the  gyro  drift  rate  error  is  correlated. 

The  results  of  this  analysis  are  shown  in  Figure  R-2  and  indicate  a high 
degree  of  correlation  Detween  gyro  drift  rate  magnitude  and  the  orientation 
of  the  gyro  output  axis  relative  to  the  specific  force  vector.  The  speci- 
fic force  vector  is  essentially  the  gravity  vector  for  the  fixed  IMU  base 
situation  and  has  been  assumed  to  be  oriented  parallel  with  the  y-axis 
of  the  graph  in  Figure  8-2.  Note  in  this  figure  that  the  correlation 
appears  to  be  especially  high  for  those  platform  test  sequences  where 
estimation  accuracy  is  high,  as  indicated  by  closely  spaced  brackets. 
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To  confirm  the  apparent  error  model  residual  gyro  output  axis 
specific  force  dependency,  three  additional  states  were  added  to  the  Kalman 
Filter  explicitly  containing  this  functional  dependency  and  an  initial 
uncertainty  of  0.4  degree  per  hour  per  g.  The  initial  uncertainty  of  the 
error  model  residual  states  was  reduced  to  the  gyro  random  noise  level  of 
0.005  degree  per  hour  and  the  recorded  IMU  measurement  data  reprocessed 
by  the  modified  Kalman  Filter.  Kalman  Filter  estimation  results  for  the 
error  model  residual  or  gyro  attitude-dependent  bias  states  are  shown  in 
Figure  8.3.  Brackets  about  the  filter's  estimates  indicate  the  initial 
uncertainty  of  0.005  degrees  per  hour  and  illustrates  that  no  estimates 
were  obtained  greater  than  this  value.  These  results  confirm  the 
existence  of  a gyro  output  axis  specific  force  dependent  error  source. 

The  results  of  Figure  8-3  along  with  the  now  obtained  accurate  calibration 
of  the  other  gyro  and  related  error  sources  showing  consistency  with  their 
initial  uncertainties  also  indicate  the  absences  of  any  additional  gyro- 
related  errors  of  magnitude  appreciably  larger  than  the  0.005  degree  per 
hour  error  model  residuals. 

The  physics  of  why  a specific  force  acting  along  the  output  axis 
of  the  KT-73  gyroflex  gyros  produced  a drift  rate  error  were  initially 
perplexing.  It  is  well  known  that  a mass  unbalance  in  conjunction  with  a 
specific  force  acting  along  a gyro  input  or  spin  axis  can  produce  a drift 
rate  error.  However,  the  vector  relationships  are  such  that  it  is  impos- 
sible for  this  to  occur  when  it  concerns  the  output  axis.  Mr.  Bud 
Bleichroth  of  AGMC  eventually  isolated  the  cause.  The  gyroflex  gyro  is 
a two-degree-of-freedom,  dry  flexure  joint  suspended,  free  rotor  gyro. 

Its  unique  feature  is  a special  rotating  flexure  suspension  on  one  end  of 
’ shaft,  supporting  and  pivoting  the  wheel.  The  other  shaft  end  is 
driven  by  a synchronous  hysteresis  motor.  If  the  flexure  suspension  joint 
is  not  machined  in  a completely  symmetrical  manner,  then  a specific  force 
acting  along  the  output  axis  acts  in  conjunction  with  the  flexure  asymmetry 
to  yield  a torque  which  produces  the  drift  rate  error.  This  error  source 
has  been  tenned  a gyro  quadrature  g error  by  the  Singer  engineers  and  has 
been  specified  at  1.0  degree  pe"  hour  per  g for  the  SRAM  KT-76  IMU. 
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8.2.2  Gyro  G-Squared  Error  Correlation  Studies 

Havinq  found  a rather  major  omission  in  the  Reference  1 I MU  error 
model  raises  the  question  as  to  the  possible  omission  of  additional  qyro 
error  sources  that  are  not  strictly  negligible  relative  to  the  inertial 
instrument  random  disturbances.  This  concern  seems  particularly  justified 
for  specific  force  related  gyro  errors  sinc.e  the  Reference  1 I MU  error 
model  was  formulated  for  navigational  error  analyses  where  accelerations 
are  generally  of  short  time  durations.  Since  the  I MU  calibration  tests  are 
conducted  with  the  platform  at  numerous  attitudes  relative  to  the  specific 
force  vector  for  extended  periods  of  time,  gyro  error  sources  insignificant 
for  navigational  accuracy  computations  may  not  be  negligible  in  consider- 
ation of  Kalman  Filter  IMU  calibration  accuracy  constraints.  For  these 
reasons,  it  was  decided  to  conduct  correlation  studies  on  the  next  most 
logical  gyro  specific  force  related  errors  omitted  from  the  Reference  1 
IMU  error  model  which  are  additional  g-squared  gyro  errors. 

The  Reference  1 gyro  error  model  contains  an  anisoelastic  drift  rate 
error  having  a magnitude  proportional  to  the  product  of  specific  forces 
acting  along  the  output  and  spin  axes  of  the  gyro.  Taking  possible  com- 
binations of  three  things,  input-spin-output  axes  soecific  forces,  taken 
two  at  a time,  leaves  five  additional  g-squared  gyro  error  sources  to  be 
evaluated. 

The  philosophy  for  evaluating  the  magnitude  of  the  gyro  g-squared 
error  sources  can  be  understood  with  the  aid  of  Figure  8-4.  The  recorded 
IMU  measurement  data  is  processed  by  two  versions  of  the  IMU  calibration 
Kalman  Filter.  The  first  is  the  nominal  or  baseline  version  while  the 
second  version  has  been  modified  to  include  three  additional  g-squared 
gyro  drift  rate  error  states.  There  is  one  error  state  per  gyro  with 
each  g-squared  error  having  the  same  gyro  axes  specific  force  dependency. 
The  initial  uncertainty  of  these  error  sources  is  found  by  iteration  such 
that  the  final  Kalman  Filter  g-squared  error  estimate  is  within  one  to 
two  standard  deviations  of  the  initial  uncertainty.  Generally,  only  one 
iteration  was  required  and  all  initial  uncertainties  were  found  to  be 
between  0.003  and  0.005  degrees  per  hour  per  g-squared. 
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The  evaluation  criteria  were  both  the  Kalman  Filter's  estimate  of 
the  g-squared  error  source  and  the  degree  to  which  the  addition  of  this 
error  source  to  the  I MU  error  model  reduced  the  filter's  estimate  for 
the  error  model  residual  or  gyro  attitude-dependent  bias  states.  The 
relative  reduction  in  gyro  attitude  dependent  bias  filter  estimates  is 
shown  in  Figure  8-4  for  the  addition  of  gyro  output  axis  g-squared  error 
states  to  the  Kalman  Filter  state  vector.  The  Kalman  Filter's  estimate 
for  these  g-squared  errors  were  0.004  and  0.003  degree  per  hour  per  g-squared 
for  the  x and  z-gyros,  respectively.  Initial  uncertainties  were  0.003 
and  0.004  degrees  per  hour  per  g-squared.  Substantial  reductions  in 
the  z-gyro  attitude-dependent  biases  are  seen  to  occur  while  little 
improvement  is  obtained  by  adding  an  output  axis  g-squared  error  source 
to  the  x-gyro  error  model . 

Of  the  five  sets  of  three  g-squared  error  sources  evaluated,  those 
showing  an  output  axis  g-squared  functional  dependency  rated  highest  in 
terms  of  the  above  criteria.  Since:  (1)  the  filter's  estimates  for  all 
the  g-squared  error  sources  were  less  than  the  correlated  gyro  random 
drift  rate  magnitudes,  (2)  only  a marginal  reduction  in  the  error  model 
residual  states  was  obtained,  and  (3)  inclusion  of  the  g-squared  errors 
in  the  Kalman  Filter  state  vector  did  not  significantly  affect  the  filter's 
estimates  for  the  I MU  error  sources,  it  was  concluded  that  the  g-squared 
errors  could  be  neglected. 

8.2.3  Heading  Sensitivity  Tests 

Prior  to  actually  trying  the  IMU  calibration  technique  in  the 
laboratory  and  demonstrating  its  feasibility,  there  was  a high  degree  of 
concern  that  potential  inertial  instrument  heating  and  magnetic  errors 
induced  by  gyro  slewing,  and  variations  in  cooling  air  and  magnetic  fields 
with  relative  platform- to- IMU  case  orientations,  might  mask  or  induce 
significant  errors  in  the  determination  of  the  IMU  calibration  parameters. 
Several  special  tests  were  conducted  to  confirm  or  negate  this  hypothesis. 

The  first  series  of  tests  were  modifications  of  the  basic  IMU  calibration 
test  and  the  second  series  constituted  heading  sensitivity  tests. 
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In  the  first  series,  d test  was  conducted  to  determine  if  rapidly 
torquing  the  platform  from  one  test  sequence  to  another  would  result  in 
significant  gyro  heating  and  an  associated  thermal  drift  rate  error. 
Usually,  a calibration  test  was  conducted  by  slewing  the  platform  to  new 
test  attitudes,  using  the  analog  torquing  mode  with  a 30  (x,  y-gyro) 
and  90  (z-gyro)  degree  per  minute  slew  rate,  and  immediately  recording 
the  I MU  measurement  data.  A special  test  was  conducted  by  allowing  the 
gyros  a fifteen  minute  cooling  down  period  between  their  time  of  arrival 
at  a new  test  attitude  and  the  time  data  recording  commenced.  Two  tests 
were  conducted,  one  using  five  and  the  other  using  ten  minute  measurement 
time  intervals  per  test  sequence.  Kalman  Filter  I MU  error  source  esti- 
mates obtained  from  these  tests  agreed  with  those  obtained  for  the  nominal 
tests  to  within  the  accuracy  predicted  by  the  filter's  covariance  matrix. 
These  results  indicate  that  gyro  heating  due  to  rapid  torquing  is  not  a 
significant  source  of  error  for  the  KT-73  I MU  gyros. 

In  another  special  test,  some  of  the  platform  test  attitudes  were 
changed  to  alternate  ones  to  determine  if  the  different  magnetic  fields 
encountered  by  the  inertial  instruments  were  sufficient  to  cause  a shift 
in  the  recovered  I MU  calibration  parameters.  The  results  of  this  test 
were  negative. 

Another  special  test  was  conducted  to  demonstrate  that  the  gyro 
drift  rate  errors  could  be  software  compensated.  Test  results  showed  a 
reduction  in  platform  total  drift  rate  from  about  1.0  to  0.01  degree  per 
hour  for  many  of  thp  test  attitudes.  Gyro  spin  and  input  axes  mass 
unbalan  e,  quadrature  g,  scale  factor,  misalignment  angles,  anisoelastic 
and  bias  error  sources  were  compensated  for  by  modification  of  the  gyro 
earth  rate  torquing  magnitudes.  I MU  calibration  parameter  values  recovered 
from  this  test  were  also  consistent  with  those  obtained  from  the  baseline 
calibration  tests. 

The  second  series  of  special  tests  were  conducted  to  determine 
heading  sensitivities.  The  term  heading  sensitivity  is  used  to  denote 
gyro  drift  rate  errors  produced  by  cooling  air,  temperature,  and  magnetic 
field  variations  associated  with  changes  in  platform-to-IMU  case  azimuth 
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orientations  for  the  situation  of  both  platform  and  I MU  case  being 
approximately  level.  To  conduct  the  heading  sensitivity  test  in  the 
laboratory,  the  IMU  case  was  maintained  in  alignment  with  "he  navigation 
frame  while  the  platform  was  indexed  one  complete  revolution  in  azimuth 
in  30-degree  increments.  The  platform  was  stabilized  relative  to  the 
navigation  frame  by  application  of  earth  rate  torquing.  IMU  measurements 
were  recorded  at  30-second  intervals  for  a duration  of  ten  minutes  at 
each  test  attitude. 

Kalman  Filter  heading  sensitivity  estimation  results  are  shown  in 
Figure  8-5.  Four  drift  rate  errors  are  depicted  to  illustrate  the  filter' 
ability  to  separate  the  various  error  sources.  Gyro  attitude-dependent 
drift  rate  error  was  modeled  as  a random  bias  and  reset  to  its  initial 
uncertainty  of  0.005  degree  per  hour  at  the  start,  of  each  new  platform 
test  sequence.  These  three  Kalman  Filter  error  states  are  the  same  as 
those  discussed  in  previous  sections  and  reflect  the  gyro  heading  sensi- 
tive drift  rate  errors  at  each  platform  orientation.  Gyro  total  drift 
rate  error  is  a composite  error  composed  of  all  gyro  errors  considered 
internal  to  the  instrument.  Platform  gyro-dependent  drift  rate  error  is  a 
composite  error  composed  of  gyro  total  drift  rate  plus  gyro  scale  factor 
and  gyro-to-pl atform  misalignment  angle  errors.  Platform  total  drift  rate 
is  a composite  error  composed  of  platform  gyro-dependent  drift  rate  plus 
IMU  case-to-navigation  frame  misalignment  angle,  g f mba 1 synchro  bias,  and 
differential  angle  errors. 

The  estimates  for  the  composite  errors  were  computed  by  summing  the 
appropriate  Kalman  Filter  estimates  for  the  individual  component  error 
states.  Estimation  errors  for  these  composite  errors  were  computed  by 
using  the  expectation  operator  and  appropriate  on  and  off  diagonal 
covariance  matrix  elements.  Estimation  error  standard  deviations  for 
all  x-gyro  and  x-axis  platform  error  quantities  plotted  in  Figure  8-5 
fall  between  0.003  and  0.006  degrees  per  hour. 

Platform  total  and  gyro-dependent  drift  rate  error  estimates  show 
similar  characteristics.  The  seemingly  abnormal  -0.01  degree  per  hour 
platform  total  drift  rate  error  at  a 90  degree  platform  heading  w;s 
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Figure  8-5.  Heading  Sensitive  Drift  Calibration  Results 
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produced  by  an  unusually  large  error  in  positioning  the  platform  via  the 
synchro  angle  readouts  at  this  particular  attitude.  Both  composite  plat- 
form drift  rate  error  character! sties  tend  to  indicate  a heading  sensi- 
tivity; however,  comparison  of  these  characteri sti cs  to  that  for  total 
gyro  drift  rate  error  and  examination  of  other  Kalman  Filter  estimation 
results  shows  the  apparent  heading  sensitivity  to  be  due  to  an  x-gyro 
scale  factor  error.  Except  for  one  point,  all  x-gyro  attitude-dependent 
bias  estimates  are  less  than  the  KT-73  I MU  heading  sensitive  drift  rate 
error  specification  of  0.004  degree  per  hour. 

The  conclusion  to  be  drawn  from  these  special  and  heading  sensitive 
test  results  is  that  thermal  and  magnetic  field  variations  over  the  dif- 
ferent platform-to-IMU  case  test  orientations  used  are  on  the  order  of 
the  gyro  random  disturbances  and  do  not  significantly  affect  gyro  error 
source  calibration  accuracy. 

8.3  ACCELEROMETER  MODELING  REFINEMENTS 

Initial  Kalman  Filter  estimation  results  for  the  accelerometer  error 
sources  were  reasonably  consistent  with  those  expected.  Three  error  model 
residual  states,  modeled  as  accelerometer  attitude-dependent  scale  factor 
errors,  had  been  included  in  the  accelerometer  error  model  to  account  or- 
al low  for  possible  error  model  discrepancies.  Functional  model inq  of  the 
residual  states  is  not  critical  and  attitude-dependent  bias  errors  would 
have  served  equally  well.  With  the  initial  variances  of  these  error  model 
residual  states  set  at  relative1'/  small  values,  the  only  detected  incon- 
sistency of  the  estimation  results  was  a non-repeatibi I ity  of  the  accel- 
erometer bias  and  scale  factor  error  estimates  to  within  the  accuracy 
predicted  by  the  filter's  covariance  matrix.  This  problem  could  be 
ameliorated,  however,  by  sufficiently  increasing  the  initial  variance  of 
the  error  model  residual  states.  Increasing  the  initial  uncertainty  of 
the  error  model  residual  states  results  in  corresponding  increases  in 
estimation  error,  as  predicted  by  the  filter's  covariance  matrix,  but 
does  not  substantially  affect  the  error  state  estimates. 
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The  above  procedure  of  adjusting  the  error  model  residual  state 
initial  uncertainties  to  obtain  error1  state  estimate  test-to-test  repeat.- 
ibili.y  to  within  the  filter's  predicted  calibration  error  was  considered 
unsatisfactory  because  it  did  not  answer  the  basic  questions.  Are  there 
thermal  and/or  magnetic  field  variations  present  at  the  various  platform 
test  attitudes  that  perturb  the  nominal  accel erometer  scale  factor  and 
bias  errors?  How  accurately  should  test-to-test  calibration  results 
repeat,  considering  normal  IMU  hardware  turn-on  variations?  In  an  attempt 
to  answer  these  questions,  it  was  decided  to  examine  in  detail  the 
repeat! bi 1 i ty  of  the  accelerometer  misalignment  angle  errors,  since  these 
error  sources  arc  essentially  static  and  not  subject  to  hardware  turn-on 
variati ons . 

8.3.1  Accelerometer  Non-Orthogonality  Angles 

An  examination  of  the  filter's  covariance  matrix  elements  corres- 
ponding to  the  platform-to-accelerometer  frame  misalignment  angle  errors 
indicates  that  the  filter  is  capable  of  reducing  the  initial  manufacturing 
tolerance  uncertainty  of  1.0  arc  minutes  to  only  0.7  arc  minutes.  In 
contrast  to  the  accelerometer  bias  and  scale  factor  error  estimates, 
the  platform-to-accelerometer  misalignment  angle  error  estimates  did 
repeat  to  within  the  accuracy  predicted  by  the  filter's  covariance  matrix. 
As  such,  the  error  state  estimates  and  estimation  errors  for  the  acceler- 
ometer misalignment  angles  as  defined  in  Section  3.3.2  are  unsuitable  for 
a detailed  study  of  'she  accelerometer  error  model  as  en1  isioned  in  the 
previous  section.  A reformulation  of  t he  platfo>'m-t.o-accelei  ometer 
misalignment  angles,  that  results  in  estimation  accuracies  limited  largely 
bv  the  accelerometer  random  disturbances,  is  necessary  to  conduct  this 
study . 


The  limitation  encountered  by  the  Kalman  Filter  in  attempting  to 
accurately  estimate  tne  angular  misalignments  of  the  accelerometers  from 
the  platform  frame  by  use  of  only  the  accelerometers  measurements  is  that 
these  measurements  indicate  a misalignment  of  the  accel erometer ' s electrical 
null  axis  from  an  assumed  orientation  relative  io  the  gravity  vector. 

The  misalignment  of  the  accelerometer's  electrical  null  axis  from  the 
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the  filter  ft  on  .uisrately  separating  tin.1  a.  < e I erorne  tfci  - to -p  I a t f orm  frame 
nn  sa  I lqn.uent  wunr-  errot  s from  the  compos  i te  accelerometer  electrical 
■ ill  axis  (.ilt  i"  . ors  Separation  of  this  error  from  the  composite  tilt, 

et rot  rejoices  m a. curate  measure  of  the  plat  form's  annular  attitude 

which  is  not  available.  Note  that  the  K 1 - / 3 I MU  synchros  are  precise  to 

‘>  arc  seconds  but.  accurate  to  only  b art.  minutes  and,  therefore,  too 

inaccurate  tor  this  measurement 

I rom  the  above  discussion,  if  van  be  rone  I tided  that  accurate 
estimation  of  some  type  ot  at  c e 1 to  omet.er  in  1 so  1 1 grimeri  l angles  can  only  be 
achieved  by  el  i mi  mil ing  to  first  order  the  g bubal  synchro  angle  bias 
errors  from  the  detei  mi  oat  ion  process.  Itiis  will  require,  however,  that, 
a measurement  re te rerun  or  its  equivalent  other  than  the  qravity  vector 
be  used.  the  most  logical  choice  tor  an  accelerometer  measurement 
reference  other  h n the  qravity  vector  is  another  accelerometer  measure- 
ment and  this  im.'it  that  .u  i eicrometer  m i sa I i qnments  from  an  orthogonal 
set  is  the  quant  1 1 ■ to  be  estimated.  To  illustrate  the  concept,  consider 
fh».'  No.  4 platter:  test,  sequeie  o shown  in  figure  / -? . Assume  that  the 
y and  / ■ acr eicrometer  time  varying  outputs  are  given  to  first  order  by: 

MO  c '®KAy  co*(d9*(0)  + ABy  4 %(*) 

MO  * co*  [ **  (^9x(»)  + <*„)]  + A0I  + ARl(f) 

where 

i is  the  non  orthogena 1 tty  error  angle  of  tne  /-accelerometer 

l X 

input  dAis  with  respect  to  the  /-accelerometer  input  axis. 

i is  produced  by  rotating  the  /-accelerometer  about  the 
Z X 

x-axis. 
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These  expressions  can  be  solved  for  a by  using  a trigonometric 
formula  and  a small  angle  approximation  to  yield: 


This  expression  demonstrates  that  accelerometer  non-orthogonality 
error  angles  can  be  determined  by  simul taneously  taking  pairs  of  acceler- 
ometer measurements  and  to  first  order,  the  solution  is  independent  of 
the  gimbal  synchro  angle  bias  errors,  which  is  the  desired  result. 
Estimation  of  the  non-orthogona 1 i ty  error  angles  by  the  Kalman  filter 
requires  that  these  errors  be  expressed  in  terms  of  the  platform-to- 
accelerorneter  misalignment  angle  error  states.  These  expressions  can  be 
derived  by  use  of  coordinate  transformations  that  "square  up"  the  acceler- 
ometer frame  and  then  rotate  it  into  the  platform  frame. 

An  orthogonal  acce I e route ter  frame  can  be  defined  in  the  following 
manner.  Rotate  the  y-acce l erometer , in  the  plane  defined  by  the  x and 
y-accelerometer  input  axes,  by  the  small  angle  a z into  orthogonality  with 
the  x-accelerometer . Rotate  the  z-a^cel erometer  by  the  small  angle  a 
about  the  x-accelerometer  axis  into  orthogonality  with  the  y-accelerometer. 
And  finally,  rotate  the  z-acceleroneter  by  the  small  angle  cx  about  the 
y-accelerometer  axis  into  orthogonal i ty  with  the  x-accelerometer.  All 
rotations  are  done  in  accordance  with  the  right-hand  rule.  The  coordinate 
transformation  rotating  the  orthogonal  accelerometer  frame  into  the  frame 

defined  by  the  ohysical  orientation  of  the  accelerometer  input  axes  is 

. _A 
given  by 

OA 
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To  relate  these  non-orthogona  I i tv  error  angles 
accelerometer  misalignment  angles  defined  in  Section 
definition  of  a olatform-to-orthogonal  accelerometer 
small  anole  coordinate  transformation  matrix  Cp 


to  the  platform- to- 
3.3.2  requires 
frame  orthogonal 
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Equating  the  platform-to-accelerometer  frame  coordinate  trans- 
formations yields: 
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The  matrix  equation  can  be  solved  for  the  expressions  relating  the 
three  accelerometer-  non-orthogonality  to  the  six  platform-to-accelerometer 
misalignment  angle  errors. 
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8.3.2  Accelerometer  Attitude  Dependent  Errors 

Kalman  Filter  estimation  results  for  the  a accelerometer 
non-orthogonal i ty  error  angle  are  shown  in  Figure  8-6.  Calibration 
results  are  depicted  in  three  sets  of  five  IMU  tests  and  illustrate  the 
history  of  improvements  made  in  estimating  this  error  source  by  refining 
the  IMU  error  model.  For  each  series  of  five  IMU  tests,  the  brackets  or 
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boundaries  superimposed  about  either  side  of  the  calibrated  average 
estimate  (average  estimate  not  shown)  represent  a standard  deviation  of 
the  estimation  error.  Total  spacing  between  the  boundaries  is,  therefore, 
two  standard  deviations.  This  data  was  computed  using  Kalman  Filter 
estimation  results. 

Ootiinal  Estimate: 

AAA 

- yX2 

Estimation  Error: 

P(“XZ , <■„ ) = E[(y>i  - /„)(/,,  - r„)] 

* p(r,z'  r,z)  + p(r,z • r„)-2P(r>l.  r„) 

The  estimation  error  shown  for  the  first  five  calibration  tests  was 
computed  using  the  accelerometer  error  modeling  given  in  Reference  1. 

This  error  model inq  resulted  in  an  estimation  error  standard  deviation  of 
0.7  arc  seconds  and  was  essentially  determined  by  the  accel erometer  random 
disturbances.  Gaussian  statistics  imply  that  68%  of  the  estimates  should 
fall  between  the  estimation  error  bounds  on  an  ensemble  average  basis  with 
virtually  all  estimates  falling  within  three  standard  deviations  of  the 
five  test  estimation  average.  Figure  8-6  clearly  shows  that  the  estimation 
results  obtained  for  the  first  five  tests  are  grossly  inconsistent  with 
these  statistical  averages  and  suggests  the  presence  of  unmodeled  error 
sources . 

As  discussed  in  Section  7.1.4,  the  real-time  data  recording  software 
suffered  from  a measurement  timing  instability  of  15.0  msec.  This 
instability  was  eventually  reduced  to  around  1.0  msec  by  the  fifteenth 
calibration  test.  It  was  necessary  to  modify  the  Kalman  F i 1 tor  by 
adding  one  error  state  to  the  state  vector  and  three  disturbances  to  the 
velocity  error  states  to  account  for  this  error  since  it  was  somewhat 
larger  than  the  permissible  0.1  msec  timing  error.  The  longer  term 
correlated  timing  error  was  mode  ed  as  a pseudo  accelerometer  scale 
factor  error,  one  error  affecting  all  three  accelerometer  outputs,  having 
an  initial  uncertainty  of  one  ppm  and  a 250  second  correlation  time. 
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The  shorter  term  disturbances  on  the  velocity  error  states  was  modeled 
as  being  proportional  to  the  amount  of  gravity  measured  by  an  acceler- 
ometer. Disturbance  modeling  is  such  that  it  produced  a velocity 
uncertainty  correspond i ng  to  one  CAPRI  pulse  in  a measurement  interval 
when  an  accelerometer  experiences  full  gravity  along  its  input  axis. 

The  estimation  results  shown  in  Figure  8-6  for  the  second  series 
of  I MU  calibration  tests  illustrate,  tests  10-14  compared  with  1-5,  how  a 
reduction  in  timing  instability  from  15.0  to  1.0  msec  and  a modification 
of  the  Kalman  Filter  to  account  for  the  residual  timing  errors  improves 
the  non-orthogonality  error  angle  estimates.  Figure  8-7  shows  the  filter's 
estimates  for  the  clock  error  state  at  the  end  of  each  platform  test 
sequence  for  I MU  calibration  Test  No.  24.  This  error  state  is  reset  to 
its  initial  uncertainty  at  the  beginning  of  each  test  seouence. 

Figure  8-6  also  shows,  however,  that  the  spread  in  accelerometer 
non-orthogonality  error  angle  estimates  is  still  inconsistent  with 
Gaussian  distribution  ensemble  averages  for  I MU  calibration  Tests  10-14. 
This  inconsistency  again  indicates  the  presence  of  unmodeled  accelerometer 
errors  and  a special  test  was  designed  to  determine  the  presence  cf 
accelerometer  attitude  dependent  errors  due  to  thermal  and  magnetic  field 
variations  with  changes  in  platforn-to-IMU  case  orientations.  Test 
philosophy  was  to  monitor  the  number  of  velocity  pulses  generated  by  an 
accelerometer  when  its  input  axis  was  parallel  to  the  gravity  vector  and 
its  case  axis  was  at  different  azimuth  attitudes  relative  to  the  IMU  case. 
Accelerometer  error  modeling  of  Reference  1 and  the  alignment  accuracy  of 
the  accelerometer  input  axis  to  the  gravity  vector  by  use  of  the  gimbal 
synchro  angle  devices  predict  no  more  then  one  velocity  pulse  variation 
over  a 30  second  measurement  interval,  lo  eliminate  the  possibility  of 
the  real-thne  software  timing  errors  affecting  the  data,  velocity  pulses 
were  accurately  summed  by  Aztec  electronic  counters  via  Singer  system  test 
rack  connectors. 

Tests  were  conducted  for  accelerometer-tc- IMU  case  orientations 
similar  to  those  used  during  the  IMU  calibration  tests.  Referring  to 
Figure  7-2,  platform  test  sequences  No.  I with  No.  3 and  No.  10  with 
No.  11  were  used  for  collecting  z-axis  velocity  pulse  data  in  both  the 
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up  and  down  orientations  after  leveling  r.ho  x and  y-accelerometers  in 
the  No.  3 and  No.  10  sequences.  Test  sequences  No.  5 and  6 were  used  to 
collect  data  for  the  y-accel erometer  in  the  up  orientation.  The 
y-accelerometer  data  similar  to  that  collected  by  the  Aztec,  electronic 
counters  is  shown  in  Figure  7-3  for  test  sequences  No.  5 and  6. 

The  x and  y-platform  axes  of  the  KT-73  I MU  are  instrumented  by  a 
two-axis  horizontal  accelerometer  employing  a force  rebalance  design. 

Test  results  indicated  that  the  number  of  velocity  pulses  generated  by 
this  accelerometer  changed  by  only  a fraction  of  a pulse,  out  of  a 
nominal  thirty  thousand  pulses  generated  in  a thirty  second  measurement 
interval,  for  the  various  platform  attitudes.  The  z-axis  accelerometer  is 
a single  axis  device  also  of  the  force  rebalance  type.  Test  results 
showed,  however,  that  the  total  velocity  pulse  count  varied  between  one 
and  four  pulses  out  of  a total  of  thirty  thousand  for  the  various  test 
attitudes.  Test  results  were  consistent  on  a day-to-day  basis  for 
specific  test  attitudes.  Test  results  also  appeared  independent  of  test 
time  with  the  total  number  of  pulses  generated  being  identical  for  the 
first  and  last  measurement  intervals  of  a particular  test  sequence. 

This  last  test  result  rules  out  a thermal  sensitivity  problem. 

Three  accelerometer  error  model  residual  states  had  previously  been 
included  in  the  Kalman  Filter  state  vector  to  account  for  a potential 
attitude  dependency.  These  states  had  been  modeled  as  random  bias 
accelerometer  scale  factor  errors  which  were  reset  to  their  initial 
uncertainty  at  the  beginning  of  each  platform  test  sequence.  These 
assumptions  are  consistent  with  the  above  test  results  and  the  only 
modification  being  required  is  an  increase  in  the  initial  uncertainty. 

The  above  data  suggests  attitude-dependent  scale  factor  errors  of  10  and 
30  ppm  for'  the  horizontal  and  vertical  accelerometers,  respectively. 

Increasing  the  initial  uncertainty  of  accelerometer  attitude- 
dependent  scale  factor  error  states  to  the  above  values  resulted  in  the 
non-orthogonality  error  angle  estimates  depicted  by  IMU  calibration  tests 
20-24  shown  in  Figure  8-6.  These  estimation  results  are  totally  consistent 
with  ensemble  averages  predicted  by  a Gaussian  distribution.  The  standard 
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deviation  of  the  estimation  error  for  the  u ^ non-orthogonality  error 
angle  predicted  by  the  filter's  covariance  matrix  is  4.3  arc  seconds. 
Accelerometer  attitude-dependent  scale  factor  estimates  are  shown  in 
Figure  8-8.  This  data  depicts  the  maximum  of  the  x,  y,  or  z-accelerometer 
error  for  each  platform  test  attitude. 

Modifying  the  Kalman  Filter  to  account  for  both  the  measurement 
timing  and  accelerometer  attitude-dependent  scale  factor  errors  also 
improved  the  repeatability  of  the  accelerometer  high  and  low  gain  mode 
bias  and  scale  factor  error  estimates.  The  Kalman  Filter  estimates  of 
these  error  sources  were  found  to  be  consistent  with  those  determined  via 
the  ASM-375  test  set.  The  repeatabi 1 i ty  of  the  estimates  was  also  found 
to  be  consistent  with  those  predicted  by  the  Singer  error  model. 

8.4  FUTURE  ERROR  MODEL  DEVELOPMENT  EFFORTS 

The  rationale  used  for  validating  both  the  IMU  calibration  Kalman 
Filter  and  error  model  has  largely  been  to  obtain  consistency  of  the 
estimation  results  with  the  known  KT-73  IMU  manufacturing  tolerances  and 
the  estimation  error  ensemble  averages  predicted  by  a Gaussian  distri- 
bution. Additional  confidence  has  been  obtained  by  cross-checking  the 
filter's  optimal  estimates  for  certain  IMU  error  sources  with  the 
calibration  results  obtained  via  Singer's  OCP  Auto-Cal  software  routines 
and  the  ASM-375  test  set.  The  IMU  error  source  software  compensation, 
heading  sensitivity,  and  thermal  and  magnetic  field  variation  special 
tests  have  provided  justification  for  some  of  the  major  assumptions  made 
while  cross-checking  baseline  IMU  calibration  test  results. 

The  present  status  of  the  IMU  calibration  Kalman  Filter  and  error 
model  is  such  that  there  is  no  known  deficiency  or  inconsistency  in  the 
calibration  results  presently  being  obtained.  This  does  not  imply  that 
undetected  IMU  error  sources  cannot  exist  or  might  not  become  evident 
under  more  exhaustive  validation  testing.  There  are  a number  of  error  model 
assumptions  (inertial  instrument  random  disturbances  modeling  is  a prime 
example)  of  fundamental  importance  in  the  IMU  calibration  methodology 
that  have  not  been  adeouately  verified.  Although  additional  studies  and 
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tests  are  considered  necessai j for  a more  complete  validation  and 
refinement  of  the  IMU  calibration  technique,  to  date  they  have  not  been 
conducted  largely  because  they  are  more  logically  and  conveniently 
performed  in  conjunction  with  the  aided-inertial  navigation  software 
development  and  evaluation  tests  discussed  in  Section  1.2. 

The  aided-inertial  navigation  software  development  and  evaluation 
tests  4,  5,  and  6 shown  in  Table  1-1,  are  expected  to  provide  a cross- 
check of  many  of  the  IMU  error  source  magnitudes  while  providing  a basis 
for  conducting  various  special  tests.  Gyro  scale  factor  and  misalignment 
angle  error  sources  can  be  accurately  estimated  by  the  Kerman  Filter 
formulated  for  the  tilt  table  static  navigation  tests.  This  test  and 
associated  filter  also  provide  a good  vehicle  for  studying  the  inertial 
instrument  random  disturbances  since  their  effects  can  be  continuously 
observed  over  several  hours.  The  MEL  perturbed  base  alignment  test  and 
associated  Kalman  Filter  will  provide  an  alternate  evaluation  of  the 
inertial  instrument  thermal  warm-up  transients.  Finally,  the  MEL  dynamic 
navigation  tests  and  associated  real-time  navigation  software  can  be  used 
to  demonstrate  the  software  compensation  of  calibrated  IMU  hardware  errors 
and  achievement  of  navigational  accuracy  which  will  be  limited  only  by 
alignment  and  system  random  disturbance  errors.  The  dat‘>  reduction/ 
reference  system  Kalman  Filter  formulated  for  the  MEL  navigation  tests 
will  be  capable  of  cross-checking  some  of  the  specific  force  dependent 
gyro  errors. 
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SECTION  IX 

I MU  CALIBRATION  RESULTS 

This  section  presents  the  major  IMU  calibration  results  obtained 
during  the  study.  These  results  are  twofold:  the  theoretical  estimation 
errors  as  obtained  from  the  Kalman  Filter's  covariance  matrix,  and  the 
filter's  estimates  for  the  magnitudes  of  the  IMU  and  associated  error 
sources.  These  results  are  conveyed  in  a pictorial  rather  than  tabular 
or  summary  manner  for  three  basic  reasons: 

(1)  This  type  of  presentation  greatly  enhances  one's  understanding 
of  the  calibration  methodology  and  associated  Kalman  Filter's  operation. 

It  also  brings  out  many  subtle  points  not  discussed  in  the  text. 

(2)  It  gives  the  reader  an  opportunity  to  view  the  consistency 
obtained  between  the  Kalman  Filter's  estimates  and  their  predicted 
uncertainty  thus  giving  credence  to  large  dimensional  Kalman  Filters  as 
viable  mathematical  tools. 

(3)  It  can  be  a valuable  basis  for  debugging  and  validating  a 
Kalatan  Filter  computer  program  for  a similar  application,  in  spite  of 
its  characteristic  of  seeming  to  be  repetitious 

Estimation  accuracies  obtained  from  the  filter's  updated  covariance 
matrix  are  shown  for  each  IMU  measurement  taken  and  plotted  as  a function 
of  time.  Platform  test,  orientations  are  superimposed  on  these  plots  t.o 
facilitate  one's  understanding  of  the  error  source  sepcr- Lion  process. 
Kalman  niter  estimates  of  the  same  error  source  are  plotted  in  groups  of 
five  and  correspond  to  the  last  five  ( 20  through  24)  IMU  calibration  tests 
conducted.  Five  data  points  were  selected  as  a compromise  between 
sufficient  data  demonstrating  estimation  consistency  and  data  reduction 
complexi ty . 
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9.1  SYSTEM-LEVEL  IMU  MEASUREMENTS 

Estimation  accuracy  for  the  system-level  IMU  differential  gimbal 
angle  and  pseudo  velocity  error  measurement  states  is  shown  in  Figure  9-1 
A typical  record  of  recorded  pseudo  velocity  pulses  was  shown  in 
Figure  7-3.  Platform  test  sequences  and  data  recording  timing  instances 
and  intervals  are  the  same  in  both  figures.  Total  measurement  time 
during  a calibration  test  is  68  minutes.  The  platform  test  sequences 
have  been  arranged  to  minimize  the  time  required  to  manually  torque  the 
platform  from  one  test  orientation  to  another.  Depending  on  the  test 
sequence  number,  the  30  and  90  degree  per  minute  x-y  and  z-gyro  slew 
rates  result  in  up  to  several  minutes  being  required  to  re-orientate 
the  platform.  Total  time  required  to  conduct  the  IMU  calibration  test, 
including  a seven  minute  warm-up  period  oefore  any  measurements  are 
taken,  is  nominally  two  to  two  and  a half  hours.  The  platform  test 
attitudes  shown  across  the  top  of  these  figures  are  the  same  as  those 
depicted  in  Figure  7-2  with  the  first  attitude  representing  the  north, 
west,  and  up  navigation  frame. 

The  initial  uncertainty  of  both  the  differential  gimbal  angle  arid 
pseudo  velocity  error  is  zero  at  the  start  of  each  test  sequence  due  to 
the  manner  in  which  these  error  states  were  defined.  The  Calcomp  plotter 
logic  used  in  generating  the  semi-log  graphs  was  to  omit  plotting  a 
30  second  propagated  uncertainty  segment  when  its  initial  value  was  less 
than  the  minimum  value  indicated  on  the  ordinate  axis.  This  logic, 
necessitated  by  the  semi-log  plots,  prevented  plotting  the  initial  steep 
rise  from  the  zero  initial  uncertainty. 

To  distinguish  the  different  sets  of  data  in  these  figures,  alphabet 
characters  are  plotted  in  the  center  of  the  dead  time  between  the 
platform  test  sequences.  Tie  vertical  positioning  of  tnese  characters 
is  determined  by  the  final  value  of  the  updated  covariance  matrix  for 
the  particular  test  sequences.  Referrring  to  the  differential  gimbal 
angle  plot  in  Figure  9-1,  only  orie  differential  synchro  angle  measurement 
is  taken  during  the  entire  test.  This  measurement  corresponds  to  the 
azimuth  gimbal  axis  at  the  end  of  the  first  test  sequence.  Estimation 
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accuracy  for  the  differential  azimuth  gimbal  angle  immediately  3ft.er  the 
measurement  is  taken  is  essentially  that  of  the  synchro  device's  precision 
of  5 arc  seconds  as  indicated  in  the  figure.  Estimation  of  the  differential 
giaibal  ingles  for  all  instances  other  than  the  above  measurement  is 
accomplished  by  the  Kalman  Filter  on  tne  basis  of  the  pseudo  velocity 
measu cements . 

In  contrast  to  the  gimbal  angle  measurements,  three  axis  pseudo 
velocity  measurements  are  utilized  by  the  Kalman  Filter  at  each 
measurement  interval.  The  accuracy  of  these  measurements  is  determined 
by  the  velocity  pulse  quantization  level  which  is  governed  by  the 
accelerometers  scale  factor.  The  accelerometer  low  gain  navigation  scale 
factor  is  nominally  used  throughout  the  test.  Tfie  high  gain  alignment 
scale  factor  is,  however,  used  by  the  x and  y-axis  accelerometers  for  the 
entire  first  and  last  minute  and  a half  of  the  third  and  tenth  test 
sequence.  Pseudo  velocity  measurement  accuracy  is  two  orders  of 
magnitude  greater  for  the  high  gain  rather  than  tne  low  gain  mode.  A 
thirty  second  measurement  intervu'  is  used  throughout  the  test  except 
during  the  first  test  sequence  where  it  is  ten  seconds. 

Differential  roll  gimbal  ang I e magni tude  estimation  results  are 
shown  in  Figure  9-2.  This  data  corresponds  to  the  total  roll  gimbal 
angular  excursion  occurring  during  each  of  the  twelve  platform  test 
sequences  and  was  obtained  as  the  final  Kalman  Filter  estimate  or  this 
error  state  for  each  test  sequence.  Superimposed  about  the  average 
estimate  for  each  test  sequence,  not  shown,  are  one  standard  deviation 
estimation  error  brackets  for  a total  bracket  separation  of  two  standard 
deviations.  These  estimation  error  standard  deviations  are  the  same  as 
those  corresponding  to  the  end  of  a test  sequence  shown  in  Figure  9-1. 

Figure  9-2  shows  that  the  roll  gimbal  angle  excursions  for  the  last  five 
calibration  tests  conducted  are  highly  consistent,  repeating  to  almost 
the  accuracy  of  the  estimate.  Note  that  test-to-test  variations  are  due 
largely  to  gyro  bias,  spin  axis  mass  unbalance  and  quadrature  c,  tu^n-on 
and  platform  positioning  variations  and  are  not  required  to  obey  the 
ensemble  averages  predicted  by  the  Gaussian  distribution  for  ■fhe 
estimation  error. 
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9.2  COMPOSITE  DRIFT  RATE  ERRORS 

Platform  and  gyro  total  drift  rate  error  estimation  accuracy  is 
shown  in  Figure  9-3  and  the  correspond!' ng  estimated  magnitudes  are 
shown  in  Figures  9-4  and  9-5,  respectively.  These  quantities  are 
composite  error  sources  and  are  estimated,  and  the  estimation  uncertainty 
is  determined  as  described  in  Section  8.2.3.  Gyro  total  drift  rate 
error  contains  all  the  error  sources  composing  the  platform  total  drift 
rate  error  except  the  gyro  scale  factor  and  misalignment  angle,  IMU 
case-to-navigation  frame  misalignment  angle  arid  the  gimbal  differential 
and  synchro  angle  bias  errors.  Both  the  absolute  value  and  estimation 
accuracy  of  these  composite  error  sources  are  of  general  interest  in 
the  development  of  Kalman  Filter  IMU  calibration  computer  programs. 

Estimation  of  platform  drift,  rate  can  be  inferred  direcfy  from  the 
accel erometer  measurements,  except  when  an  axis  is  orientated  near  vertical, 
resulting  in  the  relatively  high  estimation  accuracy  shown  at  the 
completion  of  a test  sequence.  Final  test  sequence  estimation  accuracy 
is  roughly  0.006  degrees  per  hour  and  could  be  increased  to  the  limiting 
short  correlation  time  gyro  random  disturbance  magnitude  of  0.003  degree 
per  hour  by  doubling  the  measurement  time  to  ten  minutes  per  test  sequence. 
Gyro  total  drift  rate  estimation  accuracy  is  somewhat  less  than  that  for 
the  platform  due  to  the  incomplete  separation  of  this  composite  error 
source  from  some  of  the  total  platform  drift  rate  recited  error  sources. 
Examination  of  the  estimated  magnitudes  for  these  composite  errors 
depicted  in  Figures  9-4  and  9-5  shows  that  they  are  repeatable  to  roughly 
the  estimation  error. 

9.3  INERTIAL  INSTRUMENT  WARM-UP  TRANSIENTS 

Inertial  instrument  warm-up  transients  are  not  generally  considered 
as  calibration  parameters  and  were  not  included  in  the  first  version  of 
the  Kalman  Filter.  During  the  course  of  the  development  effort, 
excessive  navigational  errors  encountered  witn  the  inertial  navigation 
systems  on  the  A-7D  aircraft  were  postulated  as  caused  by  an  incomplete 
alignment  due  to  instrument  thermal  warm-up  transients.  For  this  reason 


125 


AFAL-TR-77-,'5 


PLATFORM  TOTAL  DRIFT 
CALIBRATION  RESULTS 


PLATFORM  TOTAL  DRIFT 
CALIBRATION  RESULTS 


1.6 


j 

3 


127 


R-//-/5 


GYRO  Torn  ORff-r 
CRl.  I BRflT  I ON  RESULTS 


GYRO  TOTAL  DR  IF! 
CRL I BRRT I ON  RESULTS 


1 . 6 


1 . 2 


0.8 


0.4  I' 


o;  ct 
oi 
a:  \ 
0:0 
liJLU 


-0.4  -1 


H 

-0.8  5 


.] 

- 1 . 2 4 


-1.6  ’2 


1.6 


1 . 2 : 


G.e 


0.  4 d 


tx:  0: 

ox 

CK'v 

ItO 

UJIU 


o.o  d 


a - 


~0. 4 j 


-0.  0 4 


1.2  A 


■1.6  ^ 


PT  POS-2 


! PT  POS-1 

o.o.LitttLL 


PT  FOS-.i 


GYRO  TOTAL  DRIFT 
CAL I3RAT ION  RESULTS 


PT  POS-8 


Pi  POS-7 


1 . 6 


1 . 2 


0.8 


0.4 


GYRO  TOTAL  DRIFT 
CAI  I BRA T I ON  RESULTS 


1 . 6 


1.2 


0.8 


0.4 


PT  POS-12 


c:  :v 
tit 

0:0 

bJLU 


PT  P0S-10 


0.0  L 


cc 
cc 
< n 


■0.4.1 


-0.8  A 


-1.2.3 


■1.6 


PT  POS-11 


.1 


Figure  9-5.  Gyro  Total  Drift  Rate  Estimation  Results 


and  in  preparation  for  a demonstration  program  with  AGMC,  instrument 
warm-up  error  states  were  included  in  the  second  version  of  the  Kalman 
Filter. 

Gyro  thermal  drift  rate  error  modeling  and  fine  gyro-compass  align- 
ment timing  considerations  are  shown  in  Figure  9-6.  The  gyro  warm-up 
characteristic  consists  of  an  exponentially  decaying  bias  term  which 
can  be  assumed  to  be  the  envelope  of  a complicated  transient  function 
of  the  ambient  temperature  during  the  rapid  warm-up  period.  The  nominal 
data  depicted  in  Figure  9-6  corresponds  to  a 0°  C ambient  temperature. 

A similar  characteristic  is  obtained  for  a 35°  C ambient  temperature 
except  that  the  transient  magnitude  over  the  fine  alignment  time  interval 
is  greatly  reduced.  To  have  a negligible  effect  on  navigation  accuracy, 
instrument  thermal  transients  must  die  out  by  the  end  of  the  fine  alignment 
time  interval. 

The  inertial  instrument  exponentially  decaying  bias  modeling  has 
been  described  in  Section  3.  It  was  also  necessary  to  model  the  fine 
grain  structure  of  the  thermal  transient  to  ensure  proper  filter 
performance.  It  was  assumed  that  the  fine  grain  component  of  the  thermal 
transient  could  be  modeled  as  a random  bias  error  over  two  short  ten 
second  measurement  intervals.  The  gyro  error  model  residual  states  were 
used  to  implement  this  modeling  and  their  uncertainty  reset  to  its  initial 
value  every  twenty  seconds. 

Inertial  instrument  exponentially  decaying  bias  thermal  transient 
estimation  accuracy  is  shown  in  Figure  9-7.  These  transient  errors 
decay  rapidly  to  insignificant  levels  and  need  not  be  considered  after 
the  first  test  sequence.  The  uncertainty  of  these  er^or  states  decreases 
by  virtue  of  both  the  measurements  made  and  the  exponential  time  decay  of 
their  magnitudes.  The  decrease  in  the  x,  y-gyro  drift  rate  error- 
uncertainty  is  a factor  3.04  and  1.7  as  a result  of  the  exponential  decay 
and  measurements  taken  respectively  at  the  end  of  the  ten  minute  test 
sequence.  Corresponding  data  for  the  accelerometers  are  exponential 
decay  and  measurement  factors  of  54.5  and  only  1.24  for  the  x,y-accelerometers. 
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Inertia;  instrument  estimated  thermal  transient  magnitudes  are 
shown  in  Figure  9-8.  This  data  corresponds  to  a time  instant  seven 
minutes  after  power  turn-on  which  is  the  start  of  the  fine  gyro-compass 
alignment  time  interval.  This  data  was  obtained  by  referring  the  Kalman 
Filter  estimates  obtained  at  seventeen  minutes  after  power  turn-on, 
seven  minutes  warm-up  plus  ten  minutes  measurement  time,  backwards  in 
time  ten  minutes.  The  brackets  superimposed  about  the  average  estimates 
correspond  to  warm-up  transient  uncertainty  at  seven  minutes  after 
power  turn-on.  These  uncertainties  were  obtained  by  referring  the 
seventeen  minutes  after  power  turn-on  Kalman  Filter  estimation  errors 
backwards  in  time.  Table  9-1  lists  the  seven-minute  after-power  turn-on 
thermal  warm-up  transient-uncertainty  reduction  attributable  to  the 
Kalman  Filter. 


TABLE  9-1 

INERTIAL  INSTRUMENT  THERMAL  TRANSIENT 
UNCER1 AINTY  AT  SEVEN  MINUTES  AFTER  POWER  TURN  ON 


Inertial  Instrument 


x,  y-Gyro 
z-Gyro 

x,  y-Accelerometer 
z-Accel erometer 


Standard  Deviation  of  Thermal  Transient 
Uncertainty 

Initial  After  Kalman  Filtering 

0.014°/hr  0.0084°/hr 

0.024°/hr  0.024°/hr 


5.0  pg  4.0  pg 

200.0  pg  2,4  pg 


Table  9-1  shows  that  the  Kalman  Filter  is  capable  of  making  a 
dramatic  initial  uncertainty  reduction  for  only  the  z -acclerometer  inertial 
instrument.  This  is  possible  because  the  warm-up  transient  is  the  only 
error  source  that  gives  rise  to  the  z-accel erometer 's  output  showing  a 
time  dependency  when  oriented  vertically.  Recovered  transient  magnitudes 
of  -114,  -42,  -339,  -414,  and  -42  pg  indicate  a rather  large  turn-on 
variation  for  this  instrument. 
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A conclusion  drawn  from  the  inertial  instrument  warm-up  transient 
estimation  results  is  that  the  magnitudes  of  these  errors  decay  to 
roughly  the  instrument  noise  level  at  the  end  of  the  fine  gyro-compass 
alignment  time  interval  and  should  not,  therefore,  have  a significant 
effect  on  alignment  accuracy.  This  conclusion  applies  only  to  the 
particular  KT-73  i MU  tested  and  an  ambient  temperature  of  23°C.  While 
of  general  interest,  no  tests  were  conducted  at  other  than  with  a 23°C 
ambient  temperature. 

9.4  GYRO  ERROR  SOURCES 

Gyro  error  sources  are  considered  in  two  groups:  (1)  the  random 
disturbance  errors  which  ultimately  limit  achievable  navigation  accuracy, 
and  (2)  the  errors  commonly  termed  calibration  parameters  whose  effects 
can  be  software  compensated  for  in  the  operational  flight  computer  program. 
Gyro  instrument  angular  alignment  errors  are  discussed  in  Section  9.6.1. 

9.4.1  Random  Errors 

Gyro  random  error  estimation  accuracy  is  shown  in  Figure  9-9.  The 
platform  test  orientations  and  measurement  times  have  been  formulated  so 
that  the  Kalman  Filter  estimates  the  gyro  calibration  parameters  to  an 
accuracy  approaching  the  limit  imposed  by  these  random  errors.  No  attempt 
is  made  to  estimate  these  random  errors  by  increasing  the  IMU  measurement 
time.  Although  the  filter  does  not  significantly  reduce  the  uncertainty 
of  these  random  errors,  inclusion  of  them  in  the  filter  is  required  for 
accurate  magnitude  estimation  and  estimation  accuracy  results. 

Gyro  attitude  dependent  bias  drift  rate  error  magnitude  estimates 
are  shown  in  Figure  9-10  for  each  of  the  platform  test  orientations. 

The  estimates  for  these  error  model  residual  states  are  similar  to  those 
shown  in  Figure  8-3  except  for  being  plotted  to  a different  scale. 

Modeling  of  the  discovered  quadrature  g error  source  in  the  Kalman  Filter 
state  vector  resulted  in  a reduction  of  the  magnitude  estimates  for  these 
errors  to  the  heading  sensitive  drift  rate  error  hardware  specification 
of  0.004  degree  per  hour. 
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9.4.?  Ca  1 ibratioi)  Parameters 

Gvro  bias  drift  rate  and  scale  factor  error  estimation  accuracies 
are  shown  in  F i y u re  9-11  and  corresponding  magnitude  estimates  are  given 
in  Figure  9-12.  The  initial  uncertainty  of  the  error  sources  given  in 
Reference  1 were  not  always  found  to  be  consistent  with  those  recovered 
from  the  1 MU  tests.  In  some  instances,  these  inconsistencies  were  due 
to  discrepancies  in  the  definition  of  the  error  sources.  A good 
example  of  this  is  z-ayro  bias,  to  which  Reference  1 assigns  an  uncertainty 
of  0.02  degrees  per  hour,  implying  a turn-on  variation.  The  actual 
hardware  configuration  is  such  that  the  2-gyro  bias  is  adjusted  to  offset 
the  drift  rate  error  produced  by  the  spin  axis  mass  unbalance  when  the 
z-gyro  Input  axis  is  oriented  vertically.  This  procedure  requires  that 
the  bias  error  magnitude  be  approximately  that  of  the  spin  axis  mass 
imbalance  drift  rate  error  in  a one  g field  or  0.5  degree  per  hour. 

Regardless  of  the  reason  for  the  discrepancies  between  the  initial 
uncertainties  given  by  Reference  1 arid  those  implied  by  the  Kalman  Filter 
error  magnitude  estimates,  the  initial  covariance  matrix  was  adjusted 
so  that  its  initial  variances  were  consistent  with  the  end-of-test 
Kalman  Filter  magnitude  estimates.  This  procedure  worked  ''ell  with  no 
apparent  pitfalls;  however,  it  should  he  understood  that  the  initial  I MU 
error  source  t ’certainties  shown  by  the  estimation  accuracy  plots  of  this 
section  correspond  to  the  particular  KT-73  IMU  tested  and  do  not  represent 
an  IMU  error  model  derived  by  ensemble  IMU  testing. 

Reference  to  the  trends  of  and  final  estimation  accuracies  shown  in 
Figure  9-11  indicates  that  estimation  accuracies  could  be  further 
reduced  by  increasing  the  length  of  the  test..  A covariance  analysis 
showed  that  the  gyro  bias  drift  rate  estimation  errors  could  be  further 
reduced  to  the  limiting  0.005  degree  per  hour  by  adding  a couple  of 
additional  platform  test  att’+udes  and  doubling  the  measurement  time. 

This  was  felt  to  be  unnecessary  since  the  calibration  accuracy  obtained 
by  the  present  test  schedule  is  consistent  with  the  one  to  three  MM 
inertial  navigation  performance  of  the  KT-73  IMU.  Note  that  the  turn-on 
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gyro  bias  drift  rate  error  variations  shown  in  Figure  9-12  are  consistent 
witn  the  0.01  degree  per  hour  gyro  performance  claimed  for  the  system. 


Specific  force  dependent  gyro  drift  rate  error  estimation  accuracies 
are  shown  in  Figures  9-13  and  9-15  and  corresponding  magnitude  estimates 
shown  in  Figures  9-14  and  9-16.  The  input  axis  mass  unbalance  errors 
appears  to  be  insignificant  error  sources  for  the  I MU  while  the  recovered 
spin  axis  mass  unbalance  errors  were  found  to  be  a factor  of  three  larger 
than  given  in  Reference  1.  The  g-squared  gyro  anisoela.stic  error  also 
appears  to  be  insignificant  while  the  discovered  quadrature  g errors  can 
be  significant  for  high  velocity  and  dynamic  flight  profiles.  Note  that 
all  the  specific  force-dependent  gvro  drift  rate  errors  are  repeatable  to 
within  the  estimation  accuracy  of  the  Kalman  Filter.  No  significant  shift 
in  any  of  these  errors  was  observed  in  over  400  hours  of  I MU  operation 
accumulated  in  a year's  time.  Software  compensation  for  these  error 
sources  would  be  a straightforward  procedure. 


The  symmetry  of  construction  between  the  x and  y axes  of  the  two- 
degree-of-freedom  x,  y-gyro  provides  a ready  means  for  cross  checking 
some  of  the  Kalman  Filter  estimates.  Since  the  x and  y channels  of 
this  gyro  share  the  same  spin  axis,  mass  unbalance  of  this  axis  should 
produce  identical  mass  unbalance  drift  rate  error  coefficients  for  these 
channels.  Urisymmetrical  machining  of  the  rotating  flexure  suspension 
joint  should  also  give  rise  to  equal  quadrature  g drift  rate  error 
coefficients  for  these  channels.  Reference  to  Figures  9-14  and  9-16 
shows  that  the  filter  has  estimated  these  error  sources  to  be  of  the 
same  magnitude  for  each  channel  to  within  the  estimation  accuracy. 
Inclusion  of  the  high  degree  of  correlation  in  the  initial  covariance 
matrix  does  not,  however,  substantially  increase  the  estimation  accuracy 
for  these  error  stares. 


9.5  ACCcLEROMETER  ERROR  SOURCES 

The  accelerometer  error  sources  are  covered  in  a similar  manner  to 
those  for  the  gyro.  They  are  grouped  into  random  disturbances  and  the 
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commonly  denoted  calibration  parameters  with  the  accelerometer 
instrument  angular  misalignments  being  treated  in  Section  9.6.1. 

9.5.1  Random  Errors 

Accelerometer  random  error  estimation  accuracy  is  shown  in 
Figure  9-17  and  the  correspond! ng  magnitude  estimates  depicted  in 
Figures  9-18  and  9-19.  The  attitude-dependent  scale  factor  errors 
represent  the  accelerometer  error  model  residual  states  and  are  the 
major  random  errors  limiting  accel erometer  parameter  calibration  accuracy. 
These  errors  are  thought  to  be  related  to  magnetic  field  variations 
with  platform-to-IMU  case  attitude  changes,  since  they  do  not  show  a 
time  dependency  as  would  be  the  case  if  they  were  temperature  related. 

The  clock  stability  error  refers  to  the  real-time  software  measurement 
timing  error  and  is  modeled  as  a scale  factor  error  which  affects  ail 
three  accelerometers.  Although  the  filter  is  unable  to  significantly 
reduce  the  uncertainty  of  these  errors,  inclusion  of  their  effects  in  the 
filter  algorithm  is  essential  to  obtaining  proper  estimation  results. 

The  attitude-dependent  scale  factor  error  magnitude  estimates  shown 
in  Figure  9-18  represent  that  set  having  the  largest  recovered  value 
for  a particular  test  attitude.  Comparing  these  estimates  with  the 
platform  test  attitudes  shown  in  Figure  7-2  shows  that  discrepancies 
appear  only  after  a particular  accelerometer's  input  axis  has  been 
oriented  at  or  near  vertical  for  the  second  time,  orientations  both  up 
or  down,  as  would  be  expected.  Comparing  the  measurement  timing  with  the 
attitude-dependent  scale  factor  error  magnitude  estimates  shows  them  to 
be  of  generally  reduced  magnitudes.  Since  the  real-time  software  timing 
errors  do  not  appreciably  affect  I MU  calibration  accuracy,  as  compared 
to  the  attitude-dependent  scale  factor  errors,  efforts  to  improve  the 
real-time  data  collection  software  were  suspended. 

9.5.2  Calibration  Parameters 

Accelerometer  low  and  high  gain  mode  bias  and  scale  factor  error 
estimation  accuracies  are  shown  in  Figure  9-20  and  9-22  and  the 
corresponding  magnitude  estimates  shown  in  Figures  9-21  and  9-23. 
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Low  gain  mode  bias  and  scale  factor  errors  can  be  separated  and 
accurately  observed  by  the  Kalman  Filter  after  an  accelerometer's 
input  axis  has  bec-n  oriented  up  and  down  in  the  gravity  field.  Error 
source  separation  is  possible  using  these  attitudes  because  no  other 
I MU  errors,  except  the  low  level  random  disturbances,  have  a first-order 
effect  on  the  accelerometer's  output  when  its  input  axis  is  parallel  to 
the  gravity  vector. 

The  above  technique  cannot  be  used  to  separate  the  alignment  high 
gam  mode  bias  and  scale  factor  errors  because  the  capacitive  reset  inte- 
grators (CAPRIS)  saturate  at  a pulse  repetition  frequency  of  ten  thousand 
pulses  per  second,  which  occurs  when  an  accelerometer's  input  axis  is 
tilted  about  six  degrees  off  horizontal.  Calibration  of  the  high  gain  mode 
bias  and  scale  factor  errors,  therefore,  requires  extremely  accurate 
accelerometer  attitude  information,  including  the  angular  drift  over  the 
measurement  interval.  The  philosophy  used  in  calibrating  the  high  gain 
mode  errors  is  to  use  two  platform  test  attitudes  where  the  accelerometer's 
input  axes  are  oriented  both  up  and  down  roughly  five  degrees  off 
horizontal.  To  allow  both  an  accurate  accelerometer  attitude  and 
attitude  rate  determination , the  Kalman  Filter  employs  low  gain  acceler- 
ometer measurements  over  the  first  five  minutes  of  the  platform  test 
sequences,  which  are  No.  3 and  No.  10.  The  accelerometer  mode  is  then 
switched  from  low  to  high  gain  and  high  gain  measurements  made  for  an 
additional  minute  after  a thirty  second  switching  delay.  Complete 
separation  of  the  high  gain  errors  must,  however,  wait  until  the  filter 
is  able  to  separate  the  low  gain  errors  during  test  sequence  No.  12. 

This  is  an  efficient  procedure  for  estimating  the  high  gain  errors  because 
the  two  test  attitudes  used  were  perturbations  of  existing  ones  and  not 
additionally  added. 

Examination  of  the  bias  and  scale  factor  error  mayiiitude  estimates 
shown  in  Figures  9-21  and  9-23  indicates  that  the  repeatability  of  the 
estimate  is  on  the  order  of  the  estimation  error.  For  example,  x- 
accelerometer  low  gain  bias  estimates  fluctuate  about  20  to  25  pg  about, 
the  average  estimate,  which  is  a factor  of  3 to  4 times  the  estimation 
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error  standard  deviation.  The  y-acce'lerometer  estimate  variations  are 
somewhat  less,  being  only  1 to  2 times  the  estimation  error  standard 
deviation.  Both  variations  are  well  within  the  manufacturer's 
Reference  1 value  of  67  pg.  Whether  these  variations  are  due  to 
instrument  turn-on  fluctuations  or  measurement  error  or  combinations  of 
both  is  not  known.  Estimation  repeatability  on  the  order  of  10  p.g  was 
only  obtained  after  adding  the  random  errors  of  Section  9.5.1  to  t.he 
Kalman  Filter  state  vector. 

9.6  ANGULAR  ERROR  SOURCES 

Angular  error  sources  estimated  by  the  Kalman  Filter  fall  into  two 
classes:  inertial  instrument-to-platform  frame  misalignment  errors, 
and  angular  errors  associated  with  orienting  the  platform  relative  to 
the  navigation  frame. 

9.6.1  Inertial  Instruments 

Accelerometer  and  gyro- to- platform  frame  misalignment  angle 
estimation  accuracy  is  shown  in  Figure  9-24  and  the  corresponding 
magnitude  estimates  given  in  Figure  9-25.  As  discussed  in  Section  8.3, 
inertial  instrument  misalignment  angle  errors  relative  to  the  navigation 
frame  are  the  quantities  of  fundamental  importance  in  the  calibration 
methodology.  The  inertial  instrument-to-navigation  frame  misalignment 
angles  are  composite  error  sources  composed  of  inertial  instrument-to- 
platform,  synchro  differential  and  bias  and  I MU  case-to-navigation  frame 
angle  errors.  The  Kalman  Filter  is  precluded  from  accurately  estimating 
the  inertial  instrument-to-platform  misalignment  angle  errors  because  it 
cannot  separate  them  from  the  synchro  angle  bias  errors  which  assume  new 
values  each  time  the  platform  is  reoriented. 

In  an  attempt  to  overcome  this  estimation  accuracy  limitation,  gyro 
and  accelerometer  misalignment  angle  errors  have  been  reformulated  to 
represent  instrument  misalignments  from  orthogonal  sets  and  the  mis- 
alignment of  the  orthogonal  instrument  sets  from  each  other  rather  than 
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the  platform.  The  motivation  for  this  approach  is  to  let  the  inertial 
instruments  serve  as  relative  references  rather  than  to  use  the  gravity 
vector  as  an  absolute  reference. 

Accelerometer  and  gyro  instrument  set  non-orthogorial  i ty  angle  error 
estimation  accuracies  are  shown  in  Figure  9-26  and  the  corresponding 
magnitude  estimates  given  in  Figure  9-27.  The  estimation  accuracy  of  the 
misalignment  between  the  "squared-up"  accelerometer  and  gyro  instrument 
sets  is  shown  in  Figure  9-28  and  the  corresponding  magnitude  estimates 
given  in  Figure  9-29.  The  misalignment  of  the  instrument  sets  from  each 
other  is  approximated  by  a small  angle  orthogonal  coordinate  transformation 
matrix. 

Figure  9-26  shows  that  the  Kalman  Filter  is  capable  of  calibrating 
the  accelerometer  non-orthogonality  error  angles  with  an  error  of  only 
a few  arc  seconds  by  comparing  accelerometer  outputs  against  each  other. 
Since  the  gyros  do  not  have  an  output  analogous  to  the  accelerometers, 
the  filter  has  to  relay  on  the  equivalent  platform  drift  rate  errors 
produced  by  the  non-orthogonality  angle  errors  to  separate  these  angular 
errors  from  the  other  gyro  error  sources.  The  equivalent  drift  rate 
errors  produced  by  the  one  arc  minute  instrument  misalignments  are, 
however,  only  on  the  order  of  the  short  and  long  time  constant  gyro 
random  disturbance  of  0.003  and  0.005  degree  per  hour,  which  prevents 
the  filter  from  reducing  the  uncertainties  associated  with  these 
misalignments.  This  is  an  acceptable  filter  performance,  since  the  one 
arc  minute  gyro  misalignments  do  not  contribute  significantly  to 
KT-73  I MU  navigational  errors. 

Figure  9-27  illustrates  that  the  accelerometer  non-orthogonality 
angle  error  magnitude  estimates  are  consistent  with  Gaussian  distribution 
ensemble  averages.  The  estimation  accuracy  for  the  3^z  angle  error 
component  is  somewhat  greater  than  the  other  two  because  it  is  not 
affected  by  the  larger  attitude-dependent  scale  factor  errors  associated 
with  the  vertical  accelerometer.  Accurate  calibration  of  this  error 
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angle  can  be  important  to  navigation  accuracy  since,  if  software 
compensated  for,  it  can  ensure  that  the  x and  y-accel erometers  are 
orthogonal  in  the  horizontal  plane. 

9.6.2  Platform  Attitude  Related 

The  accuracy  with  which  the  platform  can  be  oriented  at  the  start 
of  each  test  sequence  is  determined  by  two  major  error  sources:  the 
misalignment  errors  of  the  IMU  case  relative  to  the  navigation  frame, 
and  the  errors  in  positioning  the  gimbals  (termed  the  synchro  angle 
bias  errors).  The  estimation  accuracy  for  these  angular  errors  is 
given  in  Figure  9-30  and  the  corresponding  magnitude  estimates  shown 
in  Figures  9-31  and  9-32. 

The  initial  uncertainty  of  the  roll,  pitch,  and  azimuth  control 
transmitter  and  transformer  synchro  devices  is  six  arc  minutes  and  is 
reset  to  this  value  in  t lie  covariance  matrix  at  the  start  of  each 
platform  test  sequence.  Figure  9-30  shows  that  the  filter  can  reduce 
this  uncertainty  to  about  one  arc  minute  with  the  first  pseudo  velocity 
error  measurement  for  the  latter  test  sequences.  End-of-test  one  arc 
minute  estimation  accuracy  is  also  obtained  for  the  IMU  case-to-navi gation 
frame  misalignment  angle  errors.  The  Kalman  Filter  reduces  the  uncertainty 
associated  with  these  error  sources  and  the  inertial  instrument-to-pl atform 
misalignment  angle  errors  by  statistical  averaging  rather  than  by 
separating  and  accurately  observing  them.  The  separation  process  is,  as 
previously  discussed,  precluded  by  the  introduction  of  new  synchro  angle 
bias  errors  each  time  the  platform  is  reoriented. 

Recovered  roll  synchro  angle  bias  magnitude  estimates  shown  in 
Figure  9-31  conform  to  Gaussian  distribution  ensemble  averages.  To 
cross-check  the  Kalman  Filter's  estimates  of  these  errors,  the  synchro 
devices  were  calibrated  by  making  large  angle  IMU  case  rotations  relative 
to  the  navigation  frame  and  using  the  rotary  tilt  table  with  the  IMU 
in  the  inertial  navigational  mode.  Calibration  results  for  the  rotary 
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tilt  table  and  Kalman  Filter  methods  were  consistent  to  the  estimation 
accuracy  predicted  by  the  filter.  As  an  example,  the  roll  synchro  angle 
bias  for  the  tenth  test  sequence  was  determined  as  -10.2  arc  minutes 
using  the  rotary  tilt  table. 
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SECTION  X 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  objective  of  this  study  has  been  the  design,  development,  analysis 
and  laboratory  demonstration  of  a system- 1 evel  I MU  calibration  technique. 

In  addition  to  the  calibration  function,  it  was  desirable  that  the 
formulated  methodology  afford  some  capaoility  for  conducting  I MU  error 
model  development  arid  validation  efforts.  Major  study  ground  rules 
were  that:  only  system- level  IMU  measurements  be  used  such  that  the  IMU 
would  not  have  to  be  disassembled , no  AGE  other  than  available  laboratory 
equipment  be  required,  and  the  solution  be  implementable  on  a general 
purpose  digital  computer. 

A calibration  methodology  was  formulated  on  the  basis  of  the  essential 
functions  an  IMU  performs:  realization  of  a reference  platform  coordinate 
frame  and  measurement  of  the  instantaneous  specific  force  vector. 

Physical  vector  equations  were  derived  which  embodied  the  methodology 
premises.  The  physical  vector  equations  were  manipulated  to  yield  a 
system  of  68  simultaneous  stochastic  differential  equations.  A 68-state 
Kalman  Filter  was  implemented  for  the  solution  of  this  system  of 
equations  and  optimal  estimation  of  the  IMU  calibration  parameters.  A 
computer  simulation  of  the  IMU,  Kalman  Filter,  and  laboratory  test 
procedure  indicated  the  approach  to  be  sound. 

Technical  feasibility  of  the  methodology  was  demonstrated  through 
laboratory  nesting  and  the  actual  calibration  of  a KT-73  IMU,  using  an 
off-line  data  processing  computer  program.  The  suitability  of  the 
methodology  and  associated  large  dimensional  Kalman  Filter  for  conducting 
IMU  error  model  development  and  validation  efforts  was  demonstrated  by 
the  detection,  isolation,  and  identification  of  previously  unmodeled  IMU 
error  sources.  The  experimental  calibration  results  show  complete 
consistency  with  the  statistical  ensemble  averages  predicted  by  the 
Kalman  Filter’s  covariance  matrix  and  previous  calibration  parameter 
values  obtained  via  the  ASM-375  test  set. 
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It.  is  coricl  uded  that  the  developed  I MU  calibration  technique  is 
ideally  suited  to  aid  the  Reference  Systems  Software  Group  in  the 
development  of  aided-inertial  navigation  real-time  software.  Periodic 
calibration,  of  the  IMU  will  provide  a high  degree  of  confidence  that 
the  hardware  is  operating  properly.  Software  compensation  for  the  static 
IMU  error  sources  in  the  operation  flight  computer  program  will  enhance 
the  navigational  accuracy  of  the  system  while  effectively  removing  many 
of  these  error  sources  from  consideration  in  the  design  of  the  suboptimal 
navigation  Kalman  Filter.  And  finally,  refinement  and  validation  of  the 
IMU's  error  model  by  use  of  the  calibration  Kalman  Til  ter  and  associated 
methodology  will  provide  assurance  that  the  "to  be  developed"  navigation 
Kalman  Filter  solves  the  real  problem  at  hand. 

DOD  agencies  have  observed  that  advanced  avionics  research  projects 
sometimes  reveal  techniques  that  could  result  in  significant  operational 
improvements  and  cost  reductions.  Based  on  the  results  of  the  IMU 
calibration  tests  conducted,  it  is  concluded  that  the  developed  technique 
falls  into  th.s  category.  Major  areas  of  potential  usefulness  to  the 
services  are  thought  to  be  in  the  areas  of  quality  assurance  and  IMU 
cost-of-ownership  reductions. 

Increased  quality  assurance  could  be  achieved  by  i ncorporati ng  the 
developed  Vibration  technique  into  the  services  IMU  acceptance  tests. 
The  technique  is  well  suited  to  acceptance  testing  because  it  requires 
a minimal  amount  of  age  and  test  time  and  is  accurate  and  comprehensive. 
Pure  inertial  navigational  accuracy  projections  could  be  made  by 
covariance  analysis  simulations  using  the  results  of  the  IMU  error  model 
development  effort  and  calibration  test  results. 

As  a result  of  contract  F33615-73-C-1 146  with  AFAL , The  Analytic 
Sciences  Corporation  have  identified  the  need  for  improved  diagnostic 
methods  and  more  comprehensive  IMU  tests  in  the  field  of  service 
maintenance.  The  performance,  test  and  repair  history  of  the 
A-7D  AN/ASN~?0(V)  IMS,  using  the  KT-73  IMU,  was  studied  and  the  results 
are  highlighted  in  Reference  5.  The  problems  are  thought  to  be  an 


168 


AFAL-TR- 77-75 


inadequate  intermediate  level  maintenance  ability  to  identify  out-of-spec 
tMUs,  inconsistencies  between  and  the  failure  of  depot  level  tests  to 
verify  the  results  obtained  at  the  in  termed i ate  level,  and  iterative 
depot  testing  resulting  in  c minimum  of  thirty  plus  hours  test  time  on 
an  I MU  to  gain  confidence  that  it  is  operating  properly.  AGE  acquisition 
and  ownership  costs  also  contribute  significantly  to  I MU  cost-of-ownershi p . 
It  is  anticipated  that  adoption  of  the  developed  calibration  technique 
at  both  the  intermediate  and  depot  level  will  ameliorate  these  problems 
while  resulting  in  significantly  lower  I MU  ownership  costs. 

Based  on  the  above,  it  is  recommended  that  a test  program  be 
undertaken  at  the  Aerospace  Guidance  and  Metrology  Center  (AGMC), 

Heath  AFB  to  demonstrate  the  advantages  of  the  large  dimensional  Kalman 
Filter  calibration  technique  as  opposed  to  the  contemporary  approaches 
presently  being  employed.  A study  portion  of  the  program  would  attempt 
to  quantify  the  advantages  of  employing  the  technique  for  IMU  acceptance 
tests  and  at  the  intermediate  maintenance  level. 
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APPENDIX  A:  COMPUTER  PROGRAM  LISTING 
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APPENDIX  A 

IMU  CALIBRATION  CDC-6600  FORTRAN  IV  COMPUTER  PROGRAM  LISTING 

This  appendix  contains  a complete  listing  of  the  computer  program 
used  to  calibrate  the  KT-73  IMU  off-line  by  using  recorded  measurement 
data.  The  program  is  entered  on  a CDC-6600  computer  and  requires 
900  seconds  and  221,000  octal  words  memory  for  execution.  A "production" 
version  of  this  computer  program  would  reduce  both  the  execution  time  and 
memory  requirements  to  a small  fraction  of  those  above  by  eliminating 
the  covariance  matrix  computations  and  associated  large  dimensional 
matrices.  This  would  be  accomplished  by  storing  the  time-dependent 
Kalman  gain  matrices  on  a magnetic  tape. 

The  objective  of  depicting  this  listing  is  to  provide  a reference 
for  those  attempting  to  develop  a Kalman  Filter  algorithm  for  a similar 
application.  There  is  a multitude  of  data,  algorithms,  and  logic 
associated  with  the  initialization,  resetting,  and  execution  of  the  IMU 
calibration  Kalman  Filter  that  could  not  be  adequately  portrayed  in  the 
text  of  the  report.  This  listing  provides  an  accurate  definition  and 
description  of  these  details. 

IMU  calibration  results  are  shown  at  the  end  of  the  program  listing. 
These  results  are  the  filter's  estimate  of  the  error  sources  and  its 
predicted  estimation  accuracy  immediately  before  and  after  a measurement 
is  used  to  update  the  filter.  Data  is  showr  for  the  start,  of  the  test, 
at  the  end  of  each  platform  test  sequence,  and  for  the  instances  in  the 
third  and  tenth  test  sequences,  where  the  accelerometer  modes  are 
switched.  Units  for  the  output  data  are  consistent  with  the  estimation 
results  figures  given  in  Section  IX  in  degrees  per  hour,  arc  minutes, 
micro  g,  feet  per  second,  and  numeric  values. 
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